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ABSTRACT 


Weather  radar  data  is  used  to  extrapolate  hail  fall 
energy  to  locations  where  the  energy  was  not  measured  by 
hailpads  A  formula  to  convert  radar  echo  intensities  to 
energy  fluxes  of  hail  was  developed.  These  fluxes  were 
summed  over  the  duration  of  the  storm  at  a  given  location 
to  obtain  a  raw  estimate  of  energy  density  of  hailfall  at 
the  ground  in  Jm-2.  This  areal  energy  density  was  then 
adjusted  by  the  hailpad  network  using  a  linear  regression 
of  the  energy  densities  known  from  hailpad  analyses  against 
the  raw  radar  estimates  from  the  same  locations. 

Radar  data  from  August  7  and  18,  1974,  produced 
by  the  Alberta  Hail  Project  C-band  system,  was  used  in 
this  study.  This  data  is  logged  on  magnetic  tape,  and  was 
extracted  by  computer  programs  developed  during  the  study. 
Hailpad  data  from  August  7  and  18  was  used  with  the  radar 
data.  Supplementary  hailpad  information  from  July  30  and 
August  15,  1974  was  also  used. 

The  radar -hailpad  method  of  estimating  hail  energy 
density  had  a  spatial  resolution  equal  to  that  of  the  radar 
itself,  which  was  about  1  kilometer.  The  density  of  infor¬ 
mation  that  can  be  obtained  by  this  method  is  much  higher 


v 


than  that  of  the  existing  hailpad  network.  The  error  in 
the  radar-hailpad  method  was  found  by  statistical  analysis 
to  be  substantially  less  than  the  error  in  methods  based  on 
hailpad  information  alone.  Farmers’  reports  of  hailfall 
were  also  used  in  the  analysis. 

The  storms  on  the  two  days  studied  were  different 
enough  to  produce  changes  in  the  parameters  used  in  adjust¬ 
ing  radar  estimates  with  hailpad  data.  This  implies  that 
radar  estimates  of  hailfall  should  be  adjusted  on  a  day-to- 
day  basis  using  ground-based  measurements. 

A  continuous  grid  of  estimates  of  hailfall  energy 
density  was  calculated  by  the  radar-hailpad  method  for  a 
109  square  mile  (282km2)  area  hit  by  the  storm  of  August  18, 
1974.  The  details  of  the  resulting  pattern  agreed  strongly 
with  the  known  properties  of  the  distribution  of  hailfall 
on  that  scale.  The  estimated  energy  densities  also  agreed 
with  the  observations  of  hailfall  by  farmers  in  the  area. 


ACKNOWLEDGEMENTS 


I  would  like  to  acknowledg 

e  my  sincere 

gra 

t itude 

to  Dr 

.  Robert  Charlton  in  help 

ing  to 

guide  this  p 

ro  j 

ect  to 

compl 

etion,  in  spite  of  the 

se 

tbacks 

encountered 

over  the 

last 

three  years . 

I  would  also  like 

to 

expre 

ss  thanks  to 

Dr. 

Brian 

Barge 

and  Dr.  Robert  Humphr 

ies 

for  c 

tearing  up  pr 

obi 

ems  I 

had  i 

n  the  theory  of  weathe 

r  r 

adar , 

and  to  the  other 

member 

of  the  Alberta  Hail  Project 

re 

search 

staff  for  th 

eir 

many 

useful  suggestions.  Dr.  Carlt 

on  Ulb 

rich  of  Clems 

on 

Uni- 

ver  s  i 

ty  in  South  Carolina  s 

ent 

some 

interesting  c 

orr 

espon- 

dence 

on  the  results  of  his 

wo 

rk ,  wh 

ich  was  close 

iy 

related 

to  th 

is  project. 

Thanks  also  go  to 

Dr 

.  Guy 

Goyer  and  Dr. 

D. 

L. 

MacLe i 

sh  for  serving  on  my  exam 

ining 

committee . 

I  would  also  like 

to 

expre 

ss  my  appreci 

at  i 

on  to 

Lynne 

Rowen,  who  typed  the  f 

ina 

1  draft  of  this  the 

sis 

This  study  was  amply  funded  by  a  grant  from  the 
Interim  Weather  Modification  Board  under  the  Alberta  Minis¬ 
try  of  Agriculture. 


TABLE  OF  CONTENTS 


CHAPTER  Page 

1  INTRODUCTION  .  1 

1 . 1  The  problem .  1 

1.2  Historical  Background  of 

the  Study .  4 

1.3  Description  and  Objectives  .  7 

2  SOURCES  AND  PREPARATION  OF  DATA .  9 

2.1  Theory  of  the  Hailpad .  9 

2.2  Analysis  of  the  Hailpads .  11 

2 . 3  Radar .  16 

3  THEORETICAL  BACKGROUND  .  25 

3.1  The  Radar  Reflectivity- 

Energy  Density  Relationship  .  25 

3.2  The  Effect  of  Melting  on 

Kinetic  Energy  Flux .  31 

4  METHOD  OF  RADAR- HAILPAD  INTERPOLATION 

AND  PRELIMINARY  RESULTS  .  42 

4.1  Development  of  Operational 

Formulae .  42 

4.2  Preliminary  Results:  The 

Effect  of  Varying  Quantities 
of  Information  on  the  Delinea¬ 
tion  of  Hail  Areas .  59 

5  EVALUATION  OF  THE  RADAR-HAILPAD 
ESTIMATION  SCHEME  AND  ALTERNATIVE 

METHODS .  71 

5 .  1  Triplets .  71 

5.2  Description  of  the  Estimation 

Methods .  72 

vi  i  i 


CHAPTER  Table  of  Contents,  Continued  Page 

5  5.3  Evaluation  and  Comparison .  74 

5.4  Another  Statistical  Test .  75 

6  FACTORS  AFFECTING  THE  RELATIONSHIP 
BETWEEN  THE  HAILPAD  ENERGY  DENSITY 

AND  THE  RAW  RADAR  ESTIMATE .  81 

6.1  Attenuation .  81 

6.2  Non-uniform  Beam-filling  .  88 

6.3  Error  in  Registration  of  Radar  ....  89 

6.4  The  Problem  of  Resolution- 

Scale  Mismatch .  95 

7  CONCLUSION .  97 

7.1  Summary .  9  7 

7.2  Recommendations  for  Further  Study  ...  99 

LIST  OF  REFERENCES .  102 

APPENDIX  A  SUMMARY  OF  HAILPAD  DATA  USED  IN 

THIS  STUDY .  105 

APPENDIX  B  LISTINGS  OF  MAJOR  PROGRAMS  .  Ill 

APPENDIX  C  HAILPAD  TRIPLET  DATA  .  138 

APPENDIX  D  CENTER-PAD  ESTIMATES  AND 

ERRORS  OF  TRIPLETS  . .  142 

APPENDIX  E  ANALYSES  OF  RADAR  REGISTRATION 

ERRORS .  145 

APPENDIX  F  LAND  LOCATION  NOTATION  .  150 


IX 


LIST  OF  TABLES 

Table  Page 

4.1  Linear  correlation  coefficient,  r,  of 
radar  estimate  and  hailpad  value  of  energy 
density  for  August  7  and  18,  1974  at 
radar  calculation  thresholds  of  45, 

40  ,  35  and  30  dBZ .  57 

5.1  Means  and  standard  deviations  for  the 

errors  listed  in  Appendix  D,  in  Jm“ 2 .  76 

5.2  Results  of  test  on  size  of  mean  error 

in  estimating  center-pad  energy  density, 
described  in  Section  5.4 .  79 

6.1  Two-way  attenuation  coefficients,  in 
dBkm" 1 ,  for  some  hailpads  at  C-band 

and  S-band  radar  wavelengths  .  86 

6.2  Summary  of  radar  registration  errors, 

xT0TAL .  94 


x 


LIST  OF  FIGURES 


Page 


Figure 

2.1  Hailpad  Analysis  .  13 

2.2  Output  from  PLOT .  19 

2.3  Output  from  PLOTPAD .  21 

2.4  Output  from  PLOTPAD  using  revised 

radar  data .  24 

3.1  Effect  of  melting  on  kinetic 

energy  flux  of  hailfall .  38 

3.2  Effect  of  melting  on  kinetic 

energy  density  of  hailfall  .  40 


4.1  Sample  of  PLOTPAD,  illustrating  use  of 
Equation  4.6,  based  on  (a)  original 

radar  data,  and  (b)  revised  radar  data  ...  47 

4.2  Radar  estimates  versus  hailpad  values  of 


energy  density  at  hailpad  locations  for 
(a)  August  7,  1974  and  (b)  August  18, 

1974  .  52 

4.3  Hailfall  maps,  based  on  hailpad  analyses 
only,  for  (a)  August  7,  1974  and  (b) 

August  18  ,  1974  .  60 


4.4  Hailfall  maps,  using  hailpad  data  sup¬ 
plemented  by  farmers'  reports  for  (a) 

August  7,  1974  and  (b)  August  18,  1974  ...  63 

4.5  Detail  of  August  18,  1974,  showing  hail¬ 

pad  data  supplemented  by  farmers'  re¬ 
ports  .  67 

4.6  Detail  of  August  18,  1974  showing 
hailfall  energy  density  determined 

by  hailpad- adj usted  radar  estimates  ....  68 


5.1  Hailpad  triplet .  73 

6.1  Computer- generated  PPI  data .  91 


List  of  Figures,  Continued 

Figure  Page 

6.2  Data  of  Figure  6.1  as  it  would  appear 

in  an  actual  PPI .  92 

A-l  Land  location  diagram . 151 


XI  i 


CHAPTER  1 


INTRODUCTION 


1.1  THE  PROBLEM 

The  Alberta  Hail  Project  is  involved  in  cloud 
seeding  to  reduce  crop  damage  by  hail  over  a  20,000  km2 
area  of  central  Alberta.  The  spatial  pattern  of  hailfall 
is  very  intricate  and  the  Alberta  Hail  Project  would  like 
to  detect  hail  damage  reductions  in  the  order  of  ten  per¬ 
cent.  Therefore,  it  is  desirable  to  document  the  fine- 
scale  features  of  hailfall  for  individual  storms  to  the 
greatest  extent  possible  using  all  available  information. 

The  damage  to  grain  crops  which  is  caused  by 
a  hail  storm  has  been  shown  to  be  closely  related  to  the 
kinetic  energy  per  unit  area  of  the  hailstones  which  fell 
during  the  storm.  For  instance,  100  joules  per  square 
meter  means  slight  crop  damage  and  1000  joules  per  square 
meter  usually  means  100  percent  crop  damage.  Although 
some  of  a  falling  stone's  kinetic  energy  is  tied  up  in 
its  horizontal  travel  with  horizontal  winds,  most  of  its 
energy  is  usually  incorporated  in  its  vertical  fall  speed. 
In  this  study,  the  sum  of  the  vertical  kinetic  energies  of 
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the  hailstones  which  fell  on  a  unit  area  during  a  storm 
will  be  called  the  "kinetic  energy  density"  although  it 
is  really  a  kinetic  energy  per  unit  area. 

This  study  investigates  the  use  of  weather  radar 
to  estimate  hailfall  from  individual  storms  between  points 
where  hailfall  information  is  already  available.  Since 
the  resolution  of  the  radar  is  about  one  kilometer  it  is 
theoretically  possible  to  estimate  the  hailfall  from  a 
storm  over  every  square  kilometer  of  land. 

The  Alberta  Hail  Project  relies  on  several  sour¬ 
ces  to  obtain  data  on  hail  activity:  hailpads,  radar, 
ground  surveys  by  mail  or  telephone,  field  trips  in 
specially- ins trumented  trucks,  and  crop  insurance  statis¬ 
tics. 

Telephone  and  mail  surveys,  called  farmers’ 
reports,  cover  a  wide  area  and  can  yield  a  high  areal 
density  of  hail  information.  This  information  is  rather 
imprecise,  since  observing  conditions  during  a  hailstorm 
are  hardly  ideal  for  the  amateur  observer  who  is  expected 
to  produce  quantitative  information.  Field  trips  by  pro¬ 
fessional  observers  are  expensive  and  time-consuming.  In 
the  history  of  a  typical  hailstorm,  all  the  field  trips 
undertaken  will  cover  only  about  twenty  locations.  While 
crop  insurance  statistics  reveal  a  great  deal  about  total 
storm  damage  and  seasonal  losses,  they  have  insufficient 
resolution  to  delineate  the  fine-scale  storm  effects. 

They  can,  however,  be  used  in  conjunction  with  high- 


, 
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density,  area- averaged  information  for  general  assessments 
of  storm  damage. 

Until  1974,  the  Alberta  Hail  Project  has  used 
its  radar  facilities  mainly  for  hailstorm  research.  In 
July  1974  a  PDP-11  computer  was  added  to  the  radar  sys¬ 
tem.  The  computer  makes  a  digital  record  on  magnetic 
tape  of  all  radar  echoes  on  hailstorm  days  as  well  as 
providing  instantaneous  storm- tracking  information. 

A  hailpad  is  a  piece  of  Styrofoam  about  30  cen¬ 
timeters  square  and  3  centimeters  thick.  In  the  field,  it 
is  mounted  on  a  firm,  level  surface.  Each  hailstone  im¬ 
pacting  the  Styrofoam  will  dent  it,  providing  a  record  of 
the  storm’s  passage  over  the  pad.  The  principles  of  hail¬ 
pad  analysis  are  described  in  Chapters  2  and  3.  In  the 
1974  operations  seasons,  about  600  hailpad  stations  were 
maintained  in  the  "Project"  area.  The  regular  network  con¬ 
sisted  of  491  hailpads.  The  hailpad  spacing  varied  from 
three  per  township  (36  square  miles  or  93.2  square  kilo¬ 
meters)  near  the  central  region,  to  about  one  per  town¬ 
ship  on  the  fringes.  Farmers  who  maintained  hailpad  sta¬ 
tions  also  kept  raingauges  and  collected  hail  samples. 

For  finer-scale  studies  the  Project  also  maintained  a 
dense  network  of  about  100  hailpads  over  ten  townships 
between  the  towns  of  Rimbey  and  Eckville. 

In  1973  Strong  (1975)  studied  the  use  of  hail¬ 
pads  in  providing  quantitative  information.  He  found 
that  this  information  would  be  useful  in  studying 
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entire  storms  only  if  the  density  of  the  hailpad  network 
was  made  much  higher  than  it  has  been  to  date.  His  con¬ 
clusions  are  based  on  the  use  of  hailpads  as  the  only 
source  of  information  on  spatial  distribution  of  hailfall 
energy.  The  suggested  number  of  hailpads  required  to  suc¬ 
cessfully  capture  these  variations  over  the  entire  Pro¬ 
ject  area  was  clearly  impracticable.  A  good  alternative 
is  to  supplement  the  data  from  the  existing  hailpad  net¬ 
work  with  radar  information.  The  use  of  radar  and  hail¬ 
pads  instead  of  hailpads  alone  has  two  major  advantages. 
First,  the  radar  can  cover  the  entire  area  under  observa¬ 
tion.  The  second  advantage  is  its  comparatively  low  cost. 
The  required  information  can  be  abstracted  from  data  col¬ 
lected  during  regular  radar  operation  which  involves  de¬ 
tection  of  storms  and  the  guidance  of  seeding  aircraft. 

1.2  HISTORICAL  BACKGROUND  OF  THE  STUDY 

Since  weather  radars  were  used  to  study  rainfall 
rate  before  they  were  used  to  study  hailfall,  it  is  best 
to  review  the  historical  developments  in  that  order. 

Early  in  the  history  of  radar  meteorology,  it 
was  realized  that  the  intensity  of  radar  echoes  from  rain 
could  be  correlated  with  the  rainfall  rate.  An  empirical 
relationship  between  the  radar  reflectivity  Z,  and  the 
rate  of  rainfall  R,  was  suggested  by  Marshall,  Langille 
and  Palmer  as  early  as  1947.  This  and  similar  formulae 
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have  since  been  referred  to  as  Z-R  relationships.  Using 
observations  of  raindrop- s ize  distributions,  Marshall  and 
Palmer  (1948)  developed  the  following  theoretical  Z-R 
relationship : 


Z  =  2  2  OR 1 • 6  1.1 

with  Z  in  mm6m" 3  and  R  in  mm  hr'1. 

In  the  1960's  Soviet  scientists  used  raingauges 
to  verify  Z-R  relationships.  Some  of  these  studies  were 
reported  by  Dimaksian,  Zotimov  and  Zykov  (Dimaksian  et  al .  3 
1965) . 

From  1963  onward,  studies  of  the  Z-R  relationship 
proliferated  rapidly.  Wilson  (1966)  showed  that  storm- 
to-storm  variability  exists  in  the  relationship.  Stout 
and  Mueller  (1968)  discovered  significant  regional  as 
well  as  storm- to- storm  variation.  It  became  clear  that 
a  ground  network  of  precipitation  gauges  was  needed 
to  calibrate  the  radar  for  each  storm  if  it  was  to  suc¬ 
cessfully  estimate  rainfall  rate  over  the  area  of  radar 
coverage.  There  is  an  implication  here  that  hailfall- 
radar  studies  may  also  need  ground-based  gauges. 

Raingauges  were  actively  used  in  calibrating  or 
adjusting  radar  estimates  of  rainfall  in  a  major  study  in 
Illinois  in  1964  and  1965.  Estimates  of  the  total  volume 
of  fallen  rain  from  the  existing  "climatological"  raingauge 
network  were  compared  with  estimates  from  radar  data  cali- 
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brated  by  the  raingauges.  The  known  statistical  properties 
of  the  raingauge  network  were  used  to  judge  the  relative 
accuracy  of  the  two  methods.  The  study  showed  that  the 
radar,  calibrated  by  raingauges,  was  equivalent  in  accur¬ 
acy  to  a  raingauge-only  network  of  1 5 Omi 2 /gauge .  The  cli¬ 
matological  raingauge  network  had  an  actual  density  of 
225mi2/gauge  (Huff,  1966). 

Newer  techniques  of  integrating  raingauge  and 
radar  data  improved  the  accuracy  of  areal  rainfall  esti¬ 
mates  in  the  Florida  Area  Cumulus  Experiment  (Woodley  et 
al . }  1974).  Other  studies  have  consistently  demonstrated 
the  value  of  combined  radar-ground  gauge  methods 
(Changnon,  1967;  Wilson,  1970;  Brandes,  1974,  1975). 

The  possibility  that  hailfall  could  be  analysed 
in  a  similar  way  to  rainfall  depended  on  the  development 
of  a  simple  and  reliable  instrument  for  measuring  the 
important  parameters  associated  with  hail  at  the  ground. 
These  parameters  are  mass,  kinetic  energy  or  momentum  per 
unit  area.  The  hailpad  has  become  the  favorite  instru¬ 
ment  for  this  purpose.  It  was  first  used  by  Schleusener 
and  Jennings  (1960).  They  used  the  3 0  -  cent imeter  size, 
covered  with  aluminum  foil,  which  was  in  common  use  until 
1975.  Hailpad  usage  was  first  described  by  Decker  and 
Calvin  (1961).  The  Alberta  Hail  Project  used  this  type 
of  hailpad  in  1974,  but  it  abandoned  the  use  of  aluminum 
foil  in  the  following  year.  Quantitative  methods  of  ana¬ 
lysing  foil-covered  and  "styrofoam  only"  hailpads  were 
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developed  in  Alberta  (Strong,  1974;  Strong  and  Lozowski, 
1977)  . 

The  theory  of  the  radar  detection  and  evaluation 
of  hailfall  has  been  studied  by  many  researchers.  Prac¬ 
tical  applications  were  recognized  as  early  as  1960 
(Douglas) .  A  comprehensive  review  can  be  found  in  "Radar 
Meteorology"  (Battan,  1973).  It  has  been  shown  that  the 
intensity  of  hailfall  can  be  related  to  the  strength  and 
other  characteristics  of  the  radar  echoes  in  much  the  same 
way  as  rainfall  (Barge,  1977).  It  is  possible,  therefore, 
that  hailpads  and  radar  can  be  used  to  delineate  the  pattern 
of  hailfall  in  an  analogous  way  to  the  use  of  raingauges 
and  radar  for  rainfall.  The  problem  is  somewhat  more  dif¬ 
ficult  than  rainfall  evaluation  because  of  the  generally 
smaller  dimensions  of  hailshafts  and  because  echoes  from 
hail  are  generally  masked  to  some  extent  by  those  from 
rain . 

1.3  DESCRIPTION  AND  OBJECTIVES 

The  main  objective  of  the  study  was  to  develop 
a  method  of  estimating  the  kinetic  energy  density  of  hail¬ 
fall  at  locations  where  only  radar  data  was  available.  The 
digitized  output  of  meteorological  radar  and  hailpad  in¬ 
formation  were  used.  Such  an  analysis  of  hailfall  could 
probably  be  related  to  statistics  of  crop  damage  by  hail, 
and  it  could,  therefore,  aid  in  the  assessment  of  the  ef- 
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fectiveness  of  cloud  seeding. 

This  study  uses  radar  and  hailpad  information 
from  two  hail-days  in  central  Alberta  in  1974.  August  7 
and  August  18,  1974  yielded  a  large  number  of  useful  hail- 
pads,  while  the  quality  of  the  radar  data  was  better  than 
any  available  on  other  hail-days  during  the  season.  Sup¬ 
plementary  hailpad  information  from  July  30  and  August  15 
was  also  used.  The  radar  was  the  Alberta  Hail  Project 
C-band,  operating  at  a  wavelength  of  5.5  centimeters. 

The  S-band,  10-centimeter  unit,  which  would  normally  be 
used,  was  under  repair  throughout  the  summer  of  1974. 

When  the  data  for  this  study  was  being  collec¬ 
ted,  the  method  of  digitizing  the  radar  output  was  new, 
and  defects  in  the  data  extraction  were  still  in  the  pro¬ 
cess  of  being  removed.  As  a  result,  there  are  weaknesses 
and  interruptions  in  the  radar  data  even  under  the  best 
conditions,  and  the  overall  quality  is  poor  in  comparison 
to  the  data  from  later  years.  This  dissertation  describes 
the  analysis  of  this  data  and  the  development  of  a  method 
of  radar-hailpad  interpolation.  The  results  of  this 
method  are  shown  in  Chapter  4.  Their  accuracy,  in  compari¬ 
son  to  methods  using  hailpads  alone,  is  determined  by 
elementary  statistical  tests.  The  radar  and  hailpad  data 
from  1975  and  1976  are  expected  to  provide  a  more  thorough 
test  of  the  method  described  here. 


CHAPTER  2 


SOURCES  AND  PREPARATION  OF  DATA 


2.1  THEORY  OF  THE  HAILPAD 


This  section  summarizes  how  a  hailpad  measures 
the  energy  density  of  hailfall. 

Strong  and  Lozowski  (1977)  have  shown  that  a 
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D  is  the  diameter  of  the  spherical  hailstone  and  k  is  the 
constant  of  proportionality  which  depends  on  the  factors 
mentioned  above.  The  momentum  and  kinetic  energy  of  the 
hailstone  can  be  inferred  from  this  equation.  A  count  of 
the  dents  in  particular  size  categories  will,  by  applica¬ 
tion  of  the  proper  formulae,  yield  the  kinetic  energy  per 
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unit  area  of  the  hail  that  fell  on  the  pad. 

The  hailpads  used  by  the  Alberta  Hail  Project 
are  made  of  a  brand  of  Styrofoam  (Type  *FR)  manufactured 
by  the  Dow  Chemical  Company.  In  1974  they  were  covered 
with  heavy  (.001  inch)  aluminum  foil.  A  key  assumption 
made  by  Strong  and  Lozowski  is  that  the  hailpad  exerts  a 
constant  pressure  against  the  surface  of  a  solid,  homo¬ 
geneous  sphere  as  it  undergoes  collision.  Most  of  the 
kinetic  energy  of  the  hailstone  is  assumed  to  be  used  in 
creating  the  dent  in  the  styrofoam.  The  center  of  the 
stone  does  not  penetrate  below  the  undisturbed  surface  of 
the  hailpad,  and  the  sphere  does  not  penetrate  so  deeply 
that  the  support  underneath  the  pad  affects  the  dynamics 
of  the  collision.  Fracturing  of  the  Styrofoam  and  the  ef¬ 
fect  of  deforming  the  aluminum  foil  cover  are  assumed  to 
be  unimportant.  Using  these  assumptions,  the  authors  de¬ 
rived  an  equation 
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where  R  is  the  hailstone  radius,  r  is  the  radius  of  the 
dent  it  produces,  p  is  the  pressure  the  Styrofoam  exerts 
against  the  hailstone,  g  the  gravitational  acceleration, 
Cq  the  drag  coefficient,  and  pj  and  pa  are  the  densities 
of  ice  and  air,  respectively. 

Based  on  laboratory  experiments  the  authors  de¬ 
veloped  a  polynomial  to  express  the  hailstone-dent  rela- 
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2.4 


Summing  the  kinetic  energy  given  by  Equations  2.3  and  2.4 
from  all  dents  on  the  hailpad  and  dividing  by  the  area  of 
the  pad  gives  the  energy  per  unit  area,  or  energy  density, 
of  the  hail  that  fell  at  the  hailpad  site.  To  get  the 
hailstone’s  kinetic  energy  in  units  of  joules  the  coeffi¬ 
cient  of  D4  in  Equation  2.4  is  typically  0.042  if  D  is 
in  centimeters.  The  methods  actually  used  to  calculate 
energy  density  are  described  in  the  following  section. 


2.2  ANALYSIS  OF  THE  HAILPADS 

The  original  analysis  of  the  hailpads  from  the 
1974  Hail  Project  operations  was  performed  by  part-time 
helpers  and  the  author  of  this  study  throughout  1975  and 
1976,  using  a  method  developed  by  Wojtiw  (1975). 

Each  analyst  was  supplied  with  sheets  of  clear 
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acetate  plastic  to  cover  the  exposed  faces  of  the  hail- 
pads  to  be  analysed,  and  a  transparent  template  marked 
with  circles  of  various  diameters  ranging  from  0.1  inches 
to  2.25  inches.  The  borders  of  the  dents  on  the  hailpads 
are  traced  on  the  acetate  sheets.  This  prevents  damage  to 
the  hailpads  and  provides  a  convenient  means  to  mark  off 
dents  as  the  analysis  proceeds. 

Using  the  template,  a  count  is  made  of  the  num¬ 
ber  of  dents  in  diameter  classes  at  intervals  of  0.05 
inches.  Diameter  classes  not  represented  on  the  template 
are  counted  by  the  dents  whose  diameters  lie  between  the 
two  nearest  available  circle  sizes.  Large,  irregular  dents 
are  counted  by  using  the  circles  of  approximately  the  same 
area.  The  counts  of  the  dent  sizes  are  entered  in  the 
left-hand  blank  column  of  the  analysis  sheet,  an  example 
of  which  is  shown  in  Figure  2.1.  The  kinetic  energy  per 
dent  in  each  diameter  class  is  listed  beside  the  blank 
center  column.  The  counts  of  the  dents  are  multiplied  by 
the  corresponding  kinetic  energies  per  dent  and  the  result¬ 
ing  figures,  in  joules,  are  entered  in  the  center  column. 
This  column  is  totalled  to  give  the  kinetic  energy  den¬ 
sity  in  joules  per  square  foot,  then  multiplied  by  10.76 
to  give  the  result  in  joules  per  square  meter. 

Subsequent  publication  of  the  results  of  hail- 
pad  calibration  by  Lozowski  and  Strong  (1977)  led  to  the 
discovery  of  a  discrepancy  between  their  calibration  for¬ 
mulae  (Equations  2.3  and  2.4)  and  the  ones  implicit  in 
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Figure  2.1:  A  hailpad  analysis  sheet,  based  on  actual 

hailpad  data.  The  hailpad  is  H37,  from 
August  7,  1974.  Note  that  the  entries  under 
"Total  impact  energy--by  size"  are  in  joules. 
These  are  obtained  by  multiplying  the  "Impact 
energy"  factors  by  10-7. 
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the  calculations  of  energy  done  in  the  original  analyses. 
There  is  no  difference  in  the  basic  method  of  hailpad 
analysis.  However,  calculation  of  the  hail  energy  density 
from  the  dent  counts  by  Equations  2.3  and  2.4  led  to  a 
value  smaller  than  the  value  from  the  original  analysis 
by  approximately  a  factor  of  3.  This  new  calibration,  as 
described  by  Strong  and  Lozowski  (1977),  was  judged  to  be  more  reli¬ 
able  because  of  the  care  and  thoroughness  with  which  Strong  and  Lozowski 
established  their  calibration  formulae.  The  significance  of  this  find¬ 
ing  will  be  discussed  in  Chapter  4. 

The  error  from  the  theoretical  uncertainties  in 
the  procedure  is  about  10  percent.  The  variations  in  the 
kinetic  energy  density  for  the  same  hailpad  obtained  by 
different  analysts  indicate  that  an  additional  error  of  up 
to  15  percent  can  sometimes  be  produced  by  the  analysis 
method  (Strong,  1974).  The  energy  density  of  hailfall 
determined  by  a  hailpad  analysis  will,  however,  be  ac¬ 
cepted  as  the  "true”  value.  Hailpads  have  an  operational 
range  of  0  to  about  2,000jm'z.  Some  damage  to  grain  crops 
is  usually  detected  at  energy  densities  of  50jm'2  or  more. 
Crops  are  often  totally  lost  at  400jm“2  or  greater.  Above 
2,000jm“2,  the  hailpad  is  usually  destroyed  (Strong, 

1974).  Storms  of  this  magnitude  cause  massive  damage, 
not  only  to  crops,  but  to  cars  and  buildings  as  well. 

Thirty-one  hailpads  from  July  30,  1974;  49  from  August 
7;  7  from  August  15;  and  44  from  August  18  were  analysed.  Not 
all  hailpads  within  the  areas  of  these  storms  were  avail- 
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able.  Many  pads  were  missing  or  had  faulty  or  inadequate 
documentation,  so  that  the  mean  effective  spacing  was  just 
over  one  per  standard  township  (93km2).  The  pertinent 
information  on  all  the  hailpads  used  in  this  study,  includ¬ 
ing  the  results  of  both  the  original  and  later  calibrations, 
is  given  in  Appendix  A.  Dent  size  distributions  for  all 
these  hailpads  were  compiled,  but  they  are  not  tabulated 
here.  The  hailpad  locations  are  recorded  as  the  quarter- 
section  of  land  on  which  the  hailpad  was  set  up.  This  is  a 
0.25mi2  area  (about  0.65km2). 

2 . 3  RADAR 

The  C-band  radar  data  has  undergone  considerable 
processing  to  reduce  it  to  the  form  used  here.  The  first 
stage  was  an  analog- to-digital  conversion  which  transformed 
the  standard  PPI  (plan-position  indicator)  display  into  an 
array  of  "bins".  Each  bin  occupied  1.05  kilometers  of  range 
and  0.545  degrees  of  azimuth.  The  bins  covered  a  range  of 
from  3  to  147  kilometers.  The  returned  power  of  the  radar 
in  each  bin  was  converted  into  a  numerical  value  as  the 
radar  beam  scanned  elevations  from  0°  to  21.0°  at  1.5° 
intervals . 

The  second  stage  was  the  recording  of  the  bin 
values  on  magnetic  tape.  Computer  programs  handled  such 
aspects  as  placing  the  bin  values  on  tape  and  labelling 
them  for  azimuth,  range  and  time,  indicating  changes  in 
antenna  elevation  and  the  beginning  and  end  of  scan  cycles. 
Other  programs  handled  the  calibration  of  the  radar  after 
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analog- to-digital  conversion  of  the  received  power. 

The  magnetic  tapes  containing  the  radar  data 
were  copied  and  run  on  the  University  of  Alberta  computer, 
using  programs  designed  especially  for  this  study  by 
Fred  McDougall  and  Michel  Delorme.  The  first  program, 
called  PLOT,  was  activated  by  a  series  of  land  locations 
which  were  selected  by  the  operator.  The  land  location 
description  used  was  the  standard  code  of  quarter- sect  ion , 
section,  township,  range  and  meridian  placement  (see  Ap¬ 
pendix  F) .  PLOT  accepted  any  land  location  within  range 
of  the  radar,  whether  or  not  it  contained  a  hailpad.  The 
operator  also  supplied  the  slope  and  intercept  of  the  cali¬ 
bration  curve  and  the  radar  calibration  constant  for  the 
desired  date.  The  radar  calibration  constant  was  the  dif¬ 
ference,  in  dB,  between  the  transmitted  power  and  the  radar 
constant,  both  in  dBm.  The  calibration  curve  determined 
the  relation  between  the  returned  power  of  the  radar  echo 
and  the  digital  value  stored  on  tape.  The  radar  constant 
corrects  for  the  range  of  the  target. 

The  program  PLOT  selected  the  nearest  bin,  at 
all  elevations,  to  each  land  location.  It  then  surveyed 
the  entire  history  of  the  radar  echoes  represented  by  the 
bin  set  over  a  time  span  selected  by  the  operator,  after 
converting  the  raw  bin  values  to  equivalent  radar  reflec¬ 
tivity  in  decibels  (dBZ) .  These  data  were  recorded  on  a 
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second  tape.  The  program  also  calculated  a  general  hourly 
summary  of  the  echo  history  over  each  selected  land  loca¬ 
tion.  An  example  of  the  output  from  PLOT  appears  in  Figure 

2.2. 

The  second  program, PLOTPAD ,  used  the  tape 
generated  by  PLOT  to  produce  a  graphic  time-height  sum¬ 
mary  of  the  echoes  over  land  locations  selected  by  the 
operator.  The  correct  file  created  by  PLOT  had  to  be  spe¬ 
cified  in  order  to  use  the  desired  land  locations,  and 
only  those  locations  used  in  PLOT  could  be  run,  since  a 
single  run  of  PLOT  opened  and  closed  a  single  file  on  the 
output  tape.  The  result  was  effectively  a  radar  history 
of  the  storm  as  it  passed  over  the  land  location  in  ques¬ 
tion.  Indeed,  most  of  these  time-height  summaries  re¬ 
sembled  the  thunderstorm  cross-sections  that  appear  in 
meteorology  textbooks.  Figure  2.3  shows  an  example  of  the 
output  from  PLOTPAD. 

Complete  texts  of  PLOT  and  PLOTPAD  are  given  in 
Appendix  C.  The  output  of  these  programs  had  direct  ap¬ 
plications  in  establishing  the  timing  and  local  severity 
of  the  storms.  PLOTPAD,  in  particular,  was  able  to  show 
subtle  changes  in  the  structure  and  development  of  the 
storms  when  continuous  grids  of  land  locations  were  surveyed. 
The  tape  output  from  PLOT  was  applied  to  direct  calculations 
of  radar-based  estimates  of  hail  energy  density. 

Good  radar  data  was  available  throughout  the 
afternoon  and  evening  of  August  7,  1974  so  that  all  hail- 
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Figure  2.2:  A  portion  of  the  output  from  the  program  PLOT. 

The  input  data  described  on  page  17  appear 
at  the  head  of  the  listing.  The  list  of  times 
immediately  following  the  tape-mount  lines  are 
the  times  when  the  antenna  returns  to  the 
O-degree  elevation,  0-degree  azimuth  position. 
They  provide  a  check  on  the  continuity  of  the 
data.  Limitations  of  space  prevented  showing 
the  complete  data  summary,  which  would  include 
hailpad  G39.  In  the  complete  output  listing, 
sets  of  O-degree  time  markers  and  data  sum¬ 
maries  would  be  printed  at  hourly  intervals 
through  the  time  period  specified  by  the  opera¬ 
tor. 
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Figure  2.3:  The  output  from  PLOTPAD  for  hailpad  H37,  using 

the  original  radar  data.  The  radar  record  of& 
a  single,  complete  storm  is  shown.  The  an¬ 
tenna  elevation  is  listed  in  degrees  at  the 
extreme  left,  and  the  corresponding  height  in 
kilometers  at  the  extreme  right.  The  bottom 
line  in  each  frame  is  the  time,  in  minutes, 
from  the  first  non-zero  echo  appearing  on  the 
current  data  tape.  Several  data  tapes  may 
be  used  to  produce  a  PLOTPAD  summary. 
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pads  collected  from  that  day  had  complete  radar  histories. 
The  radar  data  from  August  18,  1974  was  good  until  1832 
(6:32  p.m.),  after  which  numerous  tape  and  software  break¬ 
downs  rendered  the  radar  output  unusable.  Thus,  com¬ 
plete  radar  histories  were  available  for  only  the  eighteen 
hailpads  dented  by  the  storm  had  finished  by  1832  on  August 
18.  Neither  July  30  nor  August  15  had  good  radar  data, 
but  were  selected  mainly  for  their  good  hailpad  data. 

The  radar  data  initially  was  extracted  by  PLOT 
using  a  radar  calibration  constant  of  68.4dB.  More  re¬ 
cently  Humphries  (1978)  pointed  out  that  this  figure  did 
not  take  into  account  line  losses  in  the  waveguide  and  a 
discrepancy  between  the  way  a  radar  system  with  a  linear 
response  to  returned  power  and  one  with  a  logarithmic  res¬ 
ponse  (the  C-band)  processes  reflectivity  data.  The  com¬ 
bination  of  these  effects  adds  4.8dB  to  the  radar  cali¬ 
bration  constant.  Thus  all  the  original  radar  data  was 
re-analysed  using  a  radar  calibration  constant  of  73.2dB. 

The  significance  of  this,  along  with  the  modifications  in¬ 
troduced  to  the  hailpad  data,  will  be  examined  in  the  fourth 
chapter.  Figure  2.4  illustrates  the  output  from  PLOTPAD 
over  the  same  hailpad  as  Figure  2.3,  except  using  the  im¬ 
proved  radar  data. 
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Figure  2.4:  Same  as  figure  2.3,  except  using  the  revised 

radar  data.  Note  that  the  echoes  are  about 
5dB  higher  than  in  the  original  radar  data. 


CHAPTER  3 


THEORETICAL  BACKGROUND 


3.1  THE  RADAR  REFLECT  I VITY- ENERGY  DENSITY  RELATIONSHIP 


In  this  section  a  relationship  between  the 
radar  echo  and  the  downward  flux  of  hailfall  kinetic  energy 
is  developed.  The  t ime - integrated  flux  of  kinetic  energy 
determines  the  total  hailstone  energy  intercepted  by  a 
hailpad  during  a  hailstorm.  The  influence  of  hailstone 
melting  between  the  echo  level  and  the  ground  is  taken  into 
account.  Following  evaluation  in  Chapters  4  and  5  prob¬ 
lems  associated  with  this  method  are  discussed  in  Chapter  6 

By  using  some  simplifying  assumptions,  a  formula 
can  be  developed  that  converts  the  reflectivity  factor  of 
a  radar  echo  (in  dBZ)  to  a  quantity  expressing  the  down¬ 
ward  energy  flux  density  of  hail  (in  joules  per  square 
meter  per  second)  illuminated  by  the  radar  beam.  A  deri¬ 
vation  of  this  formula  based  on  studies  by  Barge  (1977)  , 
Smith  et  al .  (1975),  and  Ulbrich  (1977),  follows. 

The  kinetic  energy  flux  (E)  of  spherical  hail¬ 
stones  passing  through  a  unit  horizontal  area  can  be  ex¬ 
pressed  as 
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E  =  E  z^v2  •  v  3.1 

volume 

The  form  of  Equation  3.1  is  based  on  the  kinetic  energy 
formula  for  each  particle  having  mass  m  and  fall  velocity 
v.  "Volume"  refers  to  the  unit  volume  occupied  by  the 
hailstones.  The  kinetic  energy  flux  E  is  therefore  equal 
to  the  product  of  the  kinetic  energy  per  hailstone  and  its 
fall  velocity,  summed  over  all  hailstones  in  the  unit 
volume . 

Assuming  that  there  are  no  updrafts  or  downdrafts 
hailstones  will  fall  at  their  terminal  velocity  inside 
a  storm.  Numerous  studies  have  shown  that  the  terminal 
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by  substituting  Equations  3.2  and  3.3  for  m  and  v. 


In 


27 


practice,  Equation  3.4  can  be  replaced  by  an  integration: 

E  =  -fjPiitk3  £==0DMAXn(D)D4'5dD  3.5 

where  n(D)dD  is  the  number  of  hailstones  per  unit  volume 
in  a  diameter  interval  (D,D+dD).  D^x  is  the  maximum 
hailstone  diameter  in  the  unit  volume. 

Justification  for  approximating  n(D)  in  the  form 
of  a  Marshal 1 - Palmer  exponential  is  given  in  a  study  by 
Smith  et  al.  (1975),  in  which  case: 

n(D)  =  N0exp(-AD)  3.6 

N0  and  A  are  distribution  parameters.  N0  is  assumed  to 
be  constant  for  all  hailfalls  and  A  is  assumed  to  depend 
on  hailfall  rate  in  equivalent  depth  of  liquid  water  per 
unit  time.  In  earlier  studies  by  Douglas  (1963,  1964),  a 
similar  exponential  distribution  was  postulated,  except 
that  A  was  held  fixed  and  N0  left  to  vary. 

Using  a  Marshal 1 - Palmer  distribution  allows  the 
assumption  that  E^AX  is  very  large,  so  that  Equation  3.5 
takes  the  form 

E  =  -jyP|TI k3No  jj^exp  (- AD) D" ••5  dD  3.7 


The  integral  is  easily  evaluated: 
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E  =  y^-p  i7Tk3N0r  (5.  5)  A” 5  * 5  3.8 

where  r(x)  is  the  complete  gamma  function.  For  instance, 

T  (5 .  5)  has  a  value  of  52.3. 

A  relationship  between  A,  N0  and  Ze  is  now  re¬ 
quired.  Z e,  the  effective  or  equivalent  reflectivity  fac¬ 
tor,  is  the  reflectivity  factor  of  an  ensemble  of  Rayleigh- 
scattering  particles  that  would  give  the  same  returned 
radar  beam  power  as  the  precipitation  particles  under  ob¬ 
servation.  The  Rayleigh  reflectivity  factor  is  given  by 

Z  =  /“0n(D)D6dD 

Hence,  for  a  Marshall-Palmer  distribution  of  particles, 

Z  =/°°No  exp  (-AD)D6dD 
o 

Integrating  this  expression  yields 

Z  =  7  2  ON0  A" 7  3.9 

If  N0  =  3xl0~4cm'4  (Smith  et  al .  3  1975),  A  has  units  of 
cm'1  and  Z  is  in  cm3.  The  question  of  converting  N0  into 
SI  units  and  Z  into  standard  units  of  mm6m" 3  will  be  dealt 
with  shortly. 

Smith  et  al.  (1975)  adjust  the  right-hand  side 
of  Equation  3.9  by  a  power  of  0.95  to  obtain  the  effective 
reflectivity  factor,  Ze, 


at  a  10.6cm  radar  wavelength. 
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A  calculation,  similar  to  that  of  Smith  et  at.  (1975) 
showed  that  at  the  5.5cm  wavelength,  the  corresponding 
power  ranged  from  0.97  to  0.93  for  all  hailstone  diameters 
less  than  3cm.  Thus,  the  relationship  of  Smith  et  al . 
(1975)  , 


Ze  =  (720N0A"7)0'95 


3.9a 


was  sufficiently  accurate  for  the  purpose  of  this  study. 

Equation  3.8  for  energy  flux  can  now  be  written 
in  terms  of  Ze  by  using  Equation  3.9a  to  eliminate  A,  the 
hailfall  rate  parameter: 


E  =  ~Pik3r(5.5)No 


_ - _  #  n  0  •  8  3 

(72 ON 0 )  ^  6 


3.10 


To  evaluate  k  (in  Equation  3.2)  one  assumes  that 
the  hailstone  is  falling  at  terminal  velocity,  so  that  its 
weight  equals  its  drag.  Therefore, 


V 


t 


A  Pa  C[y 


2  1 

V1 


3 .11 


3.11a 


where  pa  is  the  ambient  air  density,  g  is  the  acceleration 
due  to  gravity,  and  Cp  is  the  drag  coefficient.  The  value 
of  Cp  for  a  smooth  sphere  with  a  diameter  less  than  11cm 


V 
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is  approximately  0.45.  For  hailstones  of  average  rough¬ 
ness  a  more  realistic  value  of  0.5  is  used.  The  freezing 
level  in  most  Alberta  hailstorms  is  near  750mb,  corres¬ 
ponding  to  an  air  density  of  0 . 95x10" 3kgm" 3 .  Using 

these  values  and  the  density  of  ice  as  previously  given, 

i 

k  is  found  to  be  159m2sec_1.  The  freezing  level  is  a  good 
typical  level  for  radar  observations  used  in  this  study. 
In  Alberta  a  typical  ground  level  k  based  on  a  pressure 

i 

of  900mb  and  a  temperature  of  20°C  would  be  133m2sec_1. 
This  was  confirmed  by  photographic  studies  by  Beattie 
(1975) . 

When  N0  is  expressed  as  3xl04m-4,  the  constants 
and  parameters  in  Equation  3.10  can  be  evaluated,  so  that 

E  =  2 . 625x10  9Ze  °’  8  3  3  .  1  2 


where  E  is  in  units  of  Jm~2sec"1.  As  Equation  3.12  is 
written,  Ze  has  units  of  m3,  or  m6m" 3 .  To  adjust  Equation 
3.12  so  that  Ze  has  the  usual  units  of  mm6m"3,  the  numeri¬ 
cal  coefficient  is  multiplied  by  a  factor  of  1.30xl0“15, 
leaving 


E  =  3.40xlO"6Ze0,83 


E  =  3. 40x10“6x 1012'  1 


5.13 


3.14 


e  > 


where  the  observed  effective  reflectivity  factor,  Z 


is 
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in  dBZ  (Ze[dBZ]  =10  logZeCmm6m" 3 ] ) .  Ze  will  be  given 
in  units  of  dBZ  from  here  onwards,  unless  specified  other¬ 
wise.  Equation  3.14  is  valid  in  a  hail  cloud  but  it  must 
be  adjusted  for  melting  and  air  density  changes  if  E  is 
to  be  integrated  as  ground  energy  flux,  while  Ze  is  ob¬ 
served  near  the  freezing  level. 

3.2  THE  EFFECT  OF  MELTING  ON  KINETIC  ENERGY  FLUX 

Equation  3.14  assumes  that  no  melting  occurs  in 
the  hail's  descent  from  the  freezing  level  to  the  ground. 
However,  the  diameters  of  hailstones,  and  consequently  the 
kinetic  energy  flux,  can  change  appreciably  from  melting. 
To  determine  the  effect,  a  simplified  model  of  hail  melt¬ 
ing  was  developed. 

Mason  (1971) ,  based  on  studies  by  Ludlam  (1958) 
and  Macklin  (1963) ,  gives  the  relationships  between  the 
radius  of  a  hailstone  at  the  freezing  level  and  its  radius 
upon  reaching  the  ground  after  melting  as 

R0*  -  Rg^  =  2.  23x10“  3b|’Z°Pr^KATpa^y~^dz  3.15 

where  R0  and  Rg  are  the  radii,  in  centimeters,  of  a  hail¬ 
stone  at  the  freezing  level  and  at  the  ground,  respective¬ 
ly.  b  is  a  shape  factor,  which  for  spheres  is  0.658 
(Macklin,  1963).  Z0  is  the  height  above  ground  of  the 
freezing  level  in  kilometers.  Pr  is  the  Prandtl  number, 


which  has  a  value  of  0.71  for  air.  K  is  the  thermal  con¬ 
ductivity  of  air,  in  Jm" 1  see- 1 ° K" 1 .  In  the  range  of  tem¬ 
peratures  considered  here,  K  is  approximately  2.4xl0~2. 

AT  is  the  air  temperature  in  °C  (a  function  of  Zc  pa  is 
the  ambient  air  density  in  g  m" 3 ,  and  p  is  the  dynamic 
viscosity  of  air,  which  is  approximately  1.75xl0~4g  cm_1sec~ 
for  the  range  of  temperatures  assumed  in  this  model.  To 
a  good  approximation,  the  right-hand  side  of  Equation  3.15 
is  proportional  to  the  surface  temperature,  in  °C,  and  the 
square  of  the  depth  of  the  melting  layer.  Conditions  typi¬ 
cal  of  Alberta  hailstorms  were  taken  to  be  pa  =  1.0kg  m" 3 , 

AT  =  18°C  at  the  surface  with  a  constant  lapse  rate  from 
the  surface  to  the  freezing  level,  and  z 0 =  3km.  Using  these 
values,  Equation  3.15  becomes 

7  2 

R04  -  Rq4  =  0.242  3.16 

In  terms  of  the  initial  and  final  diameters,  D0  and  Dq ,  in 
centimeters,  Equation  3.16  becomes 

D04  -  DG4  =  0.814  3.16a 

Equation  3.16a  assumes  that  the  hail  is  a  sphere  falling 
in  still,  clear  air.  The  actual  melting  may  be  different 
from  what  Equation  3.16a  implies,  when  the  falling  stone 
is  caught  in  an  updraft  or  downdraft,  but  this  effect  is 


difficult  to  assess. 
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Re-writing  Equation  3.16a  gives 


7  4 


Dg  =  (Do5  -  0.814)7 


3.17 


The  kinetic  energy  of  a  hailstone  of  diameter  D  is: 


E  =  -jCp  jk2D‘* 


3.18 


If  k  (the  air  density  dependent  terminal  velocity  coeffi¬ 
cient),  pj  and  D  are  in  cgs  units,  substituting  Equation 
3.17  into  Equation  3.18  yields 


7 

4 


1  6 


E  =  ~  p;k2G  (D  o 


0.814)“ 


3.19 


which  is  the  ground  level  kinetic  energy  in  ergs,  of  the 
hailstone  in  terms  of  its  freezing  level  diameter,  D0.  k 
takes  on  its  ground  level  value,  kG,  in  Equation  3.19. 

The  problem  of  the  change  in  k  from  the  freezing  level  to 
the  ground  as  it  influences  radar  and  hailpad  estimates  of 
E  will  be  treated  in  Chapter  4.  The  velocity  of  the  melted 
hailstone  at  the  ground,  using  Equation  3.2  and  3.17  is 


v  =  kG(D04  -  0 . 814) 7  --c.g.s.  units--  3.20 


7  2 


The  kinetic  energy  of  the  hailstone  at  cloud  level  is  given 
by  Equation  3.18  when  D  =  D0.  An  ensemble  of  hailstones 
is  assumed  to  have  a  Marshal 1 - Palmer  size  distribution. 


I 
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N0exp(-AD0),  before  melting.  N0  here  has  units  of  cm-4. 

The  kinetic  energy  flux  of  the  ensemble  at  the  ground  is 
the  product  of: 

a.  the  number  of  unmelted  hailstones  per  unit  volume 
as  a  function  of  diameter  (the  Marshall- Palmer  formula); 

b.  the  initial  kinetic  energy  per  stone, 

Y2p i k o  Do  ; 

c.  the  factor  by  which  the  energy  per  stone  changes 
due  to  melting  (the  ratio  of  Equation  3.19  to  Equation  3.18 
in  terms  of  Do)? 


2 


1  6 

-  .  8 14  )_7~ 
Do4 


and 


d.  the  velocity  at  the  ground  per  stone, 

kG(D0"4  -  0.814)7, 


integrated  over  all  diameters  of  hailstones  that  reach  the 
ground .  Thus , 


Eground  {domT^  N°exP("AD°) 


TT  i  2-n  4^3 

T2P'k°  D°  ko^ 


(D0* 


U_6. 

-  0.814)  7 

W4 


7  2 

kG(D04  -  O.814)7dD0, 


or 


35 


(Do 

'UoMIN 


7 

4 


1  8 


0.814)  7  dD0  3.21 


Eground  is  in  erSs  cm_2sec 


and  Dqmin  is  the  diameter  of 


the  smallest  hailstone  which  just  reaches  the  ground  after 
melting.  The  energy  flux  of  the  hail  before  melting  is 


E  c  I  oud  =  -j^P  jko3  N0/D°MAXexp(-AD„)Do 4 '  5dD 


3.22 


The  ratio  of  Equation  3.21  to  3.22  will  therefore  give  the 
ground  energy  flux  as  a  fraction  of  the  cloud  energy  flux 
(at  the  freezing  level)  .  D0 MAX  as  set  at  Hem  to  avoid 
problems  with  supercritical  fall  regimes.  Fortunately  for 
this  model,  hailstones  with  diameters  greater  than  11cm 
are  extremely  rare.  From  Equation  3.17,  D0^|n  is  found  to 
be  0.90cm.  Equations  3.21  and  3.22  were  evaluated  by 
Simpson’s  Rule  to  obtain  the  ratio  of  the  ground  energy  to 
the  cloud  energy  for  values  of  A  from  7.5  to  4.0cm-1.  The 
reasonableness  of  this  range  of  A  was  checked  by  the  for¬ 
mula 


which  relates  Wj,  the  water  mass  concentration  of  hail  in 
the  cloud,  to  A  when  a  Marshall- Palmer  hail  size  distribu¬ 
tion  is  assumed.  Here  Wj  and  pj  are  in  the  same  units, 
gm- 3 ,  and  N0  =  3xl04m"4.  When  A  =  4.0cm-1,  Wj  is  3.31gm-3, 


which  is  close  to  the  highest  water  mass  concentration  of 
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hail  believed  to  exist  in  hail-producing  clouds  (Charlton, 
1978).  At  a  A  of  7.5cm-1,  Wj  =  0.40gm-3.  The  ratio  of 
Equations  3.21  and  3.22  is  plotted,  ignoring  variations 
in  k,  as  a  dotted  line  in  Figure  3.1.  The  correct  form  of 
Equation  2.14  would  include  this  ratio  as  a  coefficient. 

A  true  verification  of  this  melting  model  would 
require  direct  measurements  of  the  vertical  hail  energy 
flux  at  the  freezing  level,  and  comparing  them  to  the  ener¬ 
gy  flux  recorded  by  hailpads  on  the  ground  immediately 
below.  This  procedure  was  impractical. 

A  partial  check  with  hailpads  on  the  validity  of 
the  model  was  done  by  translating  the  hailstone  diameters 
at  the  ground  into  the  corresponding  diameters  at  the  freez¬ 
ing  level  by  Equation  3.16a,  calculating  the  kinetic  energy 
at  each  diameter  by  Equation  3.18,  and  summing  over  the  size 
spectrum  observed  from  the  hailpad,  making  sure  that  the 
resulting  quantity  was  in  the  appropriate  units  (SI) .  The 
hailstone  diameters  at  the  ground  were  determined  by  the 
calibration  formulae,  Equations  2.3  and  2.4,  of  Strong  and 
Lozowski  (1977)  .  The  contribution  to  the  total  kinetic 
energy  flux  at  the  freezing  level  by  the  hailstones  that 
melted  completely  before  reaching  the  ground  was 

-pjP  |  k 3  N  o  /°'9°  exp(-AD0)D04'5  dD0 
I  z  o 

This  quantity  was  found  by  rough  calculation  to  be  negli¬ 
gible  for  all  hailpads  used  in  this  procedure.  The  kinetic 
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energy  flux  was  obtained  by  dividing  the  total  kinetic 
energy  calculated  above  by  the  duration  of  the  hailfall, 
in  seconds,  reported  by  the  farmer  operating  the  hailpad. 

A  total  of  53  hailpads  from  August  7  and  18  carried  hail¬ 
fall  duration  data  and  had  a  flux  at  the  ground  high  enough 
to  be  usefully  plotted.  Such  estimates  of  hail  energy  flux 
are  very  crude  because  the  hailfall  durations  are  given 
roughly  to  the  nearest  minute.  The  ratios  of  these  esti¬ 
mates  of  ground  flux  to  f ree z ing- level  flux  are  plotted 
as  points  in  Figure  3.1.  An  estimated  best  fit  to  these 
points  is  indicated  by  the  dashed  line.  Figure  3.1  shows 
that  the  melting  model  correctly  predicts  the  trend  (with 
increasing  kinetic  energy  flux)  and  order  of  magnitude  of 
the  ratio  of  the  flux  at  the  ground  to  the  flux  at  the 
freezing  level,  assuming  that  the  calculated  fluxes  at  the 
freezing  level  are  close  to  the  real  fluxes.  However, 
Figure  3.1  indicates  that  the  model  may  overestimate  the 
effect  of  melting  on  hail  kinetic  energy  flux.  Both  the 
model  and  the  estimates  show  that  the  fraction  of  the 
kinetic  energy  flux  retained  by  the  hail  has  increased 
with  increasing  energy  flux  at  the  ground. 

A  comparison  of  hailpad  kinetic  energy  densities 
with  the  corresponding  estimated  kinetic  energy  densities 
at  the  freezing  level  was  done  by  following  the  procedure 
used  for  the  kinetic  energy  fluxes,  except  without  dividing 
through  by  the  hailfall  durations.  This  method  was  applied 
to  all  the  hailpads  from  August  7  and  18  that  recorded 


RATIO  OF  GROUND  KINETIC  ENERGY  FLUX 
TO  CLOUD  KINETIC  ENERGY  FLUX 
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Figure  3.1:  The  ratio  of  the  ground-level  kinetic  energy 

flux  of  hail  to  cloud  (freezing  level)  kinetic 
energy  flux  for  53  hailpads  from  August  7 
(circles)  and  August  18,  1974  (triangles). 

This  ratio  is  plotted  on  the  vertical  scale 
against  ground-level  kinetic  energy  flux  in 
Jm'2sec' 1 .  Upper  curve  is  the  best  fit, 
drawn  by  eye,  to  the  hailpads.  Lower  curve 
gives  the  theoretical  ratio  for  the  given 
ground-level  energy  flux. 


' 
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measurable  hailfalls.  These  pads  are  plotted  in  Figure  3.2, 
except  for  those  with  energy  densities  too  low  to  appear 
on  the  graph.  Again,  a  trend  toward  higher  retention  of 
kinetic  energy  at  high  overall  energy  densitites  is  evident 
in  Figure  3.2.  An  estimated  best-fit  line  to  the  plotted 
points  is  drawn  in.  The  ratio  of  ground  to  freezing  level 
energy  density  reached  0.2  at  a  ground-level  energy  density 
of  about  80Jm~2.  The  value  of  0.2  was  chosen  as  the  melt¬ 
ing  coefficient  to  be  added  to  Equation  3.14.  It  should 
again  be  emphasized  that  the  effect  of  the  slowing-down  of 
the  hailstones  in  the  denser  air  near  the  surface  was  not 
taken  into  account  (i.e.,  "kM  variations)  in  producing 
Figure  3.2. 

The  choice  of  0.2  as  the  melting  factor  was  some¬ 
what  arbitrary  since  the  effect  is  markedly  non-linear  and 
depends  strongly  on  the  f ree z ing- leve 1  (or  ground- leve 1) 
kinetic  energy  flux  of  hail.  The  value  was  chosen  because 
it  is  reasonably  accurate  at  high  ground- level  energy  den¬ 
sities,  which  is  the  main  situation  of  interest  in  this 
study.  Thus,  Equation  3.14  becomes 

Ze 

E  =  0.20  x  3.40  x  10‘6  x  1012*1  3.24 


or 

Ze 


E  =  6.80  x  10" 7  x  10  12, 1 


3.25 


< 


RATIO  OF  GROUND  KINETIC  ENERGY  DENSITY 
TO  CLOUD  KINETIC  ENERGY  DENSITY 
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Figure  3.2:  The  ratio  of  the  ground-level  kinetic  energy 

density  of  hail  to  cloud  (freezing  level) 
kinetic  energy  density  for  62  hailpads  from 
August  7  (circles)  and  August  18,  1974  (tri¬ 
angles)  .  The  ratio  is  plotted  on  the  verti¬ 
cal  scale  against  the  ground-level  kinetic 
energy  density  in  Jm  2.  The  curve  is  the 
estimated  best  fit  to  the  hailpad  points. 
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with  E  in  Jm-2sec-1  and  Ze  in  dBZ.  It  will  be  shown  in 
Chapter  4  that  the  adjustment  for  the  change  in  air  den¬ 
sity  from  the  freezing  level  to  the  ground  is  a  simple 
multiplicative  constant. 


CHAPTER  4 


METHOD  OF  RADAR- HAI LPAD  INTERPOLATION 
AND  PRELIMINARY  RESULTS 

4.1  DEVELOPMENT  OF  OPERATIONAL  FORMULA 

Equation  3.25  will  be  used  to  develop  an  opera¬ 
tional  procedure  that  interpolates  the  energy  density  of 
hail  at  points  without  hailpads  when  radar  and  hailpads 
are  available.  The  procedure  is  based  on  a  time  integra¬ 
tion  of  Equation  3.25,  but  three  important  considerations 
entered  in  the  choice  of  integration  method. 

The  first  is  that  not  all  altitudes  in  the  storm 
cloud  will  provide  radar  information  that  correlates 
equally  well  with  the  hailfall  at  the  ground.  At  very  high 
altitudes,  the  hail  is  in  the  process  of  formation,  giv¬ 
ing  it  a  size  spectrum  that  may  be  different  from  the  one 
it  will  have  by  the  time  it  reaches  the  ground.  The  hail¬ 
stone’s  surface  will  furthermore  be  dry  and  cold,  rather 
than  wet  and  melting.  Also,  because  hail  is  carried  hori¬ 
zontally  to  some  extent  by  wind,  hail  at  high  altitudes  is 
less  likely  to  fall  on  the  land  location  intended  as  the 
detection  area.  For  these  reasons,  it  will  have  different 
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radar  reflectivity  properties  from  those  it  should  have  at 
lower  levels.  On  the  other  hand,  radar  information  from 
very  close  to  the  ground  will  be  affected  by  phenomena 
such  as  ground  clutter.  Time  integration  of  Equation 
3.25  was  therefore  carried  out  only  for  echoes  at  heights 
between  1.5  and  4  kilometers. 

The  radar  executes  a  complete  scan  cycle  every 
211  seconds,  making  a  complete  revolution  approximately 
every  14  seconds.  The  elevation  starts  at  0  degrees  and 
increases  1.5  degrees  with  every  revolution  up  to  21 
degrees,  and  returning  to  0  degrees  to  begin  the  next  scan 
cycle.  From  one  to  five  echoes  per  scan  cycle  are  avail¬ 
able  at  altitudes  between  1.5  and  4  kilometers,  depending 
on  the  range  of  the  land  location.  The  hail  corresponding  t 
these  echoes  was  assumed  to  be  representative  of  the  hail 
that  fell  during  the  entire  3§  minute  scan  cycle.  All 
echoes  within  the  interval  were  weighted  equally,  in  the 
absence  of  any  other  information  on  their  representative¬ 
ness  of  hailfall  at  the  ground  [see  below  for  details) . 

The  second  consideration  is  that  it  is  desirable 
to  eliminate  the  contribution  of  rain  to  the  echoes.  Rain 
contributes  most  of  the  returned  power  in  echoes  of  less 
than  about  35  dBZ  (Barge,  1972,  1977).  For  this  reason 
a  threshold  equivalent  reflectivity  was  assumed  such  that 
echoes  less  than  the  threshold  did  not  enter  the  calculation 
of  the  hail  energy  density.  Thresholds  of  30,  35,  40  and 
45  dBZ  were  tested  as  a  crude  method  of  eliminating  the 
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rain  contribution. 

Finally,  falling  hailstones  will  lose  kineti 
energy  as  a  result  of  being  slowed  down  by  the  increas 
air  density.  Recalling  Equation  3.18,  the  kinetic  ene 
of  a  hailstone  of  a  given  diameter  is  proportional  to 
Thus  we  have,  for  any  hailstone, 


c 
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rgy 
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FT 


4.1 


Er,  which  is  the  total  kinetic  energy  density  of 
the  hail  at  the  freezing  level,  estimated  by  radar,  is 
given  by 

N  . 

Ed  -  E  E;At  4.2 

i  -  1  1 


when  N  observations  are  taken,  At  is  the  mean  time  interval  be- 
tween  in-cloud  radar  observations  and  Ej  is  the  energy  flux 
of  the  hail.  Since  the  kinetic  energy  density  of  all 
hailstones  is  altered  by  the  factor  (kG2/k02),  the  total 
energy  density  at  the  ground,  estimated  by  radar  is 


Er 


Z  E  i  At 

K  o  i  =  1 


4.3 


apart  from  any  melting  effect.  It  should  be  noted  that 
energy  flux  density  varies  as  k3  between  cloud  and  ground 
(Equations  3.2  and  3.22)  but  energy  density  varies  as  k2 
because  the  time  of  hailfall  is  different  between  cloud 
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and  ground.  Since  k  near  the  freezing  level  is  159m2sec_1, 

1 

and  at  the  ground  is  about  133m2sec_1,  the  formula  becomes 


^  1332  N  • 

Er  "  I59z  j l ^ ' At 


4.4 


or 


N 

Er  =  0.70  E  E;At  4.5 

i  =  1 


The  complete  formula,  based  on  Equation  3.25  and  including 
the  melting  effect,  now  becomes 


-j-  -  -r 

Er  =  0.70  X  6  .  80  X  10-  7  tt — — -S —  E  ei0TiTr 

nSCAN  t  =  0 


or 


R 
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4.6 


where  [0,te]  is  the  time  interval  in  the  radar  record  during 
which  echoes  above  the  threshold  value  occurred.  is 

the  number  of  scans  appearing  between  1.5  and  4.0  kilo¬ 
meters.  Ts  is  the  time,  in  seconds,  it  takes  the  radar  to 
undergo  a  complete  scan  cycle,  which  in  the  case  of  the 
Alberta  Hail  Project  C-band  is  211  seconds.  Ze  is  the 
strength  of  any  echo,  in  dBZ,  within  the  height  interval 
and  equal  to  or  greater  than  the  threshold  value.  Er  has 
units  of  joule  m“  2  .  The  quantity  Er  will  from  now  onward 
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be  referred  to  as  the  "raw  radar  estimate"  of  hailfall  kine¬ 
tic  energy  density. 

A  first  attempt  at  using  Equation  4.6  was  made 
with  the  original  radar  data  (with  calibration  constant 
68.4dB),  and  comparing  the  results  with  the  hailpads  using 
the  older  calibration  scheme.  The  use  of  Equation  4.6  is 
illustrated  in  Figure  4.1(a).  This  is  the  sample  of  PLOTPAD 
made  with  the  original  data,  representing  the  echo  history 
over  hailpad  H37  of  August  7,  1974.  Here,  there  are  2 
scans  between  1.5  and  4  kilometers;  hence  N  SCAN  =  2.  A 
program  called  PADEN,  similar  to  PLOTPAD,  was  used  to  cal¬ 
culate  Er.  Using  thresholds  of  various  values,  Er  worked 
out  to  the  following: 


Id  (dBZ) 

Er  (J  m 

45 

0.0 

40 

0.099 

35 

0.321 

30 

0.437 

The  following  sample  calculation  shows  how  the  values 
listed  above  were  derived: 

In  Figure  4.1(a)  at  the  35dBZ  threshold,  the  follow- 
echoes  were  found  to  be  within  the  required  height 
i nterva I : 

35,  37,  35,  40,  39  dBZ  (indicated  by  boxes) 

Nscan  =  2;  Ts 
Er  =  4.76  x  10 

3  9  - 

+  io12 • 1 


=211  sec, 


-7 


21 


2 


3  5 


37 


3  5 


4  0 


10 


12.1 


+  10 


12.1 


+  10 


1  2 


+  10 


12.1 


Figure  4.1a,b:  This  is  the  same  PLOTPAD  output  shown  in 

Figures  2.3  and  2.4  for  (a)  the  original 
radar  data  and  (b)  the  revised  radar  data, 
except  showing  the  echoes  which  contribute 
to  the  calculation  of  equation  4.6.  In 
both  (a)  and  (b) ,  the  echoes  are  outlined 
by  boxes.  The  vertical  interval  in  which 
echoes  can  contribute  to  Equation  4.6  is 
shown  by  the  bracket  at  the  extreme  right 
of  each  frame.  In  these  examples,  two 
scans  fall  between  heights  of  1.5  and  4.0 
kilometers,  so  that  NsCAN  equals  2.  The 
threshold  reflectivity  is  35dBZ. 
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5.00  x  10” 5 [781  +  1142  +  781  +  2022  +  1671] 

=  0.321  J  rrf2. 

The  hailpad  value  of  energy  density,  according  to  the  ori¬ 
ginal  analysis,  was  27.3  Jm“  2  .  The  raw  radar  estimate  by 
the  above  calculation  is  about  85  times  smaller  than  the 
hailpad  value.  A  discrepancy  of  this  size  occurred  sys¬ 
tematically  throughout  the  radar  and  hailpad  data.  This 
prompted  a  re-evaluation  of  all  the  data  and  led  to  the 
revisions  described  in  Chapter  2. 

Figure  4.1(b)  illustrates  the  PLOTPAD  output 
based  on  the  revised  radar  data  (calibration  constant 
73.2dB)  for  hailpad  H37.  Repeating  the  calculation,  Er 
is  now  found  to  be  the  following: 


Threshold  (dBZ)  ErUitt2) 

45  0.281 
40  0.847 
35  1.143 
30  1.166 


The  sample  calculation  is  again  illustrated  for  the  35  dBZ 
threshold : 


Echoes  at  or  above  35dBZ  (see  Figure  4.1b) 
38,  40,  37,  37,  42,  37,  45,  40,  43,  35  dBZ 

nSCAN  =  2;  Ts  = 


2 1 1  sec 
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3  8  4  0  3  7  3  7 


Er  =  5.00  X 

10 

-  5[  1 0 1 2  '  1 

+ 

1012 • 1  + 

1012  •  1 

+ 

1012* 1 

+ 

42 

3  7 

4  5 

4  0 

4  3 

3  5 

1012*  1 

+ 

1012*1  + 

10 

12- 1  +  10 

12.1  + 

10 

12.1  + 

io12-1 

=  1 .  143  X  Jrrf 2 
-  1.1  J  m- 2 

The  corresponding  hailpad  energy  density,  using  the  revised 
calibration  of  Strong  and  Lozowski  (1977),  is  8.7  Jm" 2 . 

The  raw  radar  estimate  is  too  small  by  a  factor  of  8.  In 
general,  the  radar-based  calculation  still  underestimates 
the  hailpad  energy  density.  Radar  beam  attenuation  would 
account  for  this  effect,  and  will  be  discussed  in  a  later 
chapter . 

The  radar  estimates,  using  a  30dBZ  threshold, 
and  the  corresponding  hailpad  values  of  kinetic  energy  den¬ 
sity  were  plotted  for  August  7  and  18,  using  the  revised 
data  [Figures  4.2(a)  and  4.2(b)].  If  the  radar  estimation 
method  is  valid,  some  sort  of  linear  relationship  between 
the  radar  estimate  and  the  corresponding  hailpad  value  would 
be  expected.  The  desired  relationship  does  appear  in  both 
Figures  4.2(a)  and  4.2(b),  although  in  the  case  of  August  7, 
the  relationship  is  obscured  to  some  extent  by  the  strong 
scatter.  An  important  feature  of  the  August  7  data  is  the 
large  number  of  zero  or  very  small  raw  radar  estimates  as¬ 
sociated  with  appreciable  hailfalls.  The  relationship  be¬ 
tween  the  raw  radar  estimate  and  the  hailpad  value  of 
energy  density  is  much  clearer  in  the  case  of  August  18. 
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Figure  4.2a,  b:  Raw  radar  estimates  of  hail  energy  density 

are  plotted  against  the  energy  densities  as 
determined  by  the  corresponding  hailpad 
analyses  for  (a)  August  7  and  (b)  August  18, 
1974.  The  raw  radar  estimate  (ER,  in  Jm~ 2 )  is 
given  along  the  horizontal  scale,  while  the 
hailpad  value  (ER,  in  Jm“2)  is  given  along  the 
vertical  scale.  The  leas t- squares  fitting  line 
for  E|_j  and  ER  is  shown  in  both  cases  and  the 
values  of  ’a"  and  "b"  in  E  =  aER  +  b,  deter¬ 
mined  by  the. lines,  are  also  indicated.  Note 
that  the  horizontal  scale  has  been  expanded 
for  greater  clarity. 
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Occasionally,  the  raw  radar  estimates  and  hailpad 
values  did  not  conform  at  all  to  the  overall  pattern.  For 
three  cases  in  August  7  and  one  in  August  18,  heavy  hail 
fell  but  the  raw  radar  estimate  was  very  low,  a  result  of 
anomalously  low  echo  strengths.  These  points  are  circled 
in  Figures  4.2(a)  and  4.2(b).  Unusually  severe  attenuation 
of  the  radar  beam  is  the  probable  cause,  although  this  is 
difficult  to  prove.  In  one  other  case,  in  August  7,  the 
radar  estimate  was  very  high  relative  to  the  hailpad  value. 
This  anomaly  is  more  difficult  to  explain.  This  point  cor¬ 
responds  to  hailpad  G27.  The  output  of  PLOTPAD  from  the 
hailpad  site  and  from  the  eight  surrounding  quarter- sec t ions 
was  examined.  Echo-trace  maps  for  August  7  were  also  con¬ 
sulted.  These  maps  are  produced  by  the  Hail  Project,  and 
are  drawn  every  15  minutes  at  the  1.5  degree  elevation.  Gener¬ 
al  inferences  can  be  made  from  these  maps  about  the  direc¬ 
tion  of  a  storm's  movement  and  the  stage  of  its  development 
at  a  particular  time.  The  PLOTPAD  charts  and  the  echo- 
trace  maps  indicate  that  the  storm  that  affected  hailpad 
G27  was  undergoing  rapid  changes  at  the  time  it  was  passing 
over  the  hailpad.  The  storm  had  a  long,  narrow  configura¬ 
tion  and  probably  consisted  of  several  cells  in  a  line. 

Two  distinct  cells  show  in  the  PLOTPAD  charts.  The  larger 
cell,  the  first  one  to  cross  the  hailpad,  was  declining. 

If  it  still  contained  some  small  ice  particles,  it  would 
produce  a  strong  radar  echo,  although  it  would  produce 
only  light  hail  at  the  ground.  The  result  would  be  an 
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anomalously  high  raw  radar  estimate  of  hail  energy  density. 

The  linear  correlation  coefficient  r,  between  the 
hailpad  values  of  kinetic  energy  density  and  the  raw  radar 
estimates,  where 

r  =  £  (Er;  -  Er)  C E h ;  -  Er) 

s  RSH 

was  calculated  for  August  7  and  18  at  the  45,  40,  35  and 
30  dBZ  thresholds.  In  Equation  4.7,  Er;  and  Er;  are  the  raw 
radar  estimates  and  corresponding  hailpad  values  of  kinetic 
energy  density.  Er  and  Er  are  the  mean  values  of  these 
quantities,  and  sr  and  sr  are  the  standard  deviations.  r 
measures  the  consistency  of  the  relationship  between  Er 
and  Er.  A  value  of  r  of  ±1  indicates  that  there  is  an  exact 
linear  relationship,  while  a  value  of  0  indicates  that  there 
is  no  relationship  between  Er  and  Er,  i.e.,  they  are  uncor¬ 
related.  The  anomalous  hailpads  mentioned  above  were  not  in¬ 
cluded  in  the  calculation  of  the  correlation  coefficient.  The 
correlation  coefficient  was  calculated  and  listed  in  Table  4.1. 
These  results  agree  with  the  results  of  Waldvogel  et  al  (1978) ,  who  re¬ 
port  a  correlation  of  0.66  between  radar  and  hailpad  estimates  of  energy 
density  for  286  hailpads. 

If  the  threshold  echo  levels  were  really  useful 
in  improving  the  accuracy  of  the  raw  radar  estimate  by 
eliminating  rain  echoes,  the  correlation  would  be  notice¬ 
ably  lower  at  low  threshold  levels.  No  appreciable  dif¬ 
ference  was  found  in  the  correlation  coefficients  on  a  given  day  in 
any  of  the  cases,  except  for  the  45  dBZ  threshold  on  Aug¬ 
ust  7,  which  was  somewhat  lower  than  the  other  correlations 


■  '■  •  •  ■ 


Table  4.1 


LINEAR  CORRELATION  COEFFICIENT,  r,  OF  RADAR  ESTIMATE 
AND  HAILPAD  VALUE  OF  ENERGY  DENSITY  FOR  AUGUST  7  AND 
18,  1974  AT  RADAR  CALCULATION  THRESHOLDS  OF  45,  40, 

35  AND  30  dBZ 


Threshold  (dBZ) 

r 

August  7 

August  18 

45 

0.547 

0.891 

40 

0.630 

0.875 

35 

0.636 

0.870 

30 

0.637 

0.873 
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on  that  day.  It  was  concluded  that  the 
used  here  were  not  effective  in  improv 
the  radar  estimate.  However,  note  that 
between  the  radar  estimates  and  hailpad 
er  on  August  18  than  on  August  7.  Give 
tions,  the  correlations  between  the  raw 
hailpad  values  were  significant  at  all 
both  days.  Echoes  below  35  dBZ  did  not 
ficantly  to  the  total  raw  radar  estimat 
ponential  increase  in  the  contribution 
echoes  with  increasing  dBZ. 

The  final  step  in  obtaining  a 
estimator  for  hail  energy  density  is  to 
regression  of  the  raw  radar  estimates  a 
values  at  the  chosen  threshold  for  the  s 
The  final  adjusted  estimate,  E,  is  obta 
squares  fit  of  the  raw  radar  estimate  a 
for  the  given  storm-day: 


threshold  echo  levels 
ing  the  accuracy  of 
the  correlation 
values  is  strong- 
n  these  qualifica- 
radar  estimates  and 
threshold  levels  on 
contribute  signi- 
es,  due  to  the  ex- 
from  individual 

storm- ad j  us  ted 
perform  a  linear 
gainst  the  hailpad 
torm  of  interest, 
ined  from  the  least- 
gainst  hailpad  value 


E  =  a  Er  +  b 


4.8 


where  a  and  b  are  fitting  parameters.  Figures  4.2(a)  and 
4.2(b)  show  the  leas t- squares  fitting  lines  for  August  7 
and  18  at  the  30dBZ  threshold.  In  deriving  Equation  4.8 
the  four  anomalous  hailpads  from  August  7  and  the  one  from 
August  18  were  deleted.  E  can  be  calculated  for  any  loca¬ 
tion  by  Equation  4.8  provided  there  is  complete  radar  data 
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data  available  to  calculate  Er.  If  at  any  time  Equation 
4.8  yields  a  negative  value  for  E,  E  can  be  set  to  0.0  Jm“  2 
In  situations  where  there  are  not  enough  hailpads  to  set  up 
a  true  leas t- squares  fit,  average  values  of  a  and  b  based 
on  many  storms  may  be  substituted  in  Equation  4.8.  In 
general,  the  more  hailpads  that  are  available,  the  more 
reliable  will  be  the  adjusted  hail  energy  density  estimate, 
E.  Chapter  6  will  discuss  some  of  the  factors  which  af¬ 
fected  the  observed  relationship  between  the  raw  radar 
estimate  and  the  hailpad  value  of  energy  density.  The 
next  section  will  describe  the  preliminary  results  of  the 
use  of  Equations  4.6  and  4.8. 

4.2  PRELIMINARY  RESULTS:  THE  EFFECT  OF  VARYING  QUANTITY 
OF  INFORMATION  ON  THE  DELINEATION  OF  HAIL  AREAS 

Two  attempts  were  made  at  delineating  hailfall 
for  August  7  and  18  without  the  use  of  radar.  First,  the 
maps  of  Figures  4.3(a)  and  4.3(b)  were  produced  using  the 
hailpad  analyses  only.  The  contours  are  unrealistically 
smooth  and  their  placement  is  extremely  imprecise.  None 
of  the  intricate  structure  known  to  exist  in  most  hail- 
swaths  shows  up  in  these  maps. 

Figures  4.4(a)  and  4.4(b)  are  the  same  as  the 
previous  pair  of  maps,  except  farmers’  maximum  hail  size 
reports  were  also  used.  About  ten  farmers'  reports  per 
township  were  usually  available  from  the  storm  area.  The 
outer  edge  of  the  area  that  received  any  hail  at  all  (shot- 
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Figure  4.3a,b:  Hailfall  energy  densities,  in  Jm“  ,  of 

(a)  August  7  and  (b)  August  18,  as  deter¬ 
mined  by  the  hailpad  values  alone.  The 
small  squares  indicate  the  locations  of 
the  hailpads.  The  radar  site  (QF)  is  at 
X. 
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Figure  4.3b 
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Figure  4.4a,b:  Hailfall  energy  densities,  in  Jm-2,  of 

(a)  August  7  and  (b)  August  18,  determined 
by  the  hailpads  supplemented  with  informa¬ 
tion  from  the  farmers’  reports  of  maximum 
hail  size.  As  in  Figure  4.3,  small  squares 
indicate  hailpad  locations  and  the  X  marks 
the  radar  site.  The  area  inside  the  open 
square  in  Figure  4.4(b)  was  chosen  for 
more  detailed  radar  study. 
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size  being  the  largest  reported  size)  defined  the  lowest 
contour  value,  lOJrn-2,  of  hail  energy  density.  The  far¬ 
mers’  reports  of  maximum  hail  size  were  used  as  a  visual 
aid  in  drawing  successively  higher  contours  of  energy  den 
sity.  At  hailpad  locations,  the  hailpad  value  itself 
determined  the  contour  placement. 

This  use  of  finer-scale  farmers'  report  informa 
tion  already  showed  significant  smaller- scale  features  of 
the  hailswaths.  Major  changes  were  made  from  the  hailpad 
only  maps,  especially  in  the  case  of  August  7,  supporting 
the  belief  that  more  detail  could  be  discovered  at  higher 
information  densities.  A  portion  of  the  August  18  storm, 
covering  144  square  miles,  was  then  chosen  for  more  in¬ 
tensive  radar  and  hailpad  analysis.  This  area,  marked 
by  a  square  in  Figure  4.4(b),  is  shown  at  a  larger  scale 
in  Figure  4.5. 

Using  the  program  PADEN  and  Equation  4.8,  where 
a  =  6.18  and  b  =  0.5  Jm“ 2  ,  radar-based  estimates  of  hail 
energy  density  were  calculated  for  all  quarter- sect  ions 
with  complete  radar  data.  These  estimates  were  plotted 
and  then  contoured  to  give  Figure  4.6.  The  irregular 
eastern  border  of  the  plotted  area  marks  the  limit  of 
complete  radar  data  on  August  18. 

The  most  striking  feature  of  the  map  is  the  in¬ 
tricate  structure  and  strong  variations  in  the  radar- 
estimated  hail  energy  density  field.  While  the  radar- 
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4.5:  Detail  of  August  18,  which  was  indicated  in 

Figure  4.4b.  Numbers  inside  squares  are  the  kinetic 
energy  densities,  in  Jm-2,  given  by  the  hailpads 
at  those  locations.  Maximum  hail  sizes  reported 
by  farmers  are  identified  by  the  following  codes: 

X  =  no  rain  or  hail,  R  =  rain  only,  1  =  shot,  2  = 
pea,  3  =  grape,  4  =  walnut,  5  =  golfball,  6  = 
larger  than  golfball,  and  0  =  maximum  size  unknown. 
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Kinetic  energy  density  of  hailfall  in  JnT 2 , 
based  on  adjusted  radar  estimates  (E  in 
Equation  4.8)  in  the  area  covered  by  Figure 


4.5. 
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based  estimates  are  not  direct  observations,  the  structure 
and  scale  of  the  features  are  generally  in  agreement  with 
the  observations  of  the  spatial  distribution  of  hailfall  at 
the  ground  by  Strong  (1974),  Changnon  (1970)  and  others. 
Close  examination  of  Figure  4.6  shows  that,  apart  from  local 
maxima  and  minima  having  diameters  of  1.5  miles  or  less, 
there  are  some  important  larger-scale  features.  Three 
elongated  areas  of  moderate  (10-40Jm"2)  activity,  oriented 
north-south,  between  2  and  5  kilometers  wide  and  12-16 
kilometers  long  are  visible  in  the  northern  half  of  the 
map.  All  these  features  correspond  in  size  and  shape  to 
the  hailstreaks  described  by  Changnon  (1970). 

Examination  of  echo  trace  maps  of  the  kind  de¬ 
scribed  in  the  last  section  shows  that  two  distinct  storms 
crossed  the  region  of  the  map  of  Figure  4.6.  The  first 
storm  was  maturing  as  it  crossed  the  southern  half  of  the 
region  between  1530  and  1630,  moving  eas t- southeas tward . 

The  second  storm,  at  its  peak  of  development,  moved  through 
the  area  almost  due  east  between  1740  and  1830,  sweeping 
the  northern  three-quarters  of  the  region. 

The  later  storm  contributed  most  of  the  hailfall 
energy  in  the  intermediate  areas  that  experienced  both 
storms.  Superimposed  hailstreaks  from  both  storms  could 
explain  the  small  areas  of  estimated  energy  density  exceed¬ 
ing  80  Jm" 2  which  occurred  near  the  southern  ends  of  the 
north-south  oriented  streaks.  Small  cells  embedded  in  the 
second  storm,  moving  rapidly  northward,  could  explain  the 


Such  cell  movement 


unusual  orientation  of  these  streaks, 
is  commonly  observed  in  squall  lines. 

The  average  energy  density  of  hailfall  over  the 
four  townships  was  18.2  Jm"2  when  calculated  from  the 
hailpad- ad j usted  radar  estimates.  As  given  by  the  mean 
of  the  five  hailpads  in  the  area,  the  average  energy  den¬ 
sity  is  35.8  JnT 2 .  There  is  a  difference  of  a  factor  of 
2  in  these  results. 

The  sparseness  of  farmers'  reports  (about  10 
per  township)  relative  to  the  radar  data  (up  to  144  cal¬ 
culations  per  township)  made  the  consistency  of  the  two 
sources  of  information  difficult  to  ascertain.  In  most 
cases,  there  was  reasonably  good  agreement  between  the 
farmers'  reports  of  maximum  hail  size  and  the  nearby  energy 
density  estimates.  A  major  discrepancy  occurred  in  the 
south  center  of  the  map  region  where  two  farmers  reported 
rain  only,  but  the  radar  indicated  hail  in  excess  of 
80  Jm"2.  Two  land  locations  in  that  locality  yielded 
radar  energy  estimates  of  more  than  120  Jm"2.  The  owner 
of  one  of  the  two  apparently  hard  hit  land  locations,  at 
NW  S18,  T36,  R4,  W5 ,  was  interviewed  by  telephone  on 
April  1,  1977.  He  recalled  receiving  very  heavy  damage  to 
his  crops,  and  seeing  hail  on  the  ground  on  the  evening  of 
August  18,  more  than  two  hours  after  the  storm  had  passed 
over.  His  comments  tend  to  confirm  the  radar  evidence  of 
heavy  hail  at  that  land  location.  This  observation  suggest 
that  radar  data  may  be  very  reliable  in  the  detection  of  se 
vere  hail,  although  this  should  be  confirmed  by  more  study. 


CHAPTER  5 


EVALUATION  OF  THE  RADAR-HAILPAD  ESTIMATION  SCHEME 

AND  ALTERNATIVE  METHODS 


5.1  TRIPLETS 


A  quantitative  assessment  of  the  radar- hai lpad 
estimation  scheme  described  in  the  last  chapter  is  needed. 
That  scheme  will  be  compared  using  two  different  procedures. 
Each  uses  only  the  information  routinely  supplied  with  the 
hailpads.  Some  simple  statistical  tests  will  be  applied  to 
determine  whether  any  of  the  methods  was  able  to  make  valid 
estimates  of  local  energy  density. 

The  hailpads  were  grouped  into  series  of  three, 
called  "triplets'’.  The  three  pads  in  each  triplet  experi¬ 
enced  the  same  storm,  and  lay  within  thirty  degrees  of  a 
straight  line.  The  distance  to  either  end  pad  from  the 
centerpad  did  not  exceed  24  miles  (38  km) .  The  center  pad 
is  referred  to  as  Hc,  and  its  measured  hail  energy  density 
as  Ec .  The  end  pads  are  Hx  and  H2 ,  with  energy  densities  of 
Ei  and  E2,  respectively.  Similarly,  the  distances  of  Hj  and 
H2  to  the  center  pad  are  di  and  d2 .  Estimates,  called  Ec ' , 
were  made  of  the  value  of  Ec  using  Ei  and  E2  and  various 
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combinations  of  radar  information.  Figure  5.1  illustrates 
a  typical  triplet.  A  complete  listing  of  triplets  appears 
in  Appendix  C.  A  description  of  the  estimation  methods  to 
be  tested  follows  in  the  next  section. 

5.2  DESCRIPTION  OF  THE  ESTIMATION  METHODS 

Arithmetic  mean:  this  method  is  based  on  the 
assumption  that  the  energy  density  at  the  center  pad  of 
the  triplet  lies  somewhere  between  the  values  at  either 
end  pad.  No  other  information  is  used,  and  the  center- 
pad  estimate  is  the  arithmetic  mean  of  the  end-pad  values. 
Thus  , 

ECT  =  I (E !  +  E2)  5.1 

Dis tance- weighted  mean:  this  method  assumes  that 

the  energy  density  of  hailfall  between  two  end  pads  varies 
linearly.  The  implicit  assumption  of  this  scheme  is  that 
the  length  scale  of  the  most  significant  features  of  the 
hailfall  pattern  is  larger  than  the  average  distance  be¬ 
tween  hailpads.  The  formula  for  Ec '  is 

r  »  -  d2 E i  +  d ! E 2  r  o 

hc  dj  +  d2 

Radar-hailpad  estimate:  the  formula  has  already 
been  described  in  Chapter  4.  In  this  case  Ec'  is  simply 
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Figure  5.1.:  An  actual  hailpad  triplet,  drawn  to  scale. 

The  squares  represent  the  quarter- sect  ions 
containing  the  hailpads,  rather  than  the  hail- 
pads  themselves.  The  triplet  shown  is  No.  719. 
The  data  for  this  triplet  appears  in  Appendix 
C. 
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the  value  of  E  for  that  location  given 
4.8.  Note  that  this  radar  estimate  is 
hailpad  data  from  the  whole  storm,  not 
the  ends  of  the  triplet. 

Estimation  of  the  center-pad 
metic  and  dis tance -weighted  means  will 
trols  in  the  evaluation  of  the  relative 
radar - hai lpad  method. 


by  Equations  4 
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value  by  the 
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5.3  EVALUATION  AND  COMPARISON 

Estimation  by  the  arithmetic  mean,  the  distance- 
weighted  mean  and  the  radar-hailpad  method  will  be  referred 
to  as  Method  1,  Method  2  and  Method  3,  respectively.  123 
triplets  drawn  from  both  the  regular  and  the  supplementary 
(July  30  and  August  15)  hailpad  sample  were  originally 
chosen  for  the  analysis,  but  because  of  the  limited  radar 
data,  only  60  were  suitable  for  evaluating  all  three  methods. 
The  test  variable  chosen  was  the  error, 


ECT  -  Ec 


or  the  absolute  difference  between  the  estimate  and  actual 
center-pad  energy  density.  The  more  common  measure  of  er¬ 
ror,  the  fractional  error; 


t 


Eel 


could  not  be  used  because  some  of  the  center  pads  recorded 
zero  hail.  Also,  in  terms  of  crop  damage,  the  absolute  er¬ 
ror  is  more  relevant. 

The  center-pad  estimates  and  errors  are  listed 
in  Appendix  D.  The  mean  error  of  all  triplets  was  calcu¬ 
lated  for  Methods  1  and  2,  and  the  mean  error  for  60  triplet 
for  Method  3.  The  mean  over  these  60  triplets  was  also  cal¬ 
culated  for  Methods  1  and  2.  The  standard  deviations  were 
calculated  as  well.  Table  5.1  summarizes  these  results. 

A  number  of  observations  emerged  from  Table  5.1. 
The  first  is  the  high  standard  deviations  of  the  errors 
in  all  three  methods.  In  most  cases  it  is  higher  than  the 
mean  values.  The  second  is  that  the  addition  of  distance 
information  had  no  effect  on  center-pad  estimates  based  only 
on  hailpad  information,  implying  that  the  assumptions  be¬ 
hind  Method  2  are  false.  The  size  scale  of  the  significant 
variations  in  hail  activity  is  much  smaller  than  the  aver¬ 
age  triplet  length  of  about  16  miles,  or  27  km  (Changnon, 
1970)  .  Finally,  the  mean  error  of  the  estimates  by  Method  3 
is  substantially  less  than  those  of  the  other  methods,  indicating  that 
the  radar  provides  real  information  on  local  intensity  of  hail- 
f all .  This  result  requires  further  testing. 

5.4  ANOTHER  STATISTICAL  TEST 

The  procedure  adopted  for  this  final  statistical 
test  was  to  compare  the  size  of  the  mean  error  in  each  of 
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Table  5.1 

MEANS  AND  STANDARD  DEVIATIONS  FOR  THE  ERRORS  LISTED 

IN  APPENDIX  D,  IN  Jm~2 


METHOD  1 

METHOD  2 

METHOD  3 

1  .  A  1 

Ava i 1 ab 1 e  Tr i p 1 ets 

( N= 1 23 ) 

Mean 

22.5 

22.7 

Standard  Deviation 

23.4 

26.9 

2.  Tr 

p 1 ets  Ava i 1 ab 1 e  for 

Method  3 

a . 

August  7  (N=43) 

Mean 

21.5 

23.6 

13.7 

Standard  Deviation 

23.0 

30.5 

21.8 

b . 

August  18  ( N= 17) 
Mean 

16.4 

18.6 

6.4 

Standard  Deviation 

16.7 

19.6 

5.3 

c. 

August  7  and  August 
Mean 

18  (N=60 ) 

20.1 

22.2 

1  1.7 

Standard  Deviation 

21.5 

27.8 

18.7 
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the  methods  with  that  of  a  similarly  calculated  "error" 
from  a  randomized  set  of  numerical  values.  If  an  estima¬ 
tion  method  has  any  validity,  the  mean  error  associated  with 
it  should  be  significantly  smaller  than  the  corresponding 
error  from  the  randomized  data  set.  A  Monte  Carlo  method 
was  used  to  test  this  hypothesis  (McLeish,  1978)  . 

The  data  used  in  this  analysis  were  the  sixty 
triplets  used  to  calculate  the  mean  error  in  Method  3. 

The  energy  densities  observed  at  the  center  pads  of  these 
triplets  were  randomized,  first  by  numbering  them  in  the 
order  in  which  they  appeared  in  Appendix  D.  These  observa¬ 
tions  were  then  re-listed  in  the  order  in  which  their  num¬ 
bers  appeared  in  a  computer- generated  random  permutation  of 
the  integers  from  1  to  60.  The  estimates  of  the  center- 
pad  energy  density  by  each  method,  in  their  original  order, 
were  subtracted  serially  from  the  randomized  observations. 
The  absolute  values  of  the  resulting  figures  will  be  called 
the  "randomized  errors".  The  mean  of  the  actual  errors  for 
each  method  was  compared  to  the  mean  of  the  randomized  er¬ 
rors  . 

The  objective  of  the  randomization  was  to  simu¬ 
late  a  situation  in  which  there  was  total  ignorance  of  the 
true  hailfall  energy  density  pattern,  except  for  the  over¬ 
all  severity  of  the  storms.  It  thus  served  as  a  control 
data  set,  so  that  the  only  differences  between  the  random¬ 
ized  and  the  actual  errors  were  those 
application  of  the  estimation  methods. 


that  resulted  from 
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The  procedure  described  above  was  repeated  10,000 
times  using  a  different  permutation  of  the  center-pad 
values  each  time.  A  count  was  made  of  the  number  of  times 
the  mean  error  of  each  method  equalled  or  exceeded  the  mean 
randomized  error.  This  number,  as  a  fraction  of  the  10,000 
trials,  is  a  measure  of  the  level  of  significance  at  which 
the  mean  actual  error  is  less  than  the  mean  randomized 
error . 

The  "level  of  significance"  is  defined  as  the 
probability  that  a  random  results  will  appear  to  support 
the  hypothesis  under  test  at  least  as  strongly  as  the  ob¬ 
served  result.  For  example,  the  statement  "the  observed 
results  are  different  at  the  5  percent  level  of  signi¬ 
ficance"  means  that  the  probability  is  0.05,  that  results 
arising  purely  from  chance  will  have  a  difference  at  least 
as  large  as  that  between  the  observed  results.  It  follows 
that  the  lower  the  numerical  value  of  the  significance  level,  the 
stronger  is  the  probability  that  the  result  is  physically  meaningful. 

In  statistics,  a  5  percent  significance  level  is  commonly  accepted  as 
a  criterion  for  judging  a  result  as  having  physical  meaning,  instead 
of  arising  from  chance. 

By  implication,  the  fraction  mentioned  above  is  also  a 
measure  of  the  skill  of  the  estimation  procedure.  The  lower  the  frac¬ 
tion,  the  more  significant  is  the  skill.  A  computer  program,  called 
PERMS,  was  composed  by  the  author  of  this  study  to  perform  the  repeated- 
trial  statistical  test.  Table  5.2  lists  (1)  the  level  of  significance  de¬ 
scribed  above  for  the  accuracy  of  each  of  the  center-pad 
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Table  5.2 

RESULTS  OF  TEST  ON  SIZE  OF  MEAN  ERROR  IN  ESTIMATING 
CENTER-PAD  ENERGY,  DESCRIBED  IN  SECTION  5.4 

1.  Level  of  significance  in  reduction  of  mean  error  (percentage) 

2.  Randomized  error 

3.  Observed  mean  error  (as  in  Table  5.1) 

4.  Percentage  reduction  of  mean  error  from  randomized  error 


METHOD  1 

METHOD  2 

METHOD  3 

1 

3.57 

5.63 

0.00 

2 

23.  1 

25.0 

20.2 

3 

20.  1 

22.2 

11.7 

4 

13.  1 

11.1 

42.4 
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estimation  methods;  (2)  the  population  mean  randomized  er¬ 
ror.  This  is  the  error  that  can  be  expected  from  random 
guesses  of  the  center-pad  value.  It  can  be  estimated  by 
taking  the  mean  (over  10,000  trials)  of  the  mean  randomized 
error  within  each  trial.  Finally  (3)  the  observed  mean 
error  and  (4)  the  percentage  difference  between  the  ob¬ 
served  mean  error  and  the  population  mean  randomized  error 
are  listed. 

The  results  in  the  table  show  that  only  Method  3 
achieves  a  substantial  reduction  in  the  mean  error  from 
the  case  of  random  guesses.  Method  1  may  show  significant 
skill  (based  on  a  maximum  5  percent  level  of  significance), 
but  the  reduction  in  error  is  not  large  enough  to  justify 
its  routine  use.  The  0.0  percent  level  of  significance  re¬ 
corded  under  Method  3  indicates  that  the  value  is  smaller 
than  this  statistical  test  can  detect,  i.e.,  the  actual 
level  of  significance  is  probably  smaller  than  0.001  percent 
This  means  that  Method  3  showed  highly  significant  skill. 

On  the  basis  of  the  analyses  in  this  section  and 
section  5.3  it  was  concluded  that  the  combined  radar  and 
hailpad  method, outlined  in  Chapter  4,  provided  a  reliable 
estimate  of  the  energy  density  of  hail  at  places  without 
hailpads . 


CHAPTER  6 


FACTORS  AFFECTING  THE  RELATIONSHIP  BETWEEN 
THE  HAILPAD  ENERGY  DENSITY  AND  THE  RAW  RADAR  ESTIMATE 

As  noted  in  Chapter  4,  the  relationship  between 
E,  the  actual  hail  energy  density  at  the  ground,  and  Er, 
the  radar  estimate,  is  not  E  =  Er,  but  more  like  E  = 
aER  +  b,  where  a  ~  5-10  and  b  is  of  the  order  of  ±l-10Jm~2, 
in  the  case  of  5cm  radar  observations  of  hail.  The  rela¬ 
tionship  is  furthermore  not  perfect,  but  in  practice  shows 
strong  scatter.  This  chapter  will  describe  some  of  the  im¬ 
portant  factors  which  influence  the  relationship  between 
E  and  Er.  These  include  attenuation  of  the  radar  beam  by 
intervening  precipitation  [affecting  mainly  the  value  of 
"a"  mentioned  above) ,  the  statistical  bias  error  due  to 
non-uniform  beam-filling  (probably  affecting  both  "a”  and 
the  scatter),  inaccurate  registration  of  the  receiver,  and 
the  presence  of  spatial  scales  of  hail  below  the  resolution 
of  the  radar. 

6.1  ATTENUATION 

Ordinarily,  calculations  of  the  reflectivity  of 
a  radar  target  are  made  with  the  assumption  that  there  is 
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nothing  between  the  target  particles  and  the  receiver.  In 
practice,  however,  a  radar  beam  passes  through  air,  cloud 
and  precipitation.  The  result  is  a  loss  of  power  in  both 
the  incident  and  reflected  beam  through  absorption  and  scat¬ 
tering.  At  the  microwave  wavelength  involved  in  this  study 
(5.5cm),  precipitation  is  by  far  the  most  important  attenu¬ 
ating  material  (Battan,  1973).  Estimates  of  attenuation 
are  usually  given  by  the  attenuation  coefficient,  in  units 
of  decibels  per  kilometer  (dBkm-1).  The  magnitude  of  the 
attenuation  coefficient  due  to  precipitation  will  depend  on 
the  radar  wavelength  and  the  type,  composition  and  intensity 
of  the  precipitation.  The  total  attenuation  is  the  in¬ 
tegral  of  the  local  attenuation  coefficient,  between  the 
transmitter  and  the  target  and  back. 

The  value  of  the  one-way  attenuation  coefficient 
at  a  given  location  can  be  estimated  by 


kT  =  0.4343  E  Q+  6.1 

vo  I  ume 

(Battan,  1973),  where  Qt  is  the  attenuation  cross-section 
of  each  particle  in  a  unit  volume  and  ky  is  the  attentu- 
ation  coefficient.  As  Equation  6.1  is  written,  Qt  has  units 
of  cm2,  and  the  volume  is  lm3 .  This  gives  kT  in  units  of 
dBkm'1.  The  attenuation  cross-section,  Qy  ,  is  defined  as 
the  area  which,  when  multiplied  by  the  incident  beam  in¬ 
tensity,  gives  the  total  power  taken  from  the  beam  (Battan, 
1973).  Equation  6.1  can  be  written  as 


: 
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Kt  =  0.4343 


vo I ume 


Z  opapNpp 


6.2 


where  op  is  the  attenuation  cross  section  at  particle  dia¬ 
meter  D,  normalized  to  unit  area  (when  D  is  in  cm).  ap  is 
the  cross-sectional  area  of  the  particle  in  cm  ,  and  Npp 
is  the  number  of  particles  of  diameter  D  in  the  unit  volume 
at  a  given  time.  This  summation  is  over  all  diameter  clas¬ 
ses  D  in  the  volume.  Equation  6.2  enables  the  local  at¬ 
tenuation  coefficient  at  the  freezing  level  to  be  calcu¬ 
lated  from  hailpad  data.  Npp  can  be  expressed  in  terms  of 
observable  hailpad  quantities  by  the  following  relation¬ 


ship  (Barge ,  1977)  : 


6 . 3 


where  Np  is  the  total  number  of  hailstones  with  diamater 
D  passing  through  the  unit  volume  in  time  t.  t  is  the  dura¬ 
tion  of  the  hailfall  in  seconds.  Np  is  equivalent  to  the 
number  of  hailstones  of  diameter  D  intercepted  by  the  hail¬ 
pad.  A  is  the  unit  area.  If  A  is  lm2 ,  Np  must  be  expressed 
as  the  number  of  hailstones  intercepted  per  square  meter, 
or  10.76  times  the  number  per  square  foot.  vt  is  the  fall 
velocity  of  a  hailstone  as  a  function  of  its  diameter,  in 


1 


meters  per  second.  At  the  freezing  level,  vy  -  15. 9D2 
(see  Chapter  4)  . 

Equation  6.3  assumes  that  the  hailfall  is  steady - 
state,  although  for  a  rough  calculation  this  assumption  is 
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not  critical.  Substitution  of  Equation  6.3  into  Equation 
6.2  yields 


kT  =  0.4343 


D 

E 

vo I ume 


aDaO 


Nd 

Av-j-t 


6.4 


or 


NDf 
f  00 

sis 

the 


D 

k  T  =  0.4343  E 

vo I ume 

is  the  number  of  hail 
t,  a  value  directly  ob 
described  in  Chapter 
following : 


stones  o 
tainable 
2 .  Equa 


ND  f  •  10.76  ,  „ 

*  - t  6.4a 

lm2  •  15 . 9D2 t 

f  diameter  D  per  square 
from  the  hailpad  analy- 
tion  6.4a  simplifies  to 


0.2353 

t 


D 
,  E 

vo I ume 


NDf  CdD1*5 
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Equation  6.5  gives  kp ,  the  one-way  attenuation.  Two-way 
attenuation,  or  the  attenuation  in  both  the  incident  and 
returned  beam,  would  be  simply  twice  that  value.  Each  dia¬ 
meter  class  represented  in  the  hailpad  must  be  converted 
to  the  equivalent  diameter  class  at  the  freezing  level,  D, 
by  Equation  3.16a.  Values  of  ap  have  been  calculated  by 
Battan  et  al .  (1971)  for  diameters  from  0.1  to  4  centime¬ 

ters,  at  radar  wavelengths  including  5.5cm,  for  various 
proportions  of  ice  and  water  in  the  hail.  The  quantity  of 
water  in  the  hailstone  or  on  its  surface  was  expressed  as 
a  water  shell  of  a  particular  thickness  surrounding  a  core 
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of  solid  ice.  Since  the  hailstones  were  at  the  freezing 
level,  they  would  have  been  wet,  so  for  the  purpose  of 
calculation  they  were  assumed  to  be  coated  with  a  film  of 
water  0.05cm  thick. 

Using  Equation  6.5,  the  attenuation  coefficient 
at  the  freezing  level  above  some  hailpads  was  calculated 
and  listed  in  Table  6.1.  These  values  of  the  attenuation 
coefficient  are  probably  conservative,  because  they  do  not 
take  into  account  the  contribution  to  the  attenuation  from 
heavy  rain  and  from  the  hailstones  which  completely  melted 
before  reaching  the  hailpad.  Attenuation  coefficients  in 
the  S-band  for  the  same  hailpads  were  included  in  Table 
6.1  for  comparison.  Note  that  attenuation  in  the  S-band 
is  usually  an  order  of  magnitude  lower  than  in  the  C-band. 

The  values  calculated  above  for  typical  hail 
situations  under  observation  by  C-band  radar  show  that 
attenuation  was  substantial  on  the  storm-days  in  question. 

It  may  be  possible  in  principle  to  calculate  the  total 
attenuation  along  the  radar  beam.  However,  such  a  calcu¬ 
lation  would  be  too  complex  to  be  practical  and  further¬ 
more,  future  applications  of  the  radar- hai lpad  method  will 
be  made  with  the  10.6cm  S-band  system,  in  which  attenuation 
should  not  be  a  major  problem.  The  C-band  attenuation  co¬ 
efficients  listed  in  Table  6.1  indicate  that  a  4  to  8km 
depth  of  hail  between  the  hailpads  and  the  radar  could  easi¬ 
ly  account  for  the  relative  deficiency  in  echo  strength 
involved  in  the  estimation  of  energy  density  by  C-band 
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Table  6.1 

TWO-WAY  ATTENUATION  COEFFICIENTS,  IN  dB  km"1, 

FOR  SOME  HAILPADS  AT  C-BAND  AND  S-BAND  RADAR  WAVELENGTHS 


radar.  This  deficiency  is  about  6-12dB  in  most  cases  if 
the  radar  estimate  is  a  factor  of  5-10  smaller  than  the 
hailpad  value. 

One  effect  of  attenuation  that  could  contribute 
to  the  scatter  in  the  E-Er  relation  may  be  the  orientation 
of  the  storm  track  relative  to  the  radar.  If  the  storm 
tracks  tangentially  relative  to  the  radar  site,  echoes 
over  some  hailpads  will  be  much  more  strongly  affected  than 
others,  because  some  hailpads  will  always  have  a  large 
depth  of  storm  between  them  and  the  radar.  A  typical  exam¬ 
ple  would  be  hailpad  D37  of  August  7.  When  hail  was  fall¬ 
ing  at  D37,  there  was  always  well  over  20km  of  precipita¬ 
tion  between  the  pad  and  the  radar,  according  to  the  echo- 
trace  map.  D37  received  about  166  Jnr 2  of  hail  energy,  but 
yielded  a  raw  radar  estimate  of  less  than  0 . 5 JrrT  2  .  On  the 
other  hand,  if  the  storm  tracks  directly  toward  or  away 
from  the  radar,  all  the  hailpads  will,  on  the  average,  be 
looked  at  by  the  radar  through  roughly  the  same  depth  of 
storm  as  it  passes  over  each  pad.  The  echoes  over  the 
hailpads  will  therefore  tend  to  be  equally  affected  by 
attenuation . 

This  possibility  may  have  some  bearing  on  the 
difference  in  the  scatter  apparent  in  Figures  4.2a  and  4.2b. 
On  August  7,  the  storm  cells  that  affected  most  of  the  hail 
pads  never  approached  closer  than  40km  to  the  radar,  track¬ 
ing  tangentially  in  an  east-northeastward  direction.  The 
scatter  was  very  severe  in  the  E-ER  relation.  On  August  18 
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the  storms  were  heading  almost  directly  toward  the  radar, 
and  the  E-Er  relation  was  relatively  coherent. 

6.2  NON-UNIFORM  BEAM-FILLING 

Rogers  (1971)  showed  that  in  conditions  of  highly 
non-uniform  filling  of  the  beam  by  precipitation,  radars 
with  a  linear  or  logarithmic  response  will  tend  to  under¬ 
estimate  the  actual  mean  reflectivity  in  the  target  volume. 
If  the  precipitation  reflectivity  varies  by  20dB  within  the 
resolution  volume  of  the  beam  (~lkm3  in  the  case  of  the 
Hail  Project  C-band) ,  the  underestimate  by  a  logarithmic 
receiver  can  be  as  high  as  3dBZ,  and  for  a  30dB  variation, 
as  much  as  6dBZ.  A  pencil- and-paper  calculation,  using 
Equation  4.6  and  data  from  Strong’s  (1974)  hailpad  network 
show  that  such  variations  are  possible  in  a  realistic  hail- 
s  torm . 

Equation  4.6  also  shows  that  a  3dBZ  underestimate 
of  reflectivity  corresponds  to  an  underestimate  of  energy 
density  of  roughly  a  factor  of  1.8,  and  an  underestimate  of 
6dBZ  to  a  factor  of  3.2  in  the  energy  density.  In  marginal 
storms  or  hailfall  areas,  the  actual  average  echo  strength 
may  already  be  near  the  threshold  value  assumed  in  the 
energy  density  calculation.  The  statistical  bias  may  be 
enough  in  such  situations  to  reduce  the  observed  echo 
strength  below  the  threshold.  This  consideration  will 

true  for  any  other  effect,  such  as  attenuation,  which 
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introduces  a  downward  bias  in  the  observed  reflectivity 
and  in  turn  increases  the  error  in  the  energy  density  es¬ 
timate  . 

If  the  statistical  bias  effect  occurred  uni¬ 
formly  throughout  the  storm  it  may,  to  some  extent,  be 
accounted  for  automatically  by  the  regression  procedure 
used  to  form  Equation  4.8.  However,  at  hailswath  boun¬ 
daries  or  during  periods  of  rapid  change  in  the  fine-scale 
structure  of  a  storm,  radar  returns  from  some  locations  may 
be  much  more  prone  to  the  statistical  bias  effect  than  others 

6.3  ERROR  IN  REGISTRATION  OF  RADAR 


Accuracy  of  registration,  or  the  accuracy  with 
which  the  actual  position  of  the  echo  corresponds  to  the 
position  given  by  the  radar  system  output,  is  a  critical 
limitation  on  the  quality  of  any  information  derived  from 
radar  echo  data.  Echoes  in  convective  conditions  can  vary 
by  20dB  or  more  over  distances  on  the  scale  or  the  radar's 
spatial  resolution  of  about  1km.  If  the  position  given  by 
the  radar  output  is  out  by  this  distance,  errors  of  that 
magnitude  in  the  echo  strength  may  result. 

To  determine  the  accuracy  of  registration  of  the 
C-band  radar  on  the  storm  dates  under  study,  the  azimuth 
and  range  of  every  hailpad  to  within  0.1  degree  and  0.1km 
were  recorded  from  an  8-mile- to- 1- inch  scale  map  of  the 
Hail  Project  area,  for  August  7  and  18,  1974.  The  co- 


ordinates  were  measured  from  the  center  of  the  quarter- 
section  containing  the  hailpad  to  the  radar  site.  These 
co-ordinates  were  then  compared  to  the  sets  supplied  with 
the  PLOTPAD  charts,  and  which  also  appear  in  the  output  of 
PLOT.  In  addition,  14  hailpads  from  August  7,  and  14  from 
August  18  were  chosen  at  random  for  an  independent  check. 
This  was  done  by  selecting  an  echo  from  a  PLOTPAD  chart 
and  loating  it  on  a  PPI  segment  generated  by  computer  from 
the  tape  data  containing  the  echo.  The  azimuth  and  range 
of  the  echo  are  then  read  off  from  the  PPI  segment  (see 
Figures  6.1  and  6.2)  and  compared  to  the  map  co-ordinates  as  • 
before.  The  computer- generated  PPI  was  verified  by  check¬ 
ing  the  position  of  Limestone  Mountain,  a  "landmark”  echo 
whose  exact  position  was  known.  The  positions  given  for 
this  and  another  "landmark"  echo  agreed  with  the  true  posi¬ 
tion  on  both  storm  dates.  The  co-ordinates  of  the  hail¬ 
pads  from  both  the  PPI  and  PLOTPAD  were  cross-checked,  and 
in  all  instances  except  one,  the  two  computer  results 
agreed,  within  round-off  limits. 

The  results  of  the  preceding  tests  are  tabulated 
in  Appendix  E.  The  total  registration  error,  in  kilometers 
was  the  Pythagorean  sum  of  the  range  and  tangential  errors. 
The  range  error  was  the  difference,  in  kilometers,  between 
the  range  of  the  hailpad  calculated  by  computer  and  the 
actual  map  range.  The  tangential  error  is  the  difference 
between  the  computer- calculated  azimuth  and  the  map  azimuth 
It  is  expressed  in  kilometers  by  converting  to  radians  and 
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gure  6.1: 


A  PPI  segment  generated. by  the  computer  program 
PPI.  These  data  are  from  August  7,  1974  at 
16:24MDT  at  an  elevation  of  1.5  degrees.  The 
first  two  columns  on  the  left  give  the  range, 
miles  and  kilometers,  of  the  echoes.  The 
azimuths  of  the  echoes,  in  degrees,  and  the 
"ray  numbers"  are  given  along  the  top. ^  A 
"ray"  is  the  set  of  all  bins  along  a  given 
azimuth  (see  page  16  ).  Echoes  are  in  dBZ . 


Figure  6.2:  The  data  of  Figure  6.1  as  they  would  appear  in 

a  true  PPI  display.  The  data  are  outlined  by 
the  small  shaded  sector  near  250  degrees.  The 
radar  site  is  at  "0".  Ranges,  in  miles,  are 
marked  at  the  top  of  the  diagram  while  azi¬ 
muths  from  true  north  are  marked  at  the  left. 
The  curve  indicates  the  approximate  boundary 
of  echoes  lOdBZ  or  higher.  The  area  of  strong 
echoes  in  the  lower  left  of  Figure  6.1  would 
be  in  the  lower  right  corner  of  the  sector. 
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multiplying  by  the  map  range.  The  errors  are  expressed 
by  the  following  equations: 

XR  =  I r PP I  "  rM I  or  XR  =  I r OUT  "  r M I  6.6 

r  r 

XT  ~  57.3  I  a  P  P  I  "  aM  I  or  XT  =  57,3  I  aOUT  '  aM  I  6,7 

XTOTA  L  =  (XR2  +  XT 2 ) 2  6.8 

XT0TAL  tQ6al  registration  error,  in  kilometers. 

Xp  and  Xj  are  the  range  and  tangential  errors,  respective¬ 
ly.  r^  and  a^  are  the  map  range  of  the  hailpad  in  kilo¬ 
meters,  and  the  map  azimuth  in  degrees.  rpp (  and  app (  are 
the  range  in  kilometers  and  the  azimuth  in  degrees  given 
by  the  computer- generated  PPI,  while  rQy -p  and  a0UT  are  the 
range  and  azimuth  given  by  PLOTPAD .  The  errors  are  sum¬ 
marized  in  Table  6.2. 

There  is  an  uncertainty  of  {  mile  1km)  in  a 
hailpad’s  position,  since  its  location  is  given  only  by  the 
quarter- sec t ion  on  which  it  is  set  up.  PLOT  and  PLOTPAD 
selected  the  correct  radar  data  within  the  uncertainty  in 
the  hailpad  location  in  a  majority  of  cases.  The  problems 
caused  by  inaccurate  registration  become  significant  with 
appreciably  larger  errors.  Fortunately,  none  of  the  hail- 
pads  with  complete  radar  data  had  to  be  rejected  because  of 
excessively  poor  registration.  More  sophisticated  echo- 
location  routines,  such  as  those  used  in  the  handling  of 
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Table  6 . 2 


SUMMARY  OF  RADAR  REGISTRATION  ERRORS,  XT0TAL 


Breakdown  of  Errors  by  Size 

C  1  ass 

Number  in  Size  Class 

Percentage 

n  Size  Cl  ass 

Error  Size  Class  (KM) 

<1.0 

1-2 

2-3 

>3.0 

<1.0 

1-2 

2-3 

>3.0 

August  7 

PLOT-MAP  ERRORS  ONLY 

17 

1 1 

2 

0 

56.7 

36.7 

6.7 

0 

DETAILED  ANALYSES: 

PLOT-MAP 

9 

5 

0 

0 

64.3 

35.7 

0 

0 

PP 1 -MAP 

9 

3 

2 

0 

64.3 

21.4 

14.3 

0 

ALL  PLOT-MAP  ERRORS 

26 

16 

2 

0 

59.  1 

36.4 

4.5 

0 

August  18 

PLOT-MAP  ERRORS  ONLY 

18 

5 

0 

2 

72.0 

20.0 

0 

18.0 

DETAILED  ANALYSES: 

PLOT-MAP 

9 

5 

0 

0 

64.3 

35.7 

0 

0 

PP 1 -MAP 

7 

5 

2 

0 

50.0 

35.7 

14.3 

0 

ALL  PLOT-MAP  ERRORS 

27 

10 

0 

2 

69.2 

25.6 

0 

5.2 

TOTAL 

PLOT-MAP  ERRORS  ONLY 

35 

16 

2 

2 

63 . 0 

29.6 

3.7 

3.7 

DETAILED  ANALYSES: 

PLOT-MAP 

18 

10 

0 

0 

64.3 

35.7 

0 

0 

PP 1 -MAP 

16 

8 

4 

0 

57.  1 

28.6 

14.3 

0 

ALL  PLOT-MAP  ERRORS 

53 

26 

2 

2 

63.9 

3  1.3 

2.4 

2 . 4 
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the  Hail  Project  S-band  data,  will  help  to  alleviate  the 
registration  problem. 

6.4  THE  PROBLEM  OF  RESOLUT ION- SCALE  MISMATCH 

The  work  of  Strong  (1974)  and  Changnon  (1970) 
shows  that  significant  variation  in  hail  activity  can  exist 
at  scales  finer  than  the  size  of  a  quarter- sect ion .  This 
size  scale  is  similar  to  the  resolution  of  the  C-band  at 
a  range  of  40  to  80  kilometers. 

As  estimate  of  the  energy  density  of  hail  in  the 
beam  resolution  volume  at  these  ranges  is  averaged  over 
the  volume.  The  ground  hailfall  parameter  that  actually 
corresponds  to  the  radar  estimate  would  be  the  energy 
density  averaged  over  the  quarter- sec tion  under  the  beam, 
rather  than  the  hailpad  value.  The  energy  density  from  a 
hailpad  can  deviate  markedly  from  the  area  average,  causing 
strong  scattering  in  the  E-ER  relationship.  A  rough  il¬ 
lustration  of  the  fine-scale  variation  can  be  made  from  the 
data  in  the  above-mentioned  work  of  Strong.  In  1973  he 
operated  a  high-density  network  in  Section  6,  T42,  52,  W5, 
southwest  of  Rimbey,  Alberta,  consisting  of  a  5  x  5  hailpad 
array.  He  produced  a  hailfall  map  from  this  network  of  a 
storm  which  in  that  area  was  comparable  in  intensity  to  the 
ones  analysed  in  this  study.  The  map  showed  maxima  of 
hail  activity  (about  200Jm“2)  in  the  center  of  the  section 
and  near  the  southwestern  corner.  Hail  activity  varied 


' 
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widely  over  the  network.  In  the  northwest  quarter,  the 
energy  density  varied  from  171Jm'2  in  the  southeast  corner 
to  less  than  10Jm~2  in  the  northeast  corner. 

The  results  indicate  that  a  radar  estimate  of 
energy  density  may  be  more  useful  as  a  mean  figure  for  the 
hail  activity  over  an  extended  area,  such  as  a  quarter- 
section.  While  it  will  give  a  good  rough  estimate  of  the 
energy  density  on  a  hailpad  below  the  beam,  it  does  not  yet 
show  promise  of  being  a  highly  accurate  point  estimator. 


CHAPTER  7 


CONCLUSION 


7 . 1  SUMMARY 


A  method  of  estimating  kinetic  energy  density  of 
hailfall  at  locations  without  hailpads,  using  radar  data 
and  the  existing  hailpad  network  was  developed.  The  method 
was  tested  on  the  hailstorms  of  August  7  and  18,  1974  in 
central  Alberta. 

The  radar - hai lpad  estimation  method  was  compared 
to  methods  which  use  only  hailpad  information.  The  hailpads 
only  methods  showed  little  or  no  skill  at  extrapolating 
energy  density,  while  the  radar-hailpad  method  showed 
highly  significant  skill.  The  average  error  of  estimation 
by  the  radar-hailpad  method  was  about  50  percent  of  the 
average  error  by  other  methods.  45  out  of  60,  or  75  per¬ 
cent,  of  the  hailpads  tested  by  the  radar-hailpad  method 
were  estimated  to  within  10Jm'2  of  their  actual  energy  den¬ 
sities.  By  comparison,  using  an  arithmetic-mean  on  the 
hailpad  data  gave  only  43  percent  within  lOJm-2.  Using 
dis tance- weighted  means  of  hailpad  data  gave  only  40  per¬ 
cent  within  1 0 Jm“ 2 .  Similarly,  92  percent  of  the  hailpads 
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were  estimated  within  20Jm-2  by  the  radar-hailpad  method, 
as  compared  to  61  and  62  percent  by  the  two  hai lpads - only 
methods . 

Using  the  radar-hailpad  method,  109  square  miles 
(282km2)  were  intensively  surveyed.  Radar  estimates  of 
hailfall  energy  density  in  this  area  showed  variations  hav 
ing  dimensions  that  are  in  good  agreement  with  ground  ob¬ 
servations  of  small-scale  hail  patterns.  Local  estimates 
of  hailfall  generally  agreed  well  with  the  available 
farmers'  reports.  Hail  energy  density  was  first  contoured 
in  the  intensive  survey  region  using  farmers'  reports,  and 
assuming  a  doubling  of  energy  density  with  each  higher 
"category"  of  maximum  hail  size.  When  this  was  done,  the 
energy  density  at  the  farmers'  report  site  was  within  one 
countour  category  of  the  radar  estimate  at  the  site  or 
a  neighbouring  land  location  in  13  out  of  24  cases. 

The  hailpad  energy  versus  radar  estimate  rela¬ 
tionship  was  found  to  be  quite  different  on  the  two  storm 
days  studied  here.  This  implies  that  radar  estimates 
should  be  tailored  to  the  date's  ground-based  measurements 

The  method  of  estimating  local  hail  energy  den¬ 
sity  shows  promise  in  operational  applications,  but  some 
problems  must  be  overcome.  The  S-band  (10.6cm)  radar  data 
collected  since  1975  is  generally  more  stable  and  accurate 
ly  registered  than  the  C-band  data  of  1974.  The  longer 
wavelength  of  the  S-band  will  greatly  reduce  the  effect  of 
beam  attenuation.  Refining  the  programs  used  to  locate 
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the  hailpads  and  to  select  the  appropriate  radar  data  will 
also  improve  the  results.  For  data  from  1975  onward,  the 
output  from  existing  high-precision  programs  can  be  used 
as  input  for  secondary  programs  which  perform  these  tasks. 

7.2  RECOMMENDATIONS  FOR  FURTHER  STUDY 

Refined  radar-hai lpad  studies  could  readily  be 
incorporated  into  Alberta  Hail  Project  operations  and 
evaluations . 

A  more  sophisticated  and  accurate  method  is 
needed  to  filter  out  the  contribution  of  rain  to  the  radar 
echoes.  Information  from  the  circular  depolarization  ratio 
(CDR)  could,  for  example,  be  used,  since  the  CDR  is  well 
correlated  with  the  proportion  of  rain  to  hail  (Barge, 

1972)  .  The  contribution  of  the  hail  signal  to  the  total 
returned  signal  increases  with  the  magnitude  of  the  returned 
signal.  The  energy  estimation  formula  should  therefore  be 
tested  at  threshold  levels  greater  than  45dBZ. 

The  radar-hailpad  method  should  be  more  thoroughly 
tested.  A  large  amount  of  high-quality  radar  data  and 
thousands  of  hailpad  analyses  are  available  from  1975, 

1976  and  1977.  Use  of  this  data  base  for  testing  the 
radar-hailpad  method  can  more  fully  establish  its  reliabi¬ 
lity  and  feasibility. 

Programs  already  in  use  by  the  Hail  Project  can 
extract  the  estimated  rate  of  rainfall  from  the  radar  data 
for  each  storm.  A  similar  program,  using  the  energy- 
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reflectivity  relationship  appropriate  for  the  charac¬ 
teristics  of  the  radar,  could  be  used  to  produce  maps  of 
raw  hail  energy  density  estimates  for  the  entire  Hail 
Project  region  for  a  given  storm  day.  This  raw  estimate 
pattern  could  later  be  adjusted  when  the  hailpads  for  that 
date  are  analysed. 
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on  hailfall.  Unil  now,  no  fine-scale  quantitative  data 
was  available  to  measure  this  effect.  The  desired  effect 
of  seeding  is  to  shift  the  entire  size  spectrum  of  hail 
produced  by  a  storm  toward  smaller  sizes.  The  best  single 
parameter  to  measure  this  effect  is  the  energy  density  of 
hail  at  the  ground.  Calculations  can  be  made  of  the 
average  energy  density  of,  or  total  energy  yielded  by, 
seeded  and  unseeded  hailswaths,  or  seeded  and  unseeded 
portions  of  the  same  hailswath.  A  comparison,  based  on 
many  such  calculations,  could  possibly  determine  the  over¬ 
all  effect  of  seeding.  Seeding  effects  determined  by 
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analyses  of  hail  insurance  claims  in  townships  could  be 
verified  independently  by  determining  the  relationship 
between  the  average  hailfall  energy  density  and  the  total 
hail  insurance  claims  in  each  township.  These  methods, 
which  use  many  energy  density  estimates  over  large  areas, 
have  the  inherent  advantage  of  averaging  out  the  uncer¬ 
tainties  in  individual  energy  density  estimates. 
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SUMMARY  OF  HAIL PAD  DATA  USED  IN  THIS  STUDY 
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Multiple  hail  bursts:  more  than  one  hailpad  collected  at  a  site  from  a  single  storm. 

2Tota I  energy  density  from  multiple  hail  bursts. 

3See  Page  12  of  text.  Continued 


Appendix  A 

Summary  of  Hailpad  Data  Used  in  This  Study,  Continued 
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APPENDIX  B 


1  1 1 


LISTINGS  OF  MAJOR  PROGRAMS 


Notes : 

Page 

PLOT 

112 

Complete  listing,  except  for  some  special  routines  available 
only  under  University  of  Alberta  MTS. 


PLOTPAD  126 

Comp  I ete  I i st i ng . 


PADEN 

Main  program  only  listed.  PADEN  uses  the  same  subroutines 
as  PLOTPAD. 


PPI 


Main  program  only  listed.  PPI  uses  all  the  subroutines  as 
PLOT,  except  SUMARY,  LAND  and  CREATF. 


129 


131 


PERMS 


Comp  I ete  I i st i ng . 


■ 


112 


PLOT 


INTEGER  QRT ? SEC  , TNSHP , RANGE y RRNGE ( 50 ) , RAZM ( 50 ) , SPACE ( 50 ) ,11 D ( 4 ) * 
♦LTIMrHTIMr ID(4r50)  y 

#STLTI M r ST NT  I M , SPLT I M , SPHT I M , BLANK , T IMSTT ( 3 ) ? T IMSTP ( 3 ) , HR , M IN , SEC , 

* C NTNUM  y  BIN  y  HTIM  ? L T I M 
INTEGER #2  M  < 3072 ) , LEN ? IHTIM , ILTIM 
LOGICAL  DOWN r UP 
COMMON  ID r I NTEN , NUMBER 
COMMON  II 
LOGICAL#!  A ( 6 ) 

DATA  BLANK/ '' 

R  E  A  L  P  A  Z  M  (  5  0  >  y  I N  T  E  N  (  5  3  )  *  I N  C  E  P  T 

R  T I M  ( I H  T I M  y  I L.  T I M )  =  I H  T I M  *  3 2  7  6  8  + 1 L.  T I M 

B I N T I M  <  H R ?  MIN? S E C ) = H R * 2 1 6 0 0 0  *  + M I N * 3 6 0 0  *  f SEC# 6 0 * 

WRITE (6, 300) 

3  0  0  F  0  R  M  A  T  (  '  1  I N  P  U  T  B  E  G I N  N I N  G  T I M  E  A  N  D  E  N  D I N  G  T I M  E  0  F  S  U  R  V  E  Y  '  *  /  ) 

R  E  A  D  (  5 , 1 0  0  )  T I M  S  T  T  »  T I M  S  T  P 

100  FORMAT  ( 6  <  12?  :I.X)  ) 

WRITE ( 6  y  305 ) 

3 0 5  F 0  R  MAT  < '  ' ? / / ,  '  I N P  U T  S L 0  P E , I N T  E R  C  E  P T y  A  N  D  R  A D A R  C  0  N  S  T  A  NT' ?/> 

C  ERE  AD  IS  A  FORTRAN  SPECIAL  ROUTINE  USED  UNDER  MTS  y  IJ  OF  A » 

C  A  L  L  F  R  E  A  D  (  5  y  '  3  R  $  '  ?  S  L  0  P  E  y  I N  C  E  P  T  ?  R  A  D  A  R ) 

S T  A R T = B I N  T 1 M ( T IMS  T T ( 1 ) y T I M S T  T ( 2 ) ? T I M  S T  T ( 3 ) ) 

T STOP=B I NT  I M ( T IMSTP (l)r T I MSTP  <  2 ) r  T I MSTP ( 3 ) ) 

2000  WRITE ( 6, 301 ) 

3  0 1  F  0  R  M  A  T  (  '  '  v  /  /  r  '  I N  P  U  T  0  U  A  R  T  E  R  r  S  E  C  T 1 0  N  ,  T  0  W  N  S  H I P  ,  R  A  N  G  E  ,  '  / '  M  E  R I D I A  N  r 

*  A  N D  I D  E N  T I F I C  A  T 1 0  N  F 0  R  E  A C  H  L 0  C  A  T I ON* ' > 

NUMBER- 50 

DO  1  1=1 r 50 

8  R  E  A  D  <  5  y  1 0 1 )  0  R  T  y  S  E  C  ?  "f  N  S  H  P  r  R  A  N  G  E  r  M  E  R  ?  < I D  (  J  r  I )  y  ,.J  =  1  r  4  ) 

I F  >:  a  R  T  *  E  Q  *  B  L  A  N  K  )  G  0  T  0  9 

1 0 1  F  0  R  M  A  T  (  A  2  ?  4 1 .1 0  y  4  A  2  ) 

C  CONVERTS  SURVEY  COORDINATES  TO  RADAR  RANGE  (KM)  AND  AZM 
C  <631  PER  SCAN)  WITH  ZERO  AT  140  DEGREES  FROM  TRUE  NORTH* 

5  C  A  L  L  L  A  N  D  <  Q  R  T  ,  S  E  C  r  T  N  S  H  P  ?  R  A  N  G  E ?  M  E  R  ?  R  R  N  G  E  <  I )  ?  R  A  Z  M  ( 1  )  ?  P  A  Z  M  ( I  )  ) 

C  C  H  E  C  K  F  0  R  I N  V  A  L I D  I...  A  N  D  L  0  C  A  T 1 0  N  S  .  I F  S  0  ?  I N  F  0  R  M  T  H  E  W  0  R  L  D  . 

I F  C  R  R  N  G  E  ( I )  *  G  T  .  1 4  7  )  G  0  T  0  6 
I F  (  R  R  N  G  E  (  1 )  *  G  T  *  3 )  G  0  T  0  1 
UR  I TE  (  6  y  302  )  (  I D  <  J  r  I  )  r  -J=  1  r  4  ) 

3  0  2  F  0  R  M  A  T  <  '  P  A  D  '  ?  4  A  2  y  '  I S  L  E  S  S  T  H  A  N  3  K  M  *  F  R  0  M  T  H  L  R  A  D  A  R .  N  0  D  A  T  A 

*.  I S  A  V  A I L  A  B  I...  E  .  T  H  E  L  0  C  A  T 1 0  N  I S  D  E  L  E  T  E  D  '  ) 

GO  TO  8 

6  W  R I  V  E  (  6  y  3  0  3  )  (  I D  (  J  r  I  )  r  J  =  1  ?  4  ) 

3  0 3  F  0 R M  A T ( '  P A  D  ' , 4  A  2  ? '  I S  M 0 R  E  T H  A  N  1 4  7  K  M  *  F  ROM  T  H  E  R  A D  A  R *  N  0  D  A  T  A 

*  I S  A  V  A I L  A  B  L  E  *  T  H  E  L..  0  C  A  T 1 0  N  I S  D  E  L  E  T  E  D  .  '  ) 

GO  TO  8 

1  CONTINUE 

GO  TO  12 

9  NUMBER=I-1 

C  BUBBLE  SORT  LOCATIONS  TO  PREVENT  COMMON  REFERENCE  OF  SOURCE, 

12  N1 --NUMBER' --l 
DO  2  1=1 r Ml 
N2=N1-I +1 
DO  3  J  =  1 y  N 2 

I F  (  R  A  Z  M  (  J  )  *  L  "f  ,  R  A  Z  M  ( J  + 1 )  >  G  0  T  0  3 

IF  (RAZM(J) »  E 0 . R AZMCJ+l) . AND  * RRNGE ( J ) * LT . RRNGE ( J+l ) )G0  TO  3 
DO  15  K  = 1 y  4 
1 1 D  (  K  )  =  I D  (  K  y  J  > 

15  CONTINUE 

IAZM--RAZM  ( J ) 

R  P  AZM  ■-  P  A  Z  M  (  J  ) 

IR=RRNGE ( J ) 
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DO  16  K  ~  1  *  -4 

ID  <  K  f  J ) “ID ( K  p  J+l ) 

16  CONTINUE 

RAZM ( J ) “RAZM ( J  +  l  ) 

PAZM  (  J ) “PAZM  <  J+ 1 ) 

RRNGE ( J  >  =RRNGE ( J+l ) 

DO  17  K=1  y  4 
ID (K  y  J+l )  =  1  ID ( K ) 

17  CONTINUE 
RAZM< J+l )“IAZM 
PAZM (J+l )~RPAZM 

C  CALCULATE  NUMBER  OF'  RAYS  TO  SKIP  TO  GO  FROM  ONE  LOCATION  TO  THE  NEXT  * 
RRNGE ( J+l )=IR 

3  CONTINUE 

2  CONTINUE 

I SUM55!) 

WRI TE ( 6  y  400  )  NUMBER  y  (  ( I D  (  K  y  I )  yK=l  y4)  yF'AZMd  )  y  RRNGE  ( I  )  y  1  =  1  yNUMBER) 
400  F'ORMAT  (  '  PAD  IDENT  AZM  D  RNGE  KM'»I6/('  '  y  4A2  y  2X  y  F6  .  0  y  16  )  ) 

SPACE ( 1 ) =RAZM ( 1 ) 

DO  4  1=2  y  NUMBER 

SPACE ( I >  =RAZM ( I ) -RAZM (1-1 ) 

4  CONTINUE 

C  READ  ONE  RECORD y ALSO  INITIALIZE  SUM  IN  SUMARY » 

WRITE (6 y 402) 

4 0 2  F ORMAT  (  '  '  y  '  E N T E R  T A P E  M 0 IJ N T  I '  y  / ) 

2 2 2  C A L I-  F READ  (  5  y  '  I  I  S  1 1 ) 

I F  (II.GT.O. AND .II. LE .13) GO  TO  224 
UR  I TE ( 6  y  403 ) 

4  0  3  F  0  R  MAT  (  '  '  i  '  1 N  0  A  L I D  I N  P  U  T  .  .  P  L.  S  R  E  -  ENT  E  R  I N  P  LJ  T  '  / ) 

GO  TO  222 

224  CALL  TPMNT(II) 

C  CREATE  OUTPUT  FILE  ON  WHICH  DATA  WILL  BE  WRITTEN. 

CALL  CREATE 

U  R I T  E  (  7  )  N  IJ  M  B  E  R  y  1 D  y  P  A  Z  M  y  R  R  N  G  E 

C  SKIP  A  NUMBER  OF  RECORDS  (6144  IN  LENGTH)  DEPENDING  ON  INPUT  TIME. 

C  A  L  L  S  K  F'  R  E  C  ( 1 1  .*  S  T  A  R  T  y  M  y  C  N  T  N  U  M  ) 

C  GET  SOME  MORE  DATA  ON  TAPE. 

CALL  DATA ( M y  CNTNUM  y  840 , 82000 ) 

C  INl'TILAZE  SLJMARY. 

CALL  SUM ARY ( 1 ) 

DOWN55 .  TRUE  ♦ 

UP==.  FALSE. 

C  SPACE  FORWARD  TO  START  TIME. 

C  GET  DOWN  PULSE  (HEADING  MARKER55 12 )  . 

1 1  CALL  HED I NG ( M  y  CNTNUM  y DOWN  y  UP  y  8 1 1 y  840  y 82000 ) 

C  EVERY  BYTE  OF  DATA  HAS  TO  BE  SWITCHED  AROUND  ( PDF'!  1/45  y  DACZ )  . 

CALL  BYTREV ( M  y  CNTNUM+1 ) 

CALL  BYTREV  (  M  y  CNTNLJM  +  2 ) 

C  CONVERT  FROM  SYSTEM  TIME  TO  REAL  TIME. 

CLOCK^RTIM ( M ( CNTNUM+2 ) y M (CNTNUM+1 ) ) 

IF ( CLOCK . GE . START ) GO  TO  14 
CNTNUM55  CNTNUM +2 
GO  TO  11 

C  SET  ELEVATION  COUNTER  TO  ZERO  (HEADING  MARKER--12  FOLOLOUED  BY 
C  HEADING  MARKERS). 

C  SKIP  TO  NEXT  HEADING  MARKER.  IT  WILL.  INDICATE  A  ZERO  DEGREE  SCAN, 

14  CALL  HED  ING  (  M  y  CNTNUM  y  .FALSE,  y  .TRUE,  y  8 1 4  >■  840  y  82000) 

CALL  BYTREV ( M  y  CNTNUM  +  1 ) 

C  A  L..  L  B  Y  T  R  E  V  (  M  y  C  N  T  N  LJ  M  +  2  ) 

CL0CK55RTlM(M(CNTNUM+2)  y  M  ( CNTNUM+1 )  ) 

HTIM=M( CNTNUM+2) 

L.TIM-~M  ( CNTNUM+1 ) 

CNTNUM55CNTNUM+2 
22  ELEV=0 . 

CALL  TIME ( CLOCK y HR y MIN y SEC) 

C  BTD  IS  SPECIAL  MTS  ROUTINE ,( CONVERSION  OF  CHARACTERS) 
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CALL.  BID  ( HR  7  A ( 1 ) y  2  y  ND  7 ' 0 ' ) 

CALL  BTD ( MIN  y  A ( 3 )  7  2  7  ND 7 ' 0 ' ) 

CALL  BTD (SEC  7 A<5) y2yND7 'O' ) 

WRITE ( 6  7  401 ) A 

401  FORMAT ( '  ' y 2A1 > ' t ' y 2A1 7 ' l ' y 2A 1 ) 

IF  (CLOCK . GT . TSTOP ) GO  TO  40 
23  INTEN ( NUMBER+1 ) ~HTIM 
INTEN ( NUMBER  ! 2 )  -LTIM 
DO  20  1  =  1 7  NUMBER 

C  SUBTRACT  3  BECAUSE  FIRST  THREE  BINS  ARE  NOT  RECORDED 
B I N-RRNGE ( I ) -3 
N=SPACE( I ) 

IF  (N.EQ.O)GO  TO  25 

CALL  SKPRAY  (  M  y  CNTNUM  y  N  y  840  y  82000  ) 

CNTNUM=CNTNUM+2 
I B I N  ”  0 

25  CALL  DBZRAY ( M  y  CNTNUM  y  B IN  y IB  IN  y  R  y  SLOPE  y INCEPT  y  RADAR  y 840  y  82000 ) 

26  INTEN ( I ) =R 

20  CONTINUE 

INTEN  (  NUMB E R  +  3 )  =  EL.EU 
WRITEC7) INTEN 
CALL  SUMARY (2) 

CALL  HED I  NO  (  M  y  CN  I'NIJM  y  ♦  FALSE  .  ?  *  TRUE  ♦  y  827  y  840  y  82000  > 

CNTNUM=CNTNUM+2 

IF  ( ELEU  *  EQ  » 22 ♦ 5 ) GO  TO  115 

ELEV=ELEV+1 .5 

GO  TO  23 

115  CALL  HEDINGCMv CNTNUM y .FALSE. y .TRUE. , 827 y 840 y 82000 ) 

GO  TO  115 

27  ISUM=ISUM+1 

CALL  BYTREV ( M  y CNTNUM !2 ) 

CALL  BYTREU ( M  y CN TNUM+1 ) 

HT I M~M ( CNTNUM+2 ) 

LT I M~M ( CNTNUM ! 1 ) 

CLOCK=RTIM(M< CNTNUM+2) y  M ( CNTNUM+1 ) ) 

IF  ( ISUM.EQ. 17) GO  TO  41 
GO  TO  22 

40  STOP 

41  CALL  SUM ARY (3) 

ISUM=0 

GO  TO  22 
END 
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C  A  SUBROUTINE  TO  BYTE  REVERSE  INTEGERS  DATA 


SUBROUTINE  BYTREV  <  L  r  I ) 
LOGICAL*!  1.(6144)  y  LA 
L  A-L < ( 1-1 ) *2+1 ) 

L( ( I- 1 >  *2+1 ) -L ( 1*2 ) 

L< 1*2) -LA 

RETURN 

END 
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r;  A  SUBROUTINE  TO  DETECT  RAYS  AND  .SKIP  TO  THE  BEGINNING 
C  OF  A  DESIRED  RAY ♦  FORWARD  SKIPS  ONLY* 


S  U  B  R  0  U  T I N  E  S  K  P  R  A  Y  <  M  ,  C  N  T  N  U  M  >  N  y  *  »  *  ) 

INTEGER  CNTNLJM  ■>  N 
IN T EGER# 2  M<3072) 

1.2  IF  (  N  *  LT  ♦  0  )  GO  TO  10 

3  I F  (  N  *  E  Q  ♦  0  )  R  E  T  l  J  R  N 

:L  CNTNUNsCNTNUM+1 

I F  <  C  N  T  N  U  M  *  G  T  *  3  0  7  2  )  C  A  I...  I ..  D  A  T  A  (  M  y  CN  T  N  U  M  ?  %  4  r  &  5  ) 

I F  <  H  (  C  N  T  N  U  M  )  *  E  G  *  1 4  4  )  C  A I...  L_  D  A  T  A  (  M  »  C  N  T  N  U  M  t  &  4  ,  &  5  > 
I F  <  M  <  C  N  T  N  U  M  )  *  E  Q  *  0  )  G  0  T  0  2 
GO  TO  1 
2  N=N-1 

GO  TO  3 

4  RETURN 1 

5  RET URN 2 

10  CNTNUM-CNTNUM-l 

I F  (  M  (  C  N  T  N  U  M  )  *  E  Q  *  0  )  G  0  1“  0  1 1 
GO  TO  10 

11  N®N+1 

GO  TO  12 
END 


C  A  SUBROUTINE  TO  READ  DATA  FROM  TAPE  AND 
C  MODIFY  CNTNUM  (COUNTER  NUMBER)  BACK  TO  1 


SUBROUTINE  DATA  ( M i CNTNUM »*,* ) 
INTEGER  CNTNUM  f  CONTRL ( 20 ) 

INTEGER*?.  M  (  3072 )  »  LEN 

MOD-O 

LUN I T“0 

1  C A L L  R E A D ( M » L E N » M 0 D»LNR»LUNIT»S10) 
CNTNUM" I 

RETURN 

10  CALL  MOUT'P  (  1 1 1  ) 

GO  TO  (4?3r2> t III 

2  RETURN 1 

3  GO  TO  1 

4  RETURN? 

END 


SUBROUTINE  TO  DETECT  HEADING  MARKERS  AND  DOWN  PULSES 


SUBROUT I NE  HED I NG  ( M , CNTNUM > DOWN , UP y  *  y  *  ,  *  ) 
INTEGER  CNTNUM 
INTEGER *2  M(3072) 

LOGICAL  DOWN » UP 
CNTNUM-CNTNUM+l 

IF < CNTNUM. GT. 3072) CALL  DATA < M » CNTNUM »  &2 »  S4 > 

IF <M( CNTNUM) .EQ. 144) CALL  DATA < M » CNTNUM >  *2*  &4 ) 

I F ( M < CNTNUM ) ♦ EQ ♦ 0 ) GO  TO  3 

IF ( M  < CNTNUM > ♦ EQ . 12 . AND . DOWN ) RETURN 

IF ( M ( CNTNUM ) . EQ . 8 ♦ AND  *  UP ) RETURN 

I F ( M ( CNTNUM ) . EQ . 1 2 . AND . UP ) RETURN 1 

GO  TO  1 

C  N  T  N IJ  M  •-  C  N  T  N  U  M  +  2 
GO  TO  1 
RE TURN 2 
RE TURN 3 
END 
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C  SUBROUTINE  TO  GET  A  SPECIFIED  BIN  OF  DATA  FROM  CURRENT  RAY 

C 

C  B I N-RRNGE ( I ) -3 
C 

C  ~1 6385~>ERRGR  FLAG 
C  -28672=>END  OF  RECORD 
C  204 EL- >  NEW  MARKER 

C  3072=>  DOWN  PULSE 


SUBROUTINE  DBZRAY ( M  »  CNTNUM ?  BINr IBIN»R» SLOPE r INCEPT  i RADAR  r  Y  >  Y  ) 
REAL  INCEPT 
INTEGER  CNTNIJM  t  BIN 
INTEGERY2  M<3072> 

10  CNTNUM-CNTNUM+l 

I F ( CN TNUM *  GT ♦ 3072 ) GO  TO  1 5 
1 3  CALL  BYTREV  (  M  »  CNTNIJM  ) 

I F  (  M  (  CNTNIJM  )  *  GE  ♦  0  )  GO  TO  11 
IF  (  M  <  CNTNUM  )  ♦  EQ  ♦  -16385  )  GO  TO  9 
IF ( M <  CNTNUM ) ♦ EQ ♦ -28672 ) GO  TO  15 
C  SURELY  IT  IS  A  POINTER 
IB I N=M< CNTNUM) +145 
I F  < I BIN* GT ♦ B I N ) GO  TO  91 
GO  TO  10 

C  ERROR  IN  NEXT  DATA  WORD ♦  SKIP  THE  WORD ♦ 

9  CNTNUM=GNTNlJM  +  2 

GO  TO  13 

1 5  CALL  DATA ( M , CNTNUM  ? X 1 00  ? X 1 0 1 ) 

GO  TO  13 

C  NUMBER  IS  POSITIVE  SO  COULD  BE  FLAG » TIME  OR  DATA* 

1 1  I F  <  M  <  C  N  T  N  U  M  )  ♦  E  Q  ♦  0  )  G  0  T  0  1 7 

IF ( M ( CNTNUM >  *  EQ  *  2048 >  GO  TO  89 
I F  (  M  (  C N  )“ N tJM  )  ♦  E Q  *  3 0 7 2  )  G 0  T 0  8 8 
I F ( BIN* E Q  *  I B I N  > G  0  T 0  9  0 
IBIN=IBIN+ 1 
GO  TO  10 

C  NEW  PULSE* THEREFORE  M(BIN)  EQUALS  ZERO* 

17  R~0  * 

CNTNUM'-CNTNUM-l 

RETURN 

C  DOWN  PULSE t  WHAT  DO  WE  DO??? 

8 8  C  A L L  B  Y T R E  V  <  M  »  C NTNLJM  ) 

GO  TO  17  r 

C  NEW  marker:  WHAT  DO  WE  DO??? 

8 9  C  A L L  B  Y T R E V ( M  r  CNTNUM ) 

GO  TO  17 

C  CALCULATE  RANGE  IN  KM  ALSO  DBM  AND  DBZ ♦ 

90  R'RANGE”  ( IBINY1  ♦  048  )  +2  »  67 
DBM^SLOPEYM <  CNTNUM )  +  INCEPT 
R-DBM+ ( 20YAL0G 1 0 ( RRANGE ) ) + RADAR 
RETURN 

9 1  C  A  L  L-  B  Y  T  R  E  V  (M?  CNTNUM) 

GO  TO  17 

100  RETURN 1 
.101  RETURN2 
END 
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C  THIS  SUBROUTINE  CREATES  THE  HOURLY  DATA  SUMMARY 


SUBROUTINE  SUM  ARY  (ID 
REAL  INTEN ( 53 ) 

I NT EGER  SUM  <  50 , 1 0 ) , ID ( 4 , 50  > 

COMMON  ID  , INTEN  »  NUMBER 
GO  TO  <  1000  ,  2000  *  3000  )  ,  1 1 

1000  DO  1001  J=1 , 10 
DO  1002  1=1 >50 
SUM ( I  *  J ) =  0 

1002  CONTINUE 

1001  CONTINUE 
RETURN 

2 0 0 0  D 0  2 0 0 1  I==l  y  NUMBE R 

I DB- I F I X ( INTEN < I ) / 1 0 ♦ )  + 1 
IF ( IDB.GT* 10) IDB=10 
I F ( IDB.LT,0)IDB=0 
SUM ( I , I DB ) “SUM ( I t I DB ) + 1 

2001  CONTINUE 
RETURN 

3000  WRITE ( 6  , 3001 ) 

3001  FORMAT ( ' 1SUMARY  OF  DATA  FOR  THE  LAST  PERIOD'/ 

% '  NUMBER  OF  TIMES  RADAR  REFLECTIVITY  WAS  IN  RANGE  OF'/ 
#'  PAD  I DENT  0-9  10-19  20-29  30-39  40-49  50-59  60-69' 
70-79  80-89  90-99  DBZD 
DO  3002  1=1, NUMBER 

WRITE (6, 3003) ( ID( J,I ) , J  =  1 ,4) , <  SUM ( I , J ) , J=1 , 10) 

3003  FORMAT  < '  ' , 4 A2 , 1 0 1 6 ) 

3002  CONTINUE 
GO  TO  1000 
END 
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A  SUBROUTINE  TO  CONVERT  REAL  TIME  INTO  HOURS  AND  MINUTES 


SUBROUT I NE  T I ME ( CLOCK  > HR  * M I N t SEC ) 
INTEGER  HR t SEC 
H  R  =  I N  T  ( C  L  0  C  K  /  2 1 6  0  0  0  ♦  ) 

MIN  =  INT (CLOCK/36 0 0 . ) - H R  *  6  0 

SEC- I NT ( CLOCK/ 60 ) - ( HR#3600+MIN*40 ) 

RETURN 

END 


•O  «T 
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c  the:  subroutine  to  create  the  file  on  the  output  tape. 


1  0 1 


104 


10 

103 


11 

102 


SUBROUTINE  CREATE 
LOGICAL.*!  F I L  (10) 

LOG I CAL* 1  CHAR <  99 )/'$'»' M 

'•  O  '  . .  '  T  t  .  '  p  '  .  / 


9  ' 0 ' 9  ' U '  9  ' N'  9 ' 7 ' 9  ' 


'0' 


*r'0' t'  ' n ' 9 ' y ' T ' r ' P ' r '  '  9  9  ' 7 '  9  ' A'  9  ' P ' 


v  — '  y 
*■?'  Q'  >  '  L'  y 

*  '  9  'L'  y  ' R 


E 


r 

0  '  t  '  U  '  t 
r  '  C  '  f  'L. 


r 

T  '  y  '  P  ' 
_  / 


9  '*'  9  'E'  9 ' 0 
'  9  '  F  '  9  '  M  '  9 


'I 


*  y  '  1  '  9 


9-0' 9 ' 0 


/ 


y 

CALL  FTNCMD  (  '  ASSIGN  4~NELJD  I  RECTORY  '  y  21) 
CALL  FTNCMD  <  '  ASSIGN  3=NELJD  I  RECTORY  (  *L  +  1  ) 
WRITE (6 > 101 ) 


27) 


9 '0' 9 '3'  9  'A'  9 '7' 
'y'R'y'I'y'N'y' 
y  ' T '  y  ' %'  9  '  '  y'O' 

7 ' 9  '  ~  '  9  '  F  '  9  '  B  '  9  ' 
9  '  2' 9  '  E  '  9  '  -  '  9 ' 2' 


FORMAT  (  '  PL..S  ENTER  NAME  OF  FILE  TO  BE  CREATED  ON  TAPE'  »/) 

CALL  FREAD  < 5 9  ' STR I NG l '  9  CHAR ( 88 ) , 1 1 ) 

R  E  A  D  <  4  9 1 0  4  r  E  N  D  =10)  F I L. 


FORMAT ( 10A1 ) 

I F  (  L.COMC  <  1 0  »  F I L  >  CHAR  <  88  )  )  ♦  EG  ♦  0  )  GO  TO  11 
GO  TO  4 


UR I TE <  3  » 1 03  >  < CHAR ( J ) , J=88  y  98 ) 


FORMAT  <  1  J.  A 1  ) 

CALL  CMD ( CHAR  9 98 ) 

CALL  FTNCMD < ' ASSIGN  ?=*TAPE* ' t 1 5 ) 
RETURN 


WR I TE <  6 , 102 ) 

FORMAT ( '  ' 9  '  FILE  ALREADY  EXISTS  ON  TAPE...PLS  RE-ENTER  FILE  NAME. 
*  '  „  /  ) 

GO  TO  6 
END 


> 
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C  THIS  SUBROUTINE  SKIPS  A  NUMBER  OF  RECORDS  DEPENDING 
C  ON  THE  POSITION  OF  MARKER  ON  TAPE. 


SUBROUTINE  SKPREC (II* START  * M * CNTNUM ) 

INTEGER  HR  *  MI N  * SEC  * HR 1  *  MINI 
I  NT  EGER  *2  M(3072)  f  IHTIMr  IL.TIM 
INTEGER  CNTNUM 
LOGICAL  DOWN * UP 

INTEGER *2  HRS ( 1 3 ) /16 * 18 * 13 * 15*17*19*20*21 *  14*15*16* 17*18/ 
INTEGERS  MINU ( 13 ) /OS  *  32  *  10  *  55  *  27  *  40  *  18  *  55  *  43  *  35  *  55  *  37  *  36/ 
R  T I M  <  I H  T I M  *  I L  T I M  >  ~  I H  T I M  *  3  2  7  6  8  ■ + 1 L  T I M 
CLOCK-START 

C  A  L  L  T I M  E  <  C I...  0  C  K  *  H  R  *  M I N  *  SI  E  C  ) 

T0TAL“<HR*60  +  MIN>--<HRS(II  >*60  +  MINU<  II )  ) 

9  I F  (  TOTAL  .  LE  .  1 2  .  AND  .  TOTAL. .  GE  .  0  )  RETURN 
CALL  SK I P  (  0  *  I  NT  (  (  TOT  AL../6  )  #75  )  *  0  *  X4  ) 

CALL  DAT  A ( M  *  CNTNUM ) 

1 0  CALL  HED I NG ( M  *  CNTNUM * ♦ TRUE ♦ * .FALSE. *  £10) 

CALL  BYTREU ( M  *  CNTNUM+ 1 ) 

CALL  BYTREV ( M  *  CNTNUM+2 ) 

CLOCK-RTI M  (  M  (  CNTNUM  +  2  )  *M(CNTNUM+.l )  ) 

CALL  T 1 ME ( CLOCK  *  HR 1  *  M I N 1  *  SEC ) 

T  0  T  A  L  =  (  H  R  *  6  Of  MI  N  )  -  <  H  R 1  *  6  0 i  M I N 1 ) 

GO  TO  9 
4  RETURN 

END 
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C  THIS  SUBRGUT I  NEi  CHOOSES  A  RADAR  DATA  TAPE  FROM  A  NUMBERED  SERIES 
C  AND  MOUNTS  IT. 


SUBROUTINE  TPMNT(II) 

INTEGER*2  LEN/24/ 7 LEN1/80/ 

INTEGER  CONTR ( 20 ) 

GO  TO  (Ir2y3y24y25y26y27y8y9yl0yllyl2yl3)yll 
WRITE <6y 101 ) 

101  FORMAT ( '  PLS  ENTER  TAPE  MOUNT '  y  / ) 

READ  <5*1 02 ) CONTE 

102  FORMAT (80A1) 

C  MOUNT  IS  ANOTHER  SPECIAL  MTS  FORTRAN  ROUTINE. 


CALL  MOUNT ( CONTR y LEN 1 ) 
GO  TO  6 

1 

CALL  MOUNT ( '  003271 
GO  TO  6 

9  T  P 

*T* 

SIZE  =  6144 ' y  LEN ) 

0 

CALL  MOUNT ( ' 003272 
GO  TO  6 

9  TP 

*T* 

SIZE  -6144 ' r  LEN ) 

3 

CALL  MOUNT ( ' 003260 
GO  TO  6 

9  TP 

*T* 

SIZE~6144' y  LEN ) 

24 

CALL  MOUNT ( '003261 
GO  TO  6 

9  TP 

*T* 

SIZE--6144  '  y  LEN  ) 

25 

CALL  MOUNT ( '003262 
GO  TO  6 

9  TP 

*T* 

SIZE  -6144 ' y  LEN ) 

26 

CALL  MOUNT < '003263 
GO  TO  6 

9  TP 

*T* 

SIZE-6144' y  LEN ) 

27 

CALL  MOUNT ( '003264 
GO  TO  6 

9  TP 

*T* 

SIZE-6144 ' y  LEN ) 

8 

CALL  MOUNT ( '003265 
GO  TO  6 

9  TP 

*T* 

SIZE-6144' y  LEN ) 

9 

CALL  MOUNT ( '003266 
GO  TO  6 

9  TP 

*T* 

SIZE-6144' y  LEN ) 

10 

CALL  MOUNT ( '003267 
GO  TO  6 

9  TP 

•*T* 

SIZE-6144' y  LEN ) 

11 

CALL  MOUNT ( '003268 
GO  T  0  6 

9  TP 

*T* 

SIZE-6144' y  LEN) 

12 

CALL  MOUNT < ' 003269 
GO  TO  6 

9  TP 

*T)K 

SIZE-6144' y  LEN ) 

13 

CALL  MOUNT ( '003270 

9  TP 

)XT* 

SIZE-6144' y LEN) 

6 

CALL  FTNCMDC 'ASSIGN  0- 
RETURN 

END 

*T*' 

y  12) 

'  l 


1 
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C  THIS  SUBROUTINE::  CHANGES  THE  DATA  TAPE  DURING  PRINT-OUT  OF  THE  DATA 
C  SUMMARY ♦ 


SUBROUTINE  MOUTPdII ) 

REAL  INTEN(50) 

INTEGER  I D (4*50) 

COMMON  ID 1 1NTEN r  NUMBER 
INTEGER*2  L.EN 
COMMON  II 
CALL  SUMARY (3) 

WRITE < 6 t 100) 

100  FORMAT ( '  '  t  '  END-OF-FILE  ON  TAPE# ' * // r '  WHAT  DO  YOU  WANT  TO  DO 

1UITH  TAPE? ' * // 9  '  ENTER  NUMBER  CORRESPONDING  TO  THE  FOLLOWING J ' * // r 
*'  1 =>REW IND  AND  RE-USE. '»/*'  2=>UNL0AD  AND  MOUNT  NEW  TAPE. 7 » 

%/■>'  3  =  > U N L 0  AD  AN D  S T 0 P  *  '  t  /  ) 

4  CAL  L  F  R E  ADOy'I.'y 1 1 1 ) 

GO  TO  (  5  y  6  y  7  ) y  1 1 1 
WRITE (6? 102) 

102  FORMAT ( '  INVALID  INPUT ♦ * ♦ PLS  RE-ENTER  INPUT'?/) 

GO  TO  4 

5  CALL  REWIND (0) 

RETURN 

6  W R I  T E ( 6  r 103) 

1 0 3  F 0 R  M  A  T <  '  P  L  S  E N T E R  N  E W  T APE#') 

19  CALL  FREAD (Sy'I.'yll) 

I F ( 1 1 . GT . 0 . AND . 1 1  *  LE  *13) GO  TO  600 
WRITE (6? 102) 

GO  TO  19 

600  CALL  TPMNT(II) 

RETURN 

7  RETURN 
END 
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PLOTPAD 


REAL  PAZMC50)  r  INTEN<53>  y  ELEU  (.15)  ,  A  ( 2 )  y HEIGHT  ( 15  ) 

REAL  IHTIM» ILTIM»TIM1 »TIM2»CL0CK 

INTEGER  RRNGE  ( 50  )  ,  H  1ST  <1 5 , 20 )  ,A1<20>  » ID(4y50)  y WR0NG(2) 
*yPAD<4) »HR 
DATA  WRONG/ ' ST ' , ' OP ' / 

LOGICAL  EQUC 
LOGICAL*!  FIL(IO) 

LOGICAL*!  A2 ( 6 ) 


700 


1112 

1109 
407 

.1111 

401 

1 1 1 0 


RTIM(  IHTIMy  ILTIM)=sIHTIM*32768  +  ILTIM 
LOGICAL*!  CHAR ( 28 ) / ' $ ' , ' C ' y ' 0 ' y 7  N ' y '  '  y  '  *  '  , 'T ' , ' A 
*  ,  '  P  '  t  '  0  '  ,  '  S  '  t  '  N  '  y  '  =  '  / 

WRITE (6» 700) 

FORMAT ( '  'r'  PLS  INPUT  NAME  OF  OUTPUT  FILED/) 
CALL  FTNCMIK  ' ASSIGN  4-NEWD I RECTORY x  >26) 

CALL  FREAD <5* ' STRING : ' » CHAR  < 18) » 1 1 ) 

REWIND  4 


READ  <  4 , 407  v  END”  1 1 .1 1  )  F I L 
FORMAT ( 10A1 ) 

I F  ( I..COMC  <  1 0  y  F I L  r  CHAR  ( 1 9  )  )  *  EQ  .  0  )  GO  TO  111 0 
GO  TO  1.109 

WRITE:"  (  6  »  401 )  (CHARE  I)  »I  =  18»28) 

FORMAT  (  '  '  y  FILE  '  y  1 1A1  y  "DOES  NOT  EXIST. .PLS  RE-ENTER  FILENAME  '  >  /  ) 

GO  TO  1112 

CALL  F T N C M D ( C H A R y  2 8 ) 

C A L L  F T N C MD ( ' $ C 0 P Y  * T A P E *  - Y ' » 1 5 ) 

C A L L  F T N C MD ( 'ASS I G N  7--Y ' , 1 1 ) 


C  A  L  L  F  T  N  C  M  D  (  '  $  R  E  L  E  A  S  E  *  T  A  P  E  *  '  » 1 5  ) 

6  0  5  R  E  A  D  <  7  y  E  N  D  -  6  0  0  )  N  U  M  B  E  R  ,  I D  y  P  A  Z  M ,  R  R  N  G  E 
C  READ  ID  OF  DESIRED  PAD. 


WRITE (Ay  200) 

200  FORMAT ( '  INPUT  DESIRED  PAD  IDENTIFICATION  OR  “STOP" V) 
3  READ <6y 100) PAD 

100  FORMAT ( 4A2 ) 

IF  (  PAD  ( 1  )  .  EO  .  WRONG  ( .1 )  .  AND  .  PAD  (  2 )  .  EQ  .  WRONG <  2 )  )  STOP 
604  DO  1  IP  - 1  y  NUMBER' 

DO  SO  K - 1 y  4 

I F  <  P  A  D  ( K  )  .  N  E  .  I D  (  K  >  I P  )  >  G  0  T  0  1 


SO  CONTINUE 
GO  TO  2 

1  CONTINUE 

WR I TE <6*201 >  PAD 

201  FORMAT <  '  PAD  /y4A2y/  CANNOT  BE  FOUND.  PLS  RE-ENTER  I DENT'/) 
GO  TO  3 

2  N1 -NUMBER-!  1 
N 2 -NUMBER! 2 
N3- NUMBER +3 
DO  10  1=1  y 15 

HE I GHT ( 16- I ) = ( RRNGE ( I P ) ) *1 . 5* ( I- 1 ) *3 . 1 4 15927/ ISO . 

10  CONTINUE 

DO  1.1  J=1  y  20 
DO  12  K=1 y 15 
HIST  (  K  y  -J )  -0 
12  CONTINUE 

11  CONTINUE 

DO  13  ,J “  1  y  2 0 


A1 ( J ) -0 

13  CONTINUE 

6  R  E  A  D  (  7  y  E  N  D = 6  0  DINT  E  N 

14  TIM2=.INTEN  <N1 ) 

T I  Ml “IN TEN (N2 ) 

Ml  IDO 


11  =  0 
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A 1  (  1  )=0 

7  CALL  EX  TRAC  (  HIST  y  1  y  ELFJU  y  INTEN  y  N3  y  IP  y  LIH  y  1 1 ) 

I F ( L I H . GT ♦  M I H  )  M I H-L I H 

IF ( I I ♦ EG ♦ ~ 1 ) GO  TO  601 

IF ( 1 1 ♦ EG . 1 ♦ AND ♦ L IH ♦ E  G ♦ 0 ) GO  TO  14 

I F  «  L I H  *  E Q  .  0  )  G 0  TO  6 

CLOCK “RT I M ( T I M2  y T I M 1 ) 

CALL.  T I  ME  ( CLOCK  y  HR  ,  M I N  y  SEC ) 

H 1ST ( 1 , 1 ) ”0 

DO  4  1=2  y  20 

IF < 1 1 *EQ  *  0 ) READ ( 7  y END=6Q1 ) INTEN 
11=0 

A 1  (I >  =  1 F IX  < ( ( INTEN ( N 1 ) -T I M2 ) *32768 ♦  +  ( INTEN ( N2 ) ~T I M 1 ) ) /3600 . ) 

CALL  EXTRAC ( HIST  y I y ELEV  » INTEN  y  N3 y IP  y LIH  y 1 1 ) 

I F  <  L I H ♦ GT ♦ M I H ) M I H=L IH 
IF (II* EG ♦ -1 ) GO  TO  500 

4  CONTINUE 

500  IE'  (  Mil-1  ♦  EG  ♦  0  )  GO  TO  601 

C  A  L  I...  B  T  D  (  H  R  y  A  2  ( 1  >  y  2  y  N  D »  '  0  '  > 

CALL  BTD  <  M 1 N  y  A2  <  3  )  t  2  »  ND  v  -’O') 

CALL  BTD ( SEC  y  A2 ( 5 ) y  2  y ND y  '  0  '  > 

W R I T E < 6  y  2 0 2) PAD , A 2 t R R N G E ( I P ) v P A Z M  < I P ) 

2 0 2  F 0 R M A T  (  '  1  •'  y  6 X  y  '  T I M E - H E I G H T - R E F I... E C T I V I T Y  FOR  P AD  '  y  4A2  y  / 6 X  y 

*'  START  TIME="»2A1»  '  i'ySAlv  '  »2A1»  '  DISTANCE® '  y 1.3  y  '  KM  *  AZIMUTH 

*= ' y  F4 ♦ 0  y '  DEG '  /  y  ' DEG ' y 63X y " KM  * " ) 

DO  5  J~ 1 y 1 5 

LJ R 1 T E  <  6  y  2 0 3  )  E L EG  (  J )  y  ( H  1ST  (  J  y  K  )  y  K  =  .1  y  2 0  )  y  H E I G H T  ( J  > 

203  FORMAT  C  '  SF4,ly/?,y  2013  y  '  ?  "  f  F4  ♦  1  ) 

5  CONTINUE 

W R I TE ( 6  y  204 ) < A 1 ( J ) y J~ 1 y 20 ) 

204  FORMAT  (  '  s  4  X  y  +  +  >  +  +  +  +  #  +  +  *  +  +  +  +  +  y 

y  /  '  TIME  '  y  2013  y  "  MIN''  ////  ) 

M I  H=0 

DO  8  J= 1 y  20 
DO  9  K  =  .1  y  1 5 
HI ST ( K  y  J ) =0 
9  CONTINUE 

8  CONTINUE 

DO  15  J= 1 y  20 
A 1 ( J ) =0 

15  CONTINUE 

IF  (II*  EG  ♦  -•!  )  GO  TO  601 
IF ( LIH ♦ EG ♦ 0  >  GO  TO  6 
I F  <  1 1  ♦  E  G .  0  >  R  E  A  D  (  7  y  E  N  D = 6  0 1 )  I N  T  E  N 
11  =  0 

A 1(1)  =  I F I X ( ( (I NTEN ( N 1 ) -T I  M2 ) *32768+ ( I NTEN ( N2 ) -TI M 1 ) ) 73600 . ) 

GO  TO  7 

600  WRITE ( 6  y  602 ) 

602  F 0 R M  A T  <  ''  E N D - 0 F - F I L E  R E AD  W H I L E  R E A D I N G  P R E A M E L E ? ? ?  '  ) 

REWIND  7 
GO  TO  605 

601  CONTINUE 
REWIND  7 
GO  TO  60S 
END 
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C  THIS  SUBROUTINE  EXTRACTS  DBZ'S  FROM  THE  OUTPUT  TAPE 


SUBROUTINE  EXTRAC  (  HIST  1 1 ,  EL.EV  r  INTEN  »  N3  *  IP »  L IH  r  1 1  > 
REAL  ELEV< 15) » INTEN (53) 

INTEGER  HIST (IS >20) 

ELEU ( 15 ) “INTEN ( N3 ) 

IH=IFIX< INTEN ( IP) ) 

HIST  ( 15  » I )  :-IH 

LIH-IH 

DO  5  J“2  ?  15 

READ ( 7  9  END  -1 ) INTEN 

IF  ( I NTEN  <  N3  )  ♦  EQ  *0*0)  GO  TO  .1 1 

ELEU ( 16~ J ) =INTEN ( N3 ) 

IH-IFIX ( INTEN ( IP ) ) 

HIST ( ( 16- J  > t I ) - 1 H 
IF  ( IH  *  GT  ♦  L. I H  )  L I H" I H 

5  CONTINUE 

READ ( 7  ?  END” 1 ) INTEN 

IF < INTEN ( N3 ) *  EQ  *  22 . 5 ) RETURN 

HIST ( 1 y I ) “0 

1 1  1 1  =  1 

1 3  I F  (  I N TEN  (  N3  )  *  E Q  *  0  *  0  )  R E T U N 

READ ( 7 i END=1 ) INTEN 
GO  TO  13 

1  CONTINUE 

II*-1 
RETURN 
END 


j 


C  A  SUBROUTINE  TO  CONVERT  REAL  TIME  INTO  HOURS  AND  MINUTES 


SUBROUTINE  T I ME ( CLOCK  *  HR » MIN »  SEC) 
INTEGER  HR » SEC 
HR=INT < CLQCK/216000 , ) 

MIN=INT (CL0CK/3600. >-HR*60 
SEC=INT( CLOCK/60 > - < HR*3600f MIN*60 > 
RETURN 
END 


1 


PADEN 


REAL  PAZM ( 50 ) » INTEN(53) *ELEV< 15) »A<2> t HEIGHT ( 15) »ENSUM<4> 


REAL  IHTIM» ILTIM»TIM1 »TIM2» CLOCK 

INTEGER  RRNGE ( 50  > , HIST (15, 20) >ID<4»50) r WRONG (2) ? A 1<20> 
*>PAD<4) *  HR 
DATA  WRONG/ ' ST '  y ' OP ' / 

LOGICAL  EQLJC 
LOGICAL*!  FMLC10) 

LOGICAL.*!  A2  <  6 ) 


RT I M(IHTIH» I LT I M ) = I HT I M*3276B+ ILTIM 


199 


1 1 1 2 

1109 

407 


LQGICAL*1  CHAR ( 28 ) / ' $ ' »/C/»/0/>,N/»/  '  r t ' T ' r ' A' t ' P' t 

*  y  '  P  '  ?  '  0  '  »  '  S '  y  '  N  '  r  '  -  '  / 

WRITE (6? 199) 

FORMAT (  '  ENERGY  INTEGRALS  LISTED  IN  DECREASING  ORDER  OF 

*  D  B  Z  M  5  y  4  0  y  3  5  y  3  0  '  /  /  /  '  PLS  INPUT  NAME  OF  OUTPUT  FILE"  »/> 
CALL  FTNCMD ( ' ASSIGN  4-NEWD. I  RECTORY '  v  26 ) 

CALL  FREAD ( 5  y  ' STR I NG X  '  t CHAR ( 1 8 ) y 1 1 ) 

REWIND  4 


READ < 4 y 407  y END  - 1 1 1 1 ) F IL 
FORMAT (10A1) 


IF ( LCOMC ( 10  y  FIL  y  CHAR ( 18 ) ) »  EQ . 0 ) GO  TO  1110 
GO  TO  1109 


'* 


THRESHOLD 


1 1 1 1  WR I TE  <  6  f 40 1 ) ( CHAR  < I ) » I» 1 8  y 28 ) 

40 1  FORMAT ( '  ' 9 ' F I LE  ' t 1 1 A 1 r ' DOES  NOT  EXIST  *  * PLS  RE  -ENTER  F I LENAME ' ? / ) 
GO  TO  1112 


1110  CALL  FTNCMD ( CHAR , 28) 

CALL  FTNCMD ( ' $COPY  *T APE*  -Y'ylS) 

CALL  FTNCMD ( ' ASSIGN  7=-Y  '  v 1 1 ) 

CALL  FTNCMD ( ^RELEASE  *TAPE*/ v 15) 

605  READ ( 7 , END=600 ) NUMBER » ID , PAZM  ? RRNGE 
C  READ  ID  OF  DESIRED  PAD ♦ 

WRITE  <  6  9 200 ) 

2 0 0  F 0 R MAT  (  y  I N F' U T  D E S I R E D  P A D  I D E N T I F I C A T 1 0 N  0 R  "  S T 0 P  “  '  /  > 

3  READ ( 6  y 100 ) PAD 

100  FORMAT ( 4 A2) 

I F  (  PAD  ( 1 )  ♦  EC-1  ♦  WRONG  ( 1 )  ♦  AMD  *  PAD  (  2  )  *  EQ  .  WRONG  (  2  >  )  STOP 
604  DO  1  IP=1 ?  NUMBER 

DO  80  K=1 9 4 

IF ( PAD ( K ) ♦ NE ♦ ID ( K v IP ) ) GO  TO  1 
80  CONTINUE 

GO  TO  2 

1  CONTINUE 

WRITE ( 6  9 201 ) PAD 

201  FORMAT  PAD  /?4A2y'  CANNOT  BE  FOUND.  PLS  RE-ENTER  IDENTV) 


GO  TO  3 

2  Nl  -NLJMBER  f  1 

N2=NUMBER+2 
N3=NUMBER+3 
DO  10  1=1 » 15 

HE IGHT (16-1)  =  ( RRNGE ( I P ) ) * 1 ♦ 5*  < 1-1 ) *3 ♦ 1 4 15927 /ISO. 

10  CONTINUE 

C  SET  UPPER  AND  LOWER  ELEVATION  BOUNDS  FOR  DBZ  CALCULAlION. 
MAX-0 

DO  20  1=1 » 15 

IF ( HEIGHT ( I) .GE*4*0)G0  TO  20 
MAX  =  I 


GO  TO  21 

20  CONTINUE 

21  MINI-0 

DO  22  I - M A X  > 15 
I F ( HE I GHT < I ) . GE . 1 . 5 ) GO  T 0  22 
MIN I =1-1 
GO  TO  23 


j 


. 
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22  CONTINUE 

23  I NT =MINI~MAXH 
DO  11  J-1,20 

DO  12  K'=MAX  ?  MINI 
HIST  ( K  >  J )  •"<> 

12  CONTINUE 

11  CONTINUE 

6  READ ( 7  y E N D = 6 0 1 ) I N T E N 

14  T I  M2 --IN  TEN  (  N1 ) 

TIM1=INTENCN2) 

M I H  =  0 
11=0 

7  C  A  L  L  E  X  T  R  A  C  (  H I S  T  1 1  »  E  L.  E  V  y 1 N  T  E  N y  N  3 y  I P  ,  L.  I H  y  1 1 ) 

I F  (  L I H  ♦  G  T  *  M I H  )  M I H =L  I H 

IF (II* EQ . -1 ) GO  TO  601 
I F  <  1 1  *  E  G  ♦  1  ♦  A  N  D  *  L I H  *  E  Q  .  0  )  G  0  T  0  1 4 
IFr  ( LIH  ♦  EQ  *  0 ) GO  TO  6 
C  L  0  C  K = R  T I M  <  T I M  2 ,  T I M 1 ) 

C  A  L  L  T I M  E  <  t:  L.  0  C  K  t  H  R  r  M I N  t  S  E  C ) 
l-l  I  S  T  C  1  y  1  >  =0 
DO  4  1=2 y 20 

I F  ( 1 1  *  E  0  ♦  0  )  R  E  A  D  (  7 y E  N  D = 6  0 1  )  I N  T  E  N 
11  =  0 

C A L L  EXT R A C  ( H I S T  y  I  y  E L E U  ,  I N T EN  ?  N3  y  I P  y  L.  1 1-1 » 1 1 ) 

IF  CL  I  IT  ♦  G  T  ♦  M I H  )  M I H = L.  I H 
IF ( 1 1 ♦ EQ  *  -1  )  GO  TO  1300 
4  CONTINUE 

500  IF  ( MIIT  ♦  EIQ  »  0  )  GO  TO  601 

C  A  L  I...  B  T  D  C  H  R  t  A  2  C 1 )  y  2  y  N  D  r  '  0  '  ) 

CALL  BTD  ( M I N  t  A2  ( 3 )  »  2  r  ND  ,  ■'  0  '  ) 

C  A  L  L  B  T  D  -I  S  E  C  t  A  2  C  5 )  ,  2  v  N  D  *  '  0 '  ) 
l/JR  I TE  ( 6  ?  202  )  PAD  y  A2 

2 0 2  F 0 R M A T  C  '  1  '  y  '  H 0 t.J R L  Y  R A D A R - E N E R G Y  I N T E G R A L S  F 0 R  P A D  '  y  4A2  y  / '  S T A R T  T 
#  I  ME=  '  v  2  A  ly'tS  2  A 1 , 'l'  ?  2  A 1  /  ) 

THRESH=45 ♦ 0 
DO  1000  1=1/4 
ENSUMC I )=0» 

DO  1001  ,J=MAX  •/  MINI 
DO  1001  K  = 1 / 20 
FIST=FLOATCHISTC  Jy|<)  ) 

I F  C  Fr  1 3  T  ♦  G  E  <■  T  H  R  E  S  l-l )  E  N  S  U  M  C I )  =  E  N  S  U  M  <  I )  14  *  76*10**  (  F I  ST/12*  -  7  *  )  Y  (  2 
#11  * /FLOAT < I NT)  ) 

1001  CONTINUE 

W R I T  E C  6  y 2 0 3 ) E N  S U MCI) 

203  FORMAT C '  ' y  F 1 0*3) 

1 0 0 0  T  H  R E S  H  =  T  H  R E  S H -5*0 

M I  H=C) 

DO  8  J  =  1 y  2 0 
DO  9  K=MAX  f MINI 
H 1ST  C  K  y  J ) =0 
9  CONTINUE 

8  CONTINUE 
IFCII*EG.--1)GQ  TO  601 
IF  C  LIH ♦ EQ  *  0 ) GO  TO  6 

I F  C  1 1  *  E  Q  *  0  )  R  E  A  D  C  7  r  E  N  D = 6  0 1 )  1 N  T  E  N 
11=0 
GO  TO  7 

600  WRITEC 6/602) 

6 0 2  F 0 R M A T  C  '  E N D 0 F - F 1 1... E  R E AD  LIH  I L E  R E A D I N G  P R E  A M B L E???  '  > 

REWIND  7 
GO  TO  605 

601  CONTINUE 
REWIND  7 
GO  f 0  605 
END 
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C  PROGRAM  TO  GET  DATA  ON  TAPE  AND  CONVERT  IT  TO  DBZ ' S , 

INTEGERS  IHTIMy ILTIM y M < 3072 ) y RBINS < 130 , 20 ) 

REAL  INCEPT 

REAL  RAZM ( 100 ) y MILES (100) 

LOGICAL *1  A ( 6 ) 

LOGICAL  UP r  DOWN  y  SKIP/ ♦ FALSE . / 

INTEGER  MIN  y  SEC  y HR y CNTNUM 
INTEGER  TIMSTT ( 3 ) 

RT I M ( IHT I M  y ILTIM )  =  I HT I M*32768+ 1 LTIM 
BINTIM(HRyMINySEC)=HR*216000*  fMIN*3600.+SEC*60. 

COMMON  RAZMy MILES 

COMMON  SLOPE y INCEPT y RADAR y I PTR 

WRITE ( 6  y  300 ) 

300  FORMAT ( ' 1 INPUT  BEGINNING  TIME ♦ ' ) 

READ ( 5  y 101 ) TIMSTT 

101  FORMAT ( 3 ( 12 y IX )  ) 

WRITE ( 6  y  301 ) 

301  FORMAT  <.//*'  INPUT  SLOPE  y  INTERCEPT  y  RADAR  y  AND  AZIMUTH  (  DEGREES ).'  ) 
CALL  FREAD ( 5 y ' 3R y  I  J  '  y SLOPE » I NCEPT y RADAR y N ) 

N1«N 

N" ( ( N  #  1 ♦  0  )  - 1 4  0  *  0  )  *  6  3 1  .  /  3  6  0  ♦ 

IF  ( N  ♦ LT  .  0 ) N=N+631 * 

IAZ  =  N 

STAR T=B IN TIM( TIMSTT ( 1 ) y TIMSTT (2) y  TIMSTT ( 3 ) ) 

W  R I T  E  (  6  ?  3  0  4  ) 

304  FORMAT (  '  '*//'  INPUT  BEGINNING  ELEVATION  AND  ENDING  ELEVATION .  ' y / ) 

CALL  FREAD (5  y ' 2R  J y y ELEV1 , ELEV2 ) 

TELEV=ELEV1 
WR I TE ( 6  y  305 ) 

305  FORMAT ( '  y  y / / '  INPUT  S T A R T I N G  R A N G E  A N D  E N D I N G  R A N G E ♦  "  » / ) 

CALL  F R E AD  <  5  y ' 21  *  y  y I R 1 1 1 R 2 ) 

WF(ITE  <  6  y  306  ) 

306  FORMAT ('  ',//>'  START  OF  DATA? " Y ( 0 )  OR  N ( 1 ) " ' y / ) 

CALL  FREAD ( 6  y ' IX '  *  III ) 

IF ( 1 1 1 . EG . 1 ) SK I P= *  TRUE ♦ 

I RANGE" I R2~ I R 1 

C  CALCULATE  RAZM  FOR  PRINT-OUT  PURPOSES ♦ 

DO  1  1=1 y 100 

1  RAZM(I)-( (N+CI-l) )*(360*/631. ) )+140 
C  CALCULATE  MILES  FROM  KILOMETERS ♦ 

DO  2  1=1 y I RANGE 

2  MILES (I)=<IR1+(I-1 ) )/l *609 
N  =  I A  Z 

C  BY  NOW  WE  HAVE  READ  EVERYTHING  IN*  LET'S  MOUNT  THE  TAPE  AND  USE  THE 
C  WONDERFUL  SUBROUTINE  SKPREC* 

WRITE<  6  y 100  > 

100  FORMAT ( y  '  t  ' PLS  ENTER  MOUNT*' y/) 

222  CALL  FREAD  <  5  y  7 1  * ' y 1 1 ) 

C  DO  A  BIT  OF  CHECKING* 

I F  < 1 1  * GT  *  0  *  AND  *  1 1  *  LE  *13) GO  TO  224 
WRITE ( 6  y  403 ) 

403  FORMAT ('  INVALID  INPUT*. > PLS  RE-ENTER  MOUNT*'/) 

GO  TO  222 

224  CALL  TF’MNT  (II) 

C  TAPE  MOUNTED. .* READY  FOR  THE  HARD  PART. 

CALL  SKPREC ( 1 1 y START y  M  y  CNTNUM ) 

CALL  DATA ( M  y  CNTNUM  y S40 y %27 ) 

DOWN= . TRUE » 

UP=. FALSE. 

C  SPACE  FORWARD  TO  START  TIME  AND  GET  DOWN  PULbt * 

11  CALL  HED I NG ( M  y  CNTNUM  y DOWN  y UP y X 1 1  *  82 / y  827 ) 


CALL  BYTREV (My CNTNUM+1) 

CALL  BYTREV ( M , CNTNUM+2 ) 

CLQCK-RTIM ( M (  CNTNUM+2 ) y M ( CNTNUM+1 ) ) 

IF ( CLOCK . GE ♦ START ) GO  TO  14 

CNTNUM=CNTNUM+2 

GO  TO  11 

14  IF (SKIP) GO  TO  IS 

CALL  SKPRAY  <  M »  CNTNLJM  y  200  r  827  r  827  ) 

C  SET  ELEVATION  COUNTER  TO  ZERO  AND  SKIP  TO  NEXT  HEADING  MARKER ♦ 
1 6  CALL  MED ING  (  M  y  CNTNLJM  t  .  FALSE  .  y  .  TRUE  .  ,  8 1 6  t  827  t  827  ) 
CNTNUM=CNTNlJM  +  2 

1 5  I F  (  EL.EV  1  ♦  EQ  .0.0)  GO  TO  30 
DOWN- . FALSE . 

UP= . TRUE . 

C  GET  PROPER  ELEVATION. 

ELEV=0.0 

2 6  C A L L  S K P R A  Y  (  M  y  CN T N I.J M  >■  2 00  y  8 2 7*827) 

C  A  L  L  H  E  D I N  G  (  M  *  C  N  T  N  LJ  M  *  D  0  U  N  *  l  J  P  *  8  2  7  ) 

ELEV-ELEV+1 .5 

I F  ( E  L.  E  V  .  L.  T  .  E  L.  E  V 1  )  G  0  T  0  2  6 
3 0  C A L L  3 K P R A  Y  <  M  *  C N T N LJM  y  N  y  827  *  8 2 7  ) 

N--IAZ 

C  A  L  L  B  Y  T  R  E  V  ( M  » C  N  T  N  U  M  +  1  ) 

C  A  L  L  B  Y  T  R  E  V  (  M  *  C  N  T  N  LJ  M + 2  ) 

C L 0 C K - R T I M  ( M  (  C N T N LJ M  +  2  )  *  M  (  C N T  NIJM+1  )  ) 

C A L L  T I M E ( C L 0 C K  y  HR  y  MIN  y S E C ) 

C  A  L  L.  B  T  D  (  H  R  y  A  (  1  )  y  2  r  N  D  t  7  0  7  ) 

C  A  L.  L  B  T  D  (  M I N  t  A  (  3  )  r  2  y  N  D  y  7  0  7  ) 

C  A  L.  L  B  T  D  (  S  E  C  y  A  <  5  )  r  2  r  N  D  y  7  0 7  ) 

CNTNUM=CNTNUM+2 

C  TAPE  IS  NOW  POSITIONED  AT  PROPER  TIME , ELEVATION y AND  AZIMUTH. 

C  LET'S  PRINT  SOME  STATISTICS. 

W  R I T  E  <  6  y  4  0 1  )  A  t  N 1  y  I R  A  N  G  E  y  E  L  E  V 1 
4  0 1  F  0  R  MAIL  /  /  y  4  X  y  7  S  T  A  R  T I N  G  T I M  E  =  >  7  y  2  A 1  r  7  J  7  r  2  A 1  y  7  t  7  r  2  A 1  y 

*  /  9  4  X  y  7  S  T  A  R  T I N  6  A  Z I M  LJ  T  H = >  7  r  1 4  y  7  D  E  6  R  E  E  S  7  y 

*  /  y  4  X  y  7  A  Z I M  U  T  H  C  0  V  E  R  E  D  «  >  6  0  D  E  G  R  E  E  3  ( 1 0  0  R  A  Y  S  )  7  y 

*  /  y  4  X  y  7  R  A  N  G  E  C  0  V  E  R  E  D  •=  >  7  t  1 3  y  7  K  M  7 

%/  y  4  X  y  7  S  T  A  R  T I N  G  E  L.  EVA  T 1 0  N = >  7  y  F  5 . 2  >•  7  D  E  G  R  E  E  S 7  y  /  / ) 

2  3  W  E I T  E  (  6  y  4  0  0  )  <  E  A  Z  M  (  I P  T  R  )  y  I P  7'  R  =  1  t  2  0  ) 

4  0  0  F  0  R  M  A  T  <  /  y  1 X  r  7  A  Z I M  LJ  T  H  7  y  2  0  (  F  5 . 1  *  1 X  )  ) 

WR I TE  <  6  y  404 ) ( I y I - 1 y 20 ) 

4  0  4  F  0  R  M  A  T  ( 3  X  r  7  R  A  Y  *  7  1 3  X  y  ( 2  0  ( 1 3  y  3  X )  )  ) 

C  LET'S  DO  SOME  WORK  WITH  THOSE  BINS. 

CALL,  D  B  Z  R  A  Y  <  M  y  C  N  T  N  U  M  r  R  B I N  S  y  I R  A  N  G  E  y  I R 1  y  I R  2  »  8  2  7  y  8  2  7*  84  0  ) 

I F  (  E  L  E  V 1  .  E  Q  .  E  L  E  V  2 )  G  0  T  0  4  0 

C A L. L  H ED  I NG  (  M  y  C N T N U M  y  .  F A L S E  .  y  .  T R U E  .  y  8 2 7 ) 

ELEV1=ELEV1+1.5 

C  SINCE  WE  ARE  DEALING  WITH  A  NP.ySKIP  TO  PROPER  AZIMUTH. 

I F  ( E  L  E  V 1  .  G  T  .  E  L  E  V  2  )  G  0  T  0  4  0 
C  A  L.  L.  S  K  P  R  A  Y  <  M  y  C  N  T  N  LJ  M  r  N  t  8  2  7  »  8  2  7  ) 

N-IAZ 

CALL  BYTREV ( M  y CNTNUM  +  1 ) 

CALL  BYTREV <  M  y CNTNUM+2 ) 

C  L,  0  C  K  -  R  T I M  <  M  (  C  N  T  N  U  M  +  2  )  y  M  (  C  N  T  N  LJ  M  + 1  )  ) 

CNTNUM-CNTNUM+2 
C  A  L  L.  T I M  E  (  C  L.  0  C  K  y  H  R  y  M I N  y  S  E  C  ) 

C  A  L  L  B  T  D  (  H  R  y  A  <  1  )  t  2  y  N  D  y  7  0  7  ) 

CALL  BTD ( M I N  y  A ( 3 ) y  2  y  ND  y  7  0  7  ) 

CALL  BTD ( SEC y A ( 5 ) y 2  y ND  y  7  0  7  ) 

WRITE  ( 6  y  407 )  A  y  EL.EV1 

407  FORMAT < / y 7 T I ME*> 7  y  2A1 y  7 i 7  y  2A1 y  7 ♦ 7  y  2A1 y 
%/ 7  EL..EV  AT  1 0N~>  7  y  F5 . 2  y  /  ) 

GO  TO  23 

27  WRITE ( 6  y  408 ) 

408  FORMAT ( 7  TROUBLES . . . TROUBLES . . 7 ) 

STOP 
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40  CALL  HEDING  (  M  r  CNTNUM  »  ♦  TRUE .  y  *  FALSE  ♦  r  &40  *  ) 

• CNTNUM=CNTNUM+2 
ELEVlaTELEV 
GO  TO  15 
END 


co  m 
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PERMS 


DIMENSION  ERR ( 3 ) rK0UNT(3> »RER<3) »SIGN<3) ,0BS<60) ,  OBSR (60) ,EST(3rA0 
*> t PER (3) 

COMMON  /RN/IPERM ( 60 ) 

READ  <  3  » 1 )  (  (EST(IrJ)  r 1  =  1 , 3 )  r J=1 , 60 ) 

1  FORMAT  <  3F8 . 3 ) 

READ <  4  »  2 ) <OBS< I ) r 1=1 rAO) 

2  FORMAT ( 10F6  » 1 ) 

CALL  TIME(OfO) 

DO  3  K  -l >  3 

PER ( K ) =0 » 

ERR ( K ) =0 ♦ 

RER  <  K ) =0 ♦ 

3  KOUNT (  K  )  =  0 

C  CALCULATION  OF  ACTUAL  MEAN  ERRORS  OF  ESTIMATION  METHODS 
DO  4  J  =  1  r  A  0 
DO  4  1=1 » 3 

4  ERR ( I ) =ERR ( I ) TABS ( EST ( I v  J ) -QBS ( J ) )/60. 

ISEED=16807 

C  CALCULATION  OF  MEAN  RANDOMIZED  ERRORS  (REPEATED  NUMEROUS  TIMES) 

DO  5  1=1 f 10000 
CALL  PERM (I SEED) 

DO  6  J  =  1 y  60 
JJ=IPERM< J) 

6  OBSR ( J ) =QBS ( J J ) 

DO  7  J  =  1 rAO 

DO  7  K  =  1 r  3 

7  RER ( K ) =RER ( K ) TABS ( OBSR ( J ) -EST ( K  » J ) ) /AO ♦ 

C  COMPARISON  OF  MEAN  ACTUAL  ERROR  WITH  MEAN  RANDOMIZED  ERROR ♦ THIS  IS 
C  TALLIED  IF  THE  MEAN  RANDOMIZED  ERROR  IS  LESS  THAN  THE  MEAN  ACTUAL 
C  ERROR ♦ 

DO  8  L=lv3 

PER  ( L. )  =  PER  (  L  )  +  RER  (  L  )  /  1 0000 
IF ( RER ( L ) . LT . ERR ( L ) ) KOUNT ( L ) =KOUNT ( L ) + 1 
RER  <  L ) =0 ♦ 

CONTINUE 
DO  9  1  =  1  *  3 

9  SIGN  (I  )=KOUNT<  D/100* 

CALL  TIME<15»1) 

WRITE (6, 10) ( SIGN ( 1 ) > 1  =  1 »  3 ) 

10  FORMAT ( '  PERCENTAGE  OF  CASES  IN  WHICH  MEAN  ERROR  OF  ESTIMATE  BY  YA 
*RIOUS  METHODS  WAS  LESS  THAN  MEAN  ERROR  OF  RANDOMLY  PERMUTED  OBSERU 
NATIONS'//'  METHOD  1  METHOD  2  METHOD  3 ' //3 ( F5 ♦ 2  >  AX ) ) 

C  LIST  POPULATION  MEAN  RANDOMIZED  ERRORS  FOR  THE  RESPECTIVE  METHODS 
WRITE  (  Ay  11  )  (  PER  (  L. )  rL  =  l>3) 

FORMAT ( 'O' y  3 F 1 0 . 3 ) 

STOP 
END 
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SUBROUTINE  PERM < I SEED ) 
COMMON  /RN/IU ( 60 ) 

DO  20  1=1,60 
IV< I )=I 

CALL  RNDF'  ( I  SEED  ) 

RETURN 

END 
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SUBROUTINE  RNDP(ISEED) 

DIMENSION  D ( 60 ) 

COMMON  /RN/IU<60) 

N  =  60 

C  GGU3  IS  A  *IMSL  ROUTINE  THAT  GENERATES  RANDOM  NUMBERS. 
CALL  GGU3(ISEED»N»D) 

DO  20  1=2*60 
j~N-I+l 
K=J*D< I )+l 
UI  =  IO  ( K ) 

IV<K)=IV( J) 

20  IU(J)=U 

RETURN 
END 
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APPENDIX  C 


Hailpad  triplets  are  described  on  p.  of  the  text. 

In  the  following  listing,  "hailpad  I.D.’s"  designated  by  W,X,Y  or  Z 
refer  to  locations  known  to  have. recei ved  no  hail,  but  are  not  offi- 


c  i  a 

hailpad  1 ocat i ons . 

These 

are: 

W3 : 

SW6, 43, 1,5 

XI : 

SW30, 30,27,4 

Z3 

SW33, 32,2,5 

W4 : 

NE32, 45, 27, 4 

X2: 

NE 1 2, 30, 28, 4 

Z4 

SW 1 0, 33, 3, 5 

W5: 

SW3, 45,  2,  5 

X3: 

NE29, 3 1 , 28, 4 

Z6 

NW27, 34,27,4 

W6: 

SW2, 46, 26, 4 

X4: 

SE6, 31 , 1,5 

Z7 

SW30, 35, 28, 4 

W7 : 

SE 1 7, 46, 27, 4 

X6: 

SE 1 6, 4 1 , 25, 4 

Z8 

NE2 1,34,5,5 

W8 : 

SW21 ,40,5, 5 

Y4: 

SE3, 33,2,5 

Z9 

NW7, 37, 1 , 5 

Kinetic  energy  densities  are  calculated  from  Equations  2.2  and  2.3. 


Date 

Triplet 

Ha  i 

1  pad 

.D.  Ts 

Ki  net i c  Energy 
Dens i t ies ( Jm-2 ) 

D  i  stances 
(mi Ies) 

(1974) 

1  .D. 

H  I 

HC 

h2 

El 

Ec 

e2 

d  i 

d2 

July  30 

3001 

A 1  8 

B21 

C27 

23.9 

14.9 

0.5 

7.0 

14.0 

3002 

B21 

C27 

D31 

14.9 

0.5 

1  .2 

14.0 

8.5 

3003 

C27 

D31 

D36 

0.5 

1.2 

45.7 

8.5 

8.0 

3004 

E27 

D31 

D33 

1.6 

1.2 

28.6 

6.5 

5.5 

3005 

W4 

D33 

D34 

0 

28.6 

6.0 

1.0 

3.5 

3006 

B33 

W7 

D33 

43.4 

0 

28.6 

9.0 

2.0 

3007 

B21 

B33 

B36 

14.9 

43.4 

0.9 

24.0 

5.0 

3008 

B36 

B39 

C39 

0.9 

1.2 

1  .  1 

4.5 

5.5 

3009 

B39 

C36 

W6 

1.2 

10.8 

0 

4.0 

3.0 

3010 

B39 

C36 

D34 

1.2 

10.8 

6.0 

4.0 

8.0 

301  1 

D33 

D34 

D36 

28.6 

6.0 

45.7 

3.5 

3.0 

3012 

W5 

C27 

B33 

0 

0.5 

43.4 

7.5 

16.0 

3013 

C27 

W7 

C36 

0.5 

0 

10.8 

10.0 

9.0 

3014 

W7 

C36 

C39 

0 

10.8 

1.1 

9.0 

6.5 

3015 

B39 

C39 

B46 

1.2 

1 . 1 

1  .9 

5.5 

9.0 

3016 

C36 

C39 

B46 

10.8 

1  .  1 

1.9 

6.5 

9.0 

3017 

D34 

W6 

C39 

6 . 0 

0 

1  .  1 

6.0 

8.0 

3018 

D34 

C36 

C39 

6.0 

10.8 

1  .  1 

8.5 

6.5 

3019 

C39 

B46 

B51 

1.1 

1.9 

9.6 

9.0 

14.5 

3020 

W2 

E27 

D36 

0 

1.6 

45.7 

CO 

3.5 

it i nuec 

13.0 
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July  30 

3021 

G37 

G36 

F42 

0.8 

19.2 

12.5 

1.5 

8. 

0 

(Con+ ’ d ) 

3022 

F42 

F46 

G42 

12.5 

52.3 

27.5 

7.0 

8. 

0 

3023 

G42 

G43 

H48 

27.5 

27.7 

48.7 

4.5 

7. 

5 

3024 

D48 

D49 

G42 

10.0 

10.0 

27.5 

3.0 

14. 

5 

3025 

D49 

D48 

B51 

10.0 

10.0 

9.6 

3.0 

10. 

5 

3026 

D49 

D51 

E57 

10.0 

8.9 

86.7 

4.0 

1  1 

0 

3027 

D48 

D51 

E57 

10.0 

8.9 

86.7 

5.5 

1  1 

0 

3028 

B51 

C48 

E57 

9.6 

3.  1 

86.7 

6.0 

16. 

5 

3029 

B46 

D48 

D51 

1.9 

10.0 

8.9 

14.5 

4. 

0 

3030 

B46 

D48 

E57 

1.9 

10.0 

86.7 

14.5 

15. 

0 

3031 

A 1  8 

C27 

D31 

23.9 

0.5 

1  .2 

21.0 

8. 

5 

3032 

A 1 9 

B21 

C27 

76.6 

14.9 

0.5 

10.0 

14. 

0 

3033 

E27 

D31 

W4 

1  .6 

1.2 

0 

6.5 

5. 

5 

3034 

B36 

C36 

W6 

0.9 

10.8 

0 

7.5 

3. 

0 

3035 

1  12 

1  14 

W8 

126.3 

64.  1 

0 

6.0 

2. 

5 

3036 

D33 

C36 

C39 

28.6 

10.8 

1 . 1 

9.5 

6 

5 

3037 

W4 

C36 

C39 

0 

10.8 

1.1 

10.0 

6 

5 

3038 

W3 

E27 

D31 

0 

1  .6 

1.2 

10.0 

6. 

5 

August  7 

701 

S9 

SI  1 

X2 

23.3 

4.2 

0 

3.  5 

7. 

5 

702 

S  1  2 

XI 

R20 

1.8 

0 

0.1 

6.0 

5. 

5 

703 

R4 

X4 

SI  1 

5.4 

0 

4.2 

5.5 

6. 

5 

704 

S 1  2 

X3 

Q 1 6 

1.8 

0 

2.9 

6.5 

6 

5 

705 

X3 

Q 1 6 

P 1 7 

0 

2.9 

58.1 

6.5 

7. 

0 

706 

R20 

Q  1  8 

P21 

0.  1 

0.5 

19.1 

10.0 

7. 

5 

707 

Q4 

R4 

X4 

8.2 

5.4 

0 

7.0 

4. 

5 

708 

R4 

Q9 

A  1 6 

5.4 

0.3 

2.9 

4.5 

9. 

5 

709 

00 

P 1 2 

Q 1 6 

3.0 

0 

2.9 

6.0 

13. 

5 

710 

00 

P 1 2 

Q4 

3.0 

0 

8.2 

6.0 

15 

5 

71  1 

NO 

Q 1 6 

R20 

2.0 

2.9 

0.1 

22.5 

14 

0 

712 

NO 

P 1  7 

Q 1 8 

2.0 

58.1 

0.5 

19.5 

16 

5 

714 

F 1 8 

G28 

H37 

1.8 

23.9 

8.7 

19.0 

14 

5 

715 

F 1  8 

G27 

H37 

1.8 

4.5 

8.7 

18.0 

15 

5 

716 

H37 

X6 

G43 

8.7 

0 

16.5 

7.5 

9 

0 

717 

G27 

F30 

E33 

4.5 

19.5 

22.4 

7.0 

6 

0 

718 

G28 

F30 

E33 

23.9 

19.5 

22.4 

9.0 

6 

.0 

719 

F30 

F33 

F36 

19.5 

17.3 

8.2 

3.5 

3 

0 

720 

F34 

F36 

D36 

15.1 

8.2 

50.2 

3.5 

7 

0 

721 

F33 

E33 

D33 

17.3 

22.4 

1.6 

5.5 

6 

.0 

722 

D37 

E33 

F36 

166.1 

22.4 

8.2 

5.5 

6 

.0 

723 

D31 

E33 

F36 

25.6 

22.4 

8.2 

5.5 

6 

.0 

724 

D33 

D36 

E36 

1.6 

50.2 

4.5 

5.0 

4 

.5 

725 

F33 

F36 

F39 

17.3 

8.2 

21 . 1 

3.0 

8 

.5 

726 

D36 

E36 

F39 

50.2 

4.5 

21  .  1 

4.5 

7 

.0 

727 

E33 

D36 

C42 

22.4 

50.2 

18.9 

cc 

4.0  14 

:>nt  i  nued 

.5 
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August  7 
(Cont f d ) 


August  15 
August  18 


728 

F34 

E36 

C42 

15.1 

4.5 

18.9 

8.5 

14.0 

729 

F39 

G39 

X6 

21  .  1 

59.3 

0 

7.0 

8.0 

730 

E36 

G39 

G43 

4.5 

59.3 

16.5 

13.0 

12.5 

731 

C42 

D45 

D49 

18.9 

0.4 

8.5 

7.5 

9.5 

733 

F30 

E33 

D36 

19.5 

22.4 

50.2 

6.0 

4.0 

734 

D37 

E33 

E36 

166.1 

22.4 

4.5 

5.5 

6.5 

735 

D31 

E33 

E36 

25.6 

22.4 

4.5 

5.5 

6.5 

736 

D37 

D36 

E36 

166.1 

50.2 

4.5 

7.5 

4.5 

737 

D31 

D36 

E36 

25.6 

50.2 

4.5 

8.0 

4.5 

738 

F30 

E36 

D45 

19.5 

4.5 

0.4 

10.5 

17.5 

739 

F33 

E36 

D45 

17.3 

4.5 

0.4 

7.5 

17.5 

740 

F30 

D36 

C42 

19.5 

50.2 

18.9 

9.5 

14.5 

741 

F33 

E36 

C42 

17.3 

4.5 

18.9 

7.5 

14.0 

742 

D37 

D33 

C42 

166.1 

1  .6 

18.9 

4.5 

16.5 

743 

D31 

D33 

C42 

25.6 

1.6 

18.9 

5.0 

16.5 

744 

D31 

D37 

D36 

25.6 

166.1 

50.2 

0.5 

7.5 

745 

F36 

E36 

C42 

8.2 

4.5 

18.9 

5.0 

14.0 

153 

Q82 

Q80 

Y4 

2.9 

0.4 

0 

10.0 

4.0 

1801 

K21 

L2 

L4 

4.5 

26.8 

63.8 

8.0 

2.5 

1802 

J  6 

K21 

L2 

0.0 

4.5 

26.8 

12.5 

8.0 

1803 

L2 

L4 

LI  8 

26.8 

63.8 

33.0 

2.5 

6.0 

1804 

L0 

L2 

L4 

16.1 

26.8 

63.8 

13.0 

2.5 

1805 

MO 

L 1 8 

K6 

14.5 

33.0 

0.1 

8.5 

5.5 

1806 

K21 

K6 

K8 

4.5 

0.  1 

0 

17.5 

2.5 

1807 

L4 

L 1 8 

M4 

63.8 

33.0 

101  .9 

6.0 

5.0 

1808 

K6 

M4 

N80 

0.  1 

101.9 

41 .6 

8.0 

7.0 

1809 

L 1  8 

Ml 

N2 

33.0 

39.8 

17.5 

6.0 

7.0 

1810 

Ml 

N2 

04 

39.8 

17.5 

22.4 

7.0 

3.5 

1811 

N2 

04 

05 

17.5 

22.4 

33.3 

3.5 

4.0 

1812 

04 

05 

R1 

22.4 

33.3 

41  .9 

4.0 

12.5 

1813 

Q 1 

R 1 

R8 

13.9 

41.9 

13.7 

8.5 

7.0 

1814 

R8 

Q9 

Q 1 6 

13.7 

10.5 

14.1 

8.5 

9.5 

1815 

Q 1 6 

Q 1  8 

Q31 

14.1 

0.4 

0.6 

15.5 

23.0 

1816 

P 1 7 

017 

018 

0.5 

3.2 

2.3 

2.0 

3.0 

1817 

N9 

016 

018 

18.0 

13.1 

2.3 

5.0 

6.5 

1818 

N9 

016 

017 

18.0 

13.  1 

3.2 

5.0 

5.0 

1819 

N9 

016 

P 1  7 

18.0 

13.  1 

0.5 

5.0 

5.5 

1820 

K8 

Ml  3 

07 

0 

127.7 

15.  1 

1  1.0 

1  1.5 

1822 

K8 

M 1 3 

N 

0 

127.7 

7.3 

1  1.0 

7.0 

1823 

L 1  8 

M4 

M 1 3 

33.0 

101  .9 

127.7 

5.5 

6.5 

1824 

M4 

Ml  3 

N9 

101.9 

127.7 

18.0 

6.5 

12.5 

1825 

M4 

M 1  3 

07 

101  .9 

127.7 

15.1 

6.5 

1  1  .5 

1826 

Ml 

N80 

N 

39.8 

41.6 

7.3 

5.5 

4.5 

1827 

N2 

N80 

Ml  3 

17.5 

41  .6 

127.7 

4.0 

6.0 

conti nued 
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August  18 

1828 

L4 

Ml 

N80 

63.8 

39.8 

41.6 

10.0 

5.5 

(Cont ' d ) 

1829 

N2 

N 

07 

17.5 

7.3 

15.1 

6.0 

8.5 

1831 

MO 

N22 

04 

14.5 

56.8 

22.4 

5.0 

9.0 

1832 

N22 

Ml 

M4 

56.8 

39.8 

101  .9 

7.5 

5.0 

1833 

Ml 

M4 

K8 

39.8 

101  .9 

0 

4.5 

10.5 

1834 

L4 

Ml 

N 

63.8 

39.8 

7.3 

10.0 

5.5 

1835 

MO 

N22 

05 

14.5 

56.8 

33.3 

5.0 

12.0 

1837 

05 

Z3 

Q9 

33.3 

0 

10.5 

1  1.0 

6.5 

1838 

Z4 

N 

M 1  3 

0 

7.3 

127.7 

12.5 

7.5 

1839 

07 

016 

Z6 

15.1 

13.  1 

0 

7.5 

6.5 

1840 

07 

N9 

Z7 

15.  1 

18.0 

0 

5.0 

2.0 

1841 

Z8 

N22 

LI  8 

0 

56.8 

33.0 

6.5 

1  1.5 

1842 

M4 

M13 

11 

101  .9 

127.7 

0 

6.5 

12.0 

1843 

N80 

Ml  3 

19 

41.6 

127.7 

0 

6.0 

7.5 

Appendix  D 


CENTER-PAD  ESTIMATES  AND  ERRORS  OF  TRIPLETS  (IN  Jm" 2  ) 


The  following  is  a  list  of  the  center-pad  estimates  and 
the  errors  of  these  estimates  described  on  pp.  72  and  74  of  the  text. 


Triplet 

1  .D. 

Method  1 
Est i mate 

lEc’-Eol 

Method  2 
Est i mate 

IEc’-EcI 

Method  3 
Estimate 

1 ECT-E 

-c  1 

701 

1  1.7 

7.5 

15.9 

11.7 

5.7 

1.5 

702 

1  .0 

1.0 

0.9 

0.9 

5.0 

5.0 

703 

4.8 

4.8 

4.8 

4.8 

5.0 

5.0 

704 

2.4 

2.4 

2.4 

2.4 

5.0 

5.0 

705 

29.1 

26.2 

28.0 

25.1 

5.0 

2. 1 

706 

9.6 

9.1 

1 1  .0 

10.5 

6.7 

6.2 

707 

4.  1 

1.3 

3.2 

2.2 

5.0 

0.4 

708 

4.2 

3.9 

4.6 

4.3 

5.0 

4.7 

709 

3.0 

3.0 

3.0 

3.0 

5.8 

5.8 

710 

5.6 

5.6 

4.6 

4.6 

5.8 

5.8 

711 

1  .  1 

1.9 

0.8 

2.  1 

5.0 

2.1 

712 

1.3 

56.9 

1.2 

56.9 

6.4 

51.7 

714 

5.3 

18.7 

5.7 

18.2 

21  .3 

2.6 

715 

5.3 

0.8 

5.5 

1.0 

65.3 

60.8 

716 

12.6 

12.6 

12.2 

12.3 

5.0 

5.0 

717 

13.5 

6.  1 

14.1 

5.4 

12.5 

7.0 

718 

23.2 

3.7 

23. 1 

3.6 

12.5 

7.0 

719 

13.9 

3.5 

13.4 

3.9 

19.9 

2.6 

720 

32.7 

24.5 

26.8 

18.6 

13.8 

5.6 

721 

9.5 

13.0 

9.8 

12.6 

24.4 

2.0 

722 

87.2 

64.8 

90.6 

68.2 

24.4 

2.0 

723 

16.9 

5.5 

17.3 

5.  1 

24.4 

2.0 

724 

3.  1 

47.2 

3.  1 

47.  1 

36.3 

14. 

725 

19.2 

1  1  .0 

18.3 

10.1 

13.8 

5.6 

726 

35.7 

27.5 

38.8 

30.6 

7.8 

5 . 

727 

20.7 

29.6 

21  .6 

28.6 

5.0 

14. 

728 

17.0 

12.5 

16.5 

12.0 

5.0 

5.5 

729 

10.6 

48.8 

1 1.8 

48.0 

6 . 2 

54.3 

730 

10.5 

48.8 

10.6 

48.7 

24.4 

54.1 

5 

731 

13.7 

13.3 

14.3 

13.9 

24.4 

5.  / 

733 

34.9 

12.5 

37.9 

15.5 

24.4 

cont 

ZA 

nued 
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734 

85.3 

62.9 

735 

15.  1 

7.4 

736 

85.3 

35.1 

737 

15.  1 

35.2 

738 

10.0 

5.5 

739 

8.9 

4.4 

740 

19.2 

31.0 

741 

18. 1 

13.6 

742 

92.5 

90.9 

743 

22.3 

20.7 

744 

108.2 

82.6 

745 

13.6 

9.1 

153 

1.5 

1.1 

1801 

34.2 

7.4 

1802 

13.4 

8.9 

1803 

29.9 

33.9 

1804 

40.0 

13.2 

1805 

7.3 

25.7 

1806 

2.3 

2.2 

1807 

82.9 

49.9 

1808 

20.9 

81 . 1 

1809 

25.3 

14.6 

1810 

31 . 1 

13.6 

181  1 

25.4 

3.0 

1812 

32.2 

1.2 

1813 

13.8 

28.  1 

1814 

13.9 

3.4 

1815 

7.4 

7.0 

1816 

1.4 

1.8 

1817 

10.2 

3.0 

1818 

10.6 

2.5 

1819 

9.3 

3.9 

1820 

7.6 

120.2 

1822 

3.7 

124.  1 

1823 

80.4 

21.6 

1824 

60.0 

67.8 

1825 

58.5 

69.2 

1826 

23.6 

18.1 

1827 

72.6 

31  .0 

1828 

52.7 

12.9 

1829 

16.3 

9.0 

1831 

18.5 

38.4 

1832 

79.4 

39.6 

1833 

19.9 

82.0 

1834 

35.6 

4.3 

.0 

69.6 

24.0 

2.0 

.9 

6.5 

24.4 

2.0 

.  1 

14.9 

36.3 

14.1 

.  1 

38.1 

36.3 

14.1 

.3 

7.8 

7.8 

3.3 

.2 

7.7 

7.8 

3.3 

.3 

30.9 

36.3 

14.1 

.9 

13.4 

7.8 

3.3 

.6 

133.0 

10.0 

8.4 

.0 

22.4 

10.0 

8.4 

.9 

133.3 

166.1 

155.0 

.0 

6.5 

7.8 

3.3 

.8 

0.4 

.7 

22.9 

28.7 

1.9 

.3 

1  1  .8 

0.8 

3.7 

.6 

35.2 

44.9 

18.9 

.  1 

29.3 

28.7 

1  .9 

.8 

27.2 

24.5 

8.5 

.  6 

1.5 

5.5 

5.4 

.  6 

51.6 

24.5 

8.5 

.2 

79.7 

.8 

14.0 

.2 

10.7 

.9 

2.5 

.1 

6.2 

.8 

28.  1 

.9 

3.4 

.5 

2.  1 

0.0 

0.4 

.2 

2.0 

1  .7 

1.5 

.2 

1  .9 

0.0 

13.  1 

.  6 

2.5 

0.0 

13.  1 

.7 

3.4 

0.0 

13.  1 

.4 

120.3 

.5 

123.2 

.4 

25.5 

.2 

54.5 

.6 

57.  1 

.9 

19.7 

47.5 

5.9 

.  6 

20.0 

47.5 

5.9 

.5 

9.7 

.5 

9.2 

.3 

39.5 

.9 

44.  1 

.9 

74.0 

.3 

12.5 

cont  i  nued-- 

92 

15 

65 

12 

12 

12 

19 

17 

134 

24 

158 

1  1 

0 

49 

16 

28 

56 

5 

1 

84 

22 

25 

28 

24 

27 

13 

13 

2 

1 

1 1 

10 

9 

7 

4 

76 

73 

70 

21 

61 

49 

16 

17 

83 

27 

27 
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1835 

23.9 

32.9 

20.0 

36.8 

1837 

21  .9 

21.9 

19.0 

19.0 

0.0 

0.0 

1838 

63.9 

56.6 

79.8 

72.5 

10.3 

3.0 

1839 

7.6 

5.6 

7.0 

6.1 

0.0 

13.  1 

1840 

7.6 

10.5 

4.3 

13.7 

1841 

16.5 

40.3 

1 1  .9 

44.9 

1842 

51.0 

76.8 

66 .  1 

61.6 

1843 

20.8 

106.9 

23.  1 

104.6 

3001 

12.2 

2.7 

16.  1 

1.2 

3002 

8.  1 

7.6 

6.4 

5.9 

3003 

23.  1 

21.9 

23.8 

22.6 

3004 

15.1 

13.9 

16.2 

15.0 

3005 

3.0 

25.6 

1  .3 

27.3 

3006 

36.  1 

36.  1 

31  .4 

31.4 

3007 

7.9 

35.5 

3.3 

40.  1 

3008 

1  .  1 

0.2 

1  .0 

0.2 

3009 

0.6 

10.2 

0.5 

10.3 

3010 

3.6 

7.2 

2.8 

8.0 

3011 

37.2 

31.2 

37.8 

31.8 

3012 

21.7 

21  .2 

13.9 

13.4 

3013 

5.7 

5.7 

5.9 

5.9 

3014 

0.6 

10.2 

0.7 

10.1 

3015 

1.6 

0.5 

1.5 

0.4 

3016 

6.4 

5.3 

7.  1 

6.0 

3017 

3.6 

3.6 

3.9 

3.9 

3018 

3.6 

7.3 

3.2 

7.6 

3019 

5.4 

3.5 

4.4 

2.5 

3020 

22.9 

21.3 

9.7 

8.  1 

3021 

6.7 

12.6 

2.6 

16.6 

3022 

20.0 

32.3 

19.5 

32.8 

3023 

38.6 

10.9 

35.8 

8.  1 

3024 

18.8 

8.8 

13.0 

3.0 

3025 

9.8 

0.2 

9.9 

0.1 

3026 

48.4 

39.5 

30.5 

21.6 

3027 

48.4 

39.5 

35.6 

26.7 

3028 

48.2 

45.  1 

30.2 

27.  1 

3029 

5.4 

4.6 

7.4 

2.6 

3030 

44.3 

34.3 

43.6 

33 . 6 

3031 

12.6 

12.  1 

7.7 

7.2 

3032 

38.6 

23.7 

44.9 

30.0 

3033 

0.8 

0.4 

0.7 

0.5 

3034 

0.5 

10.4 

0.3 

10.5 

3035 

63.2 

1  .0 

37.  1 

27.0 

3036 

14.9 

4.  1 

12.3 

1  . 5 

3037 

0 . 6 

10.3 

0.7 

10.  1 

3038 

0 . 6 

1.0 

0.7 

0.9 
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Appendix  E 


The  table  on  the  following  pages  lists  all  the  analyses 
used  in  determining  the  accuracy  of  registration  of  the  C-band  radar 
output.  The  results  of  these  analyses  are  summarized  in  Secion  6.3 
of  the  text.  The  radar  data  from  the  hail-free  locations  Z1  to  Z5 
and  Z7  to  Z9  were  studied,  and  used  for  registration  analyses.  The 
land  locations  Z3,  Z4,  Z5,  Z7,  Z8  and  Z9  have  already  been  noted  in 
Appendix  C.  The  two  remaining  locations  are: 

Z 1 :  SW26, 33, 1 , 5  Z2:  NW36,33, 1,5 


INDEX  TO  TABLE  HEADINGS 


rm 

am 

rPLOT 
aPLOT 
rPP  I 

aPp| 

PLOT 


X 

X. 

X 

X 

X. 

X 


R 

PLOT 

T 

PPI 

R 

PPi 

T 

PLOT 

TOTAL 

PPI 

TOTAL 


map  range  in  kilometers 

map  azimuth  in  degrees 

range  given  by  program  PLOT  kilometers) 

azimuth  given  by  PLOT  (degrees) 

range  given  by  PPI 

azimuth  given  by  PPI 

range  error  in  kilometers  (PLOT  relative  to  map) 

tangential  error  in  kilometers  (PLOT  relative  to  map) 

range  error  (PPI  relative  to  map) 

tangential  error  (PPI  relative  to  map) 

total  error,  in  kilometers  (PLOT  relative  to  map) 

total  error  (PPI  relative  to  map) 


AUGUST  7,  1974  "PLOT"  vs.  MAP  ANALYSES 


PAD 

1  .D. 

MAP 

rm  am 

PLOT 

rPLOT  a PLOT 

PLOT  ERROR 

XrL0T 

t— 

o 

_i 

Q-l — 

X 

PLOT 

aT0TAL 

B51 

115.5 

32.5 

115 

33 

0.5 

1  .0 

1  .  1 

C42 

89.6 

18.2 

89 

18 

0.6 

0.3 

0.7 

C48 

112.8 

36.6 

1  1  1 

360 

1.8 

1.2 

2.2 

D33 

82.5 

0.5 

82 

360 

0.5 

0.7 

0.9 

D45 

88.0 

32.5 

87 

33 

0.5 

1.0 

1  .  1 

D49 

94.6 

35.4 

93 

35 

1.6 

0.6 

1.7 

F 1 8 

70.6 

327.5 

70 

328 

0.6 

0.6 

0.8 

F39 

66.5 

18.8 

66 

19 

0.5 

0.2 

0.8 

G27 

53.9 

350.3 

53 

350 

0.9 

0.3 

0.9 

G39 

56.3 

24.3 

56 

23 

0.3 

1  .3 

1.3 

G43 

52.3 

45.2 

52 

44 

0.3 

1  .  1 

1  .  1 

G45 

68.4 

40.5 

69 

39 

1.0 

1.8 

2.1 

H37 

42.5 

16.4 

42 

14 

0.5 

1.8 

1  .9 

154 

80.0 

68.6 

79 

68 

1.0 

0.8 

1  .3 

J  6 

83.7 

284.7 

83 

285 

0.7 

0.4 

0.8 

M10 

78.3 
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AUGUST  18,  1974:  "PLOT"  vs.  MAP  ANALYSES 
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Appendix  F 


LAND  LOCATION  NOTATION 

Most  inhabited  areas  of  Alberta  are  parcelled  out  accord¬ 
ing  to  the  Dominion  of  Canada  Survey.  The  United  States  border  marks 
the  southern  limit  of  the  survey.  The  largest  divisions  in  the  Survey 
are  delineated  by  meridians,  which  are  longitude  lines  spaced  at 
four-degree  intervals.  The  eastern  border  of  Alberta  coincides  with 
the  fourth  meridian,  while  the  fifth  meridian  lies  about  13  kilometers 
west  of  Red  Deer,  Alberta.  The  Survey  divides  the  land  further  into 
townships,  sections  and  quarter-sections.  Smaller  divisions  exist, 
but  they  are  not  important  here.  The  location  of  a  township,  which 
normally  covers  36  square  miles,  is  identified  by  the  township  num¬ 
ber,  which  increases  northward  from  1  at  the  U.S.  border,  and  the 
range  number,  which  increases  westward  from  the  meridian,  beginning 
with  Range  1.  Because  of  the  geometry  of  the  survey  system,  some  town¬ 
ships  immediately  to  the  east  of  the  meridian  line  will  not  cover  a 
full  36  square  miles.  With  these  exceptions,  each  township  contains 
36  sections,  each  covering  1  square  mile,  while  each  section  is  di¬ 
vided  into  4  quarter-sect i ons .  The  following  diagram  shows  the 
labelling.  Referring  to  the  diagram,  a  typical  land  location  might 
be  NW1 5, 39, 2, 5  or  NW1 5,T39,R2, W5.  This  is  read  as  "Northwest  quarter, 
section  15,  township  39,  range  2,  west  of  the  Fifth  Meridian."  Sec¬ 
tions  or  q ua rter- sect i ons  are  sometimes  referred  to  exp  I  icitly  in  the 

text . 
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Figure  A-1-:  Land  locations.  The  township  shown  in  detail 

is  Township  39,  Range  2.  The  Fifth  Meridian 
is  at  the  eastern  border  of  Township  39, 

Range  1.  The  heavily  outlined  quarter  -  section 
is  described  as  "Northwest  Quarter,  Section  15, 
Township  39,  Range  2,  west  of  the  Fifth  Meri¬ 
dian."  In  the  text,  this  would  be  written  as 
"NW1 5 ,39,2,5". 
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